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Chapter 1 
Introduction 



1.1 History 

1.1.1 J/ifj and the November revolution 

Simultaneously discovered in November 1974 by Ting et al. [1] at the Brookhaven National 
Laboratory^ and by Richter et al. [2] at the Stanford Linear Accelerator^, the J /if) directly 
stands out from the crowd of other particles by its peculiar name composed of two symbols: 
J by Ting in comparison^ to the electric current i^{x), and '0 because of the typical pattern 
of its decay in tt+tt" at SPEAR (see Fig. 1.1). 




Figure 1.1: Left: the sound [dTng] in the Chinese alphabet; Right: typical Mark I (SPEAR) event 
display for J/ip decay (from [3]). 

Ting's experiment was based on the high-intensity proton beams of the Alternating Gra- 
dient Synchrotron (AGS) working at the energy of 30 GeV, which bombarded a fixed target 
with the consequence of producing showers of particles detectable by the appropriate ap- 
paratus. More precisely, they used a beryllium target and a beam with a momentum in 



1 BNL 

2 SLAC 

^ See however the left part of Fig. 1.1 
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the laboratory frame of 19-20 GeV. Taking advantage of a double-arm spectrometer, they 
detected a strong peak in the energy distribution of the produced electron-positron pairs at 
a value of 3.1 GeV as can be seen on Fig. 1.2. On the other hand, Richter's experiment 
used the electron-positron storage ring SPEAR, whose center-of-momentum energy could be 
tuned at the desired value. With the Mark I detector, they discovered a sharp enhancement 
of the production cross section in different channels: e^e~, H^H', 7r'''7r~,. . . The resulting 
cross sections are shown in Fig. 1.3. 




2.50 2.75 3.00 3.2i 3.50 

total eneigr IGeV] 



Figure 1.2: Energy distribution of the produced electron-positron pairs for Ting's experiment 
(from [4]). 

In the following weeks, the Frascati group (Bacci et al. [5]) confirmed the presence of this 
new particle whose mass was approximately 3.1 GeV. The confirmation was so fast that it 
was actually published in the same issue of Physical Review Letters, ie. vol. 33, no. 23, 
issued the second of December 1974. In the meantime -ten days later according Richter's 
Nobel Lecture-, his group discovered another resonant state with a slightly higher mass 
which was called"^ ip' . 

In 1976, Ting and Richter were awarded the Nobel prize for this simultaneous discovery. 
Indeed, considering the state of comprehension of particle physics at that time, this discovery 
was a shock for many. During the two years between the discovery and the award of the 



^ We shall also make us of the name 'ip{2S). 
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Nobel prize, an intense research activity in this domain took place. Were these new states 
the sign of a new generation of quarks ? Was the quark model to be finally trusted? 

An important fact to bear in mind here is that further studies based on the on- and 
off-peak production of muons and on forward/backward asymmetry in the observed angular 
distribution promptly established that the quantum numbers of the J/ip were the same as 
those of the photon , i.e. 1 . In a sense, this is not surprising at all since it is much more 
likely to produce a resonance with the same quantum number as the photon in e^e~ colliding 
experiments than other resonances. 



T 




I -U^J \ \ \ I 1 L 

200 " 'cos 91 - 0-6 



3 100 r 
S - 

b : 




20 Lj_^^j L L \ I \ \ I J U 

3.050 3.090 3.100 3.110 3.120 3.130 
ENERGY Ec.m.(GeV) 

Figure 1.3: Various pair production cross sections as functions of the center-of-momentum energy 
(top: hadrons; middle: tt+tt", yU"''yU~ and K^K~; bottom: e^e~) [2]. 

Moreover, since the following ratio 

^ cross section for e"^e~ — > hadrons 
cross section for e^e~ — > 

was much larger on-resonance than off, it was then clear that the J/ip did have direct 
hadronic decays. The same conclusion held for the ip' as well. The study of multiplicity 
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in pion decays indicated that decays were restricted by a specific selection rule, called 
G-parity conservation, known to hold only for hadrons. Consequently, J/tp and ip' entered 
the family of hadrons with isospin and G-parity -1. 

Particles with charge conjugation C different from -1, i.e. with quantum numbers differ- 
ent from those of the photon, were found later. Indeed, they were only produced by decay 
of ip" and the detection of the radiated photon during the (electromagnetic) decay was then 
required. The first to achieve this task and discover a new state was the DASP collabora- 
tion [6] based at DESY (Deutsches Elektronen-Synchrotron), Hamburg, working at an e^e~ 
storage ring called DORIS. This new particle, named^ Pc, had a mass of approximately 3.500 
GeV. At SPEAR other resonances at 3.415, 3.450 and 3.550 GeV were discovered, the 3.500 
GeV state was confirmed. Later, these states were shown to be C = +1. 

Coming back to the ratio R, let us now analyse the implication in the framework of the 
quark model postulated by Gell-Mann^ and Zweig in 1963. In 1974 at the London conference, 
Richter presented [7] the experimental situation as in Fig. 1.4. 




Figure 1.4: Experimental status of R as of July 1974 (from [7]). 



In the framework of the Gell-Mann-Zweig quark model with three quarks, a plateau was 

^ The name chosen here reveals that physicists already had at this time an interpretation of this state as 
bound-state of quarks c. Indeed, the symbol P follows from the classical spectroscopic notation for a £ = 1 
state, where i stands for the angular momentum quantum number. 

^ The name "quark" was chosen by Gell-Mann. He found this word in the James Joyce's novel Finnegans 
Wake where the following enigmatic sentence appears: 'Three quarks for Muster Mark'. This name comes 
from the standard English verb quark, meaning "to caw, to croak" , and also from the dialectal verb quawk, 
meaning "to caw, to screech like a bird" . 
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expected with a value of 2/3 or 2 if the quark were considered as "coloured" - a new concept 
introduced recently then -. The sign of a new quark - as the charm quark first proposed by 
Bjorken and Glashow in 1964 [8]- would have been another plateau from a certain energy 
(roughly its mass) with a height depending on its charge. Retrospectively, one cannot blame 
anyone for not having seen this "plateau" on this kind of plots. Quoting Richter, "the 
situation that prevailed in the Summer of 1974" was "vast confusion" . 

This charm quark was also theoretically needed following the work of Glashow, Iliopoulos 
and Maiani [9], which suggested that a symmetry between leptons and quarks should exist 
in order to cancel the anomaly in weak decays. And following this symmetry proposal, the 
charm quark was expected to exist and to have an electric charge 2/3. R was therefore to 
be 10/3 in the coloured quark model, still not obvious in Fig. 1.4. This explains why the 
discovery of such sharp and easily interpreted resonances in November 1974 was a revolution 
in particle physics. 

It became quite rapidly obvious that the J/ip was the lowest-mass cc system with the 
same quantum numbers as photons - explaining why it was so much produced compared 
to some other members of its system. These cc bound states were named "charmonium" , 
firstly by Appelquist, De Riijula, Politzer and Glashow [10]^, in analogy with positronium, 
which has a similar bound-state level structure. 

At that time, the charm had nevertheless always been found hidden, that is in charm-anti- 
charm bound states. On the other hand, it became quickly clear that the mass of the lightest 
possible state constituted of a charm quark and a light quark, thus an explicit charmed 
meson, named was to lie around 1900 MeV. A more precise lower bound, 1843 MeV, was 
inferred from the fact that the if)' {m^i = 3684 MeV) was too narrow or, in other words, had 
too long a lifetime to be above the threshold to produce a pair of -D's, that is 2m/). The 
upper bound was, on the other hand, deduced from the point where R started to rise before 
reaching the plateau with 4 active quark fiavours. 

In order to study these D mesons, the investigations were based on the assumption that 
the D was to decay weakly and preferentially into strange quarks. The weak character of 
the decay motivated physicists to search for parity-violation signals. The first resonance 
attributed to meson was found in K~^7t~^ decay by Goldhaber et al. [11] in 1976. A 
little later, and D~ were also discovered as well as an excited state, D* , with a mass 
compatible with the decay ip'" — D^D*. And, finally, the most conclusive evidence for the 
discovery of charmed meson was the observation of parity violation in D decays [12]. To 
complete the picture of the charm family, the first charmed baryon was discovered during the 
same year [13]. The quarks were not anymore just a way to interpret symmetry in masses 
and spins of particles, they had acquired an existence. 



^ To be complete, the initial title of the preprint "Charmonium spectroscopy" was immediately refused by 
the senior editor Pasternak of Physical Review Letters because of "frivolous new terminology" . A consensus 
eventually emerged thanks to Glashow: the term was accepted in the text but not in the title! 
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1.1.2 The truth: a third generation 

In 1975, another brand new particle was unmasked at SLAC by Perl et al. [14], for which 
he was awarded the Nobel prize twenty years later. This was the first particle of a third 
generation of quarks and leptons, the r, a leptor?. Following the standard model, which was 
more and more a reference trusted by most, two other quarks were expected to exist. Their 
discovery would then be the very proof of the theory that was developed since the sixties. 
Two names for the fifth quark were already chosen: "beauty" and "bottom", and in both 
cases represented by the letter 6; the sixth quark was as well already christened with the 
letter t for the suggested names^ "true" or "top" . 

The wait was not long: two years. After a false discovery^° of a resonance at 6.0 GeV, a 
new dimuon resonance similar to J /ip and called T was brought to light at Fermilab, thanks 
to the FNAL proton synchrotron accelerator, by Herb et al. [16], with a mass of 9.0 GeV; 
as for charmonia, the first radial excited state (T(2S')) was directly found [17] thereafter. 
Various confirmations of these discoveries were not long to come. The SS* state was then 
found [18] at Fermilab as well as evidence that the 45" state was lying above the threshold 
for the production of B mesons. The latter was confirmed at the Cornell e^e~ storage ring 
with the CLEO detector. The first evidence for B meson and, thus, for unhidden h quark, 
was also brought by the CLEO collaboration [19] in 1980. One year later, "the first evidence 
for baryons with naked beauty" (szc) was reported by CERN physicists [20]"*^^. 

Another decade was needed for the discovery of the sixth quark which was definitely 
christened "top". Indeed, in 1994, the CDF Collaboration observed it at the Tevatron 
collider in Fermilab [21]. This accelerator - which reaches now a cm. energy of 1960 GeV 
- and this collaboration will accompany us throughout this thesis, but for a different signal. 
Unfortunately, due to its very short lifetime, this quark cannot bind with its antiquark to 
form the toponium. To conclude this historical prelude, we give the spectra of the cc and 
hh systems as well as two tables summing up the characteristics of the observed states as of 
today. 



Meson 






Mass (GeV) 


r,.M (keV) 




1 


0^+ 


2.980 


N/A 




1 


1~ 


3.097 


5.40 


XcO, Xcl, Xc2 


1 ^-Po,l,2 


0++,l++,2++ 


3.415,3.511,3.556 


N/A 


he 




1+- 


3.523 


N/A 


Vc{2S) 


2^50 


0-+ 


3.594 


N/A 


^' 


2 3^1 


1-- 


3.686 


2.12 



Table 1.1: Properties of charmonia (cf. [57]). 



^ in contradiction with the etymology because it was heavier than most of known baryons at that time. 
^ The name "truth" was given as a reference to the fact that its discovery should be the proof of the 
truth of the standard model. 

Nonetheless, this paper [15] first suggested the notation T for any "onset of high-mass dilepton physics" . 
Here is an unambiguous sign that European editors had a much softer editorial approach concerning 
titles of articles. . . 
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¥(2S) 



Xc2(1 P) 



hadrons 




hadrons hadrons radiative 
jPC = 0-+ 1— 0++ 



1++ 



2++ 



Figure 1.5: Spectrum and transitions of the charmonium family (from [57]). 



T (11 020) 



T (10860) 



T(4S) 



BB threshold 



ilf,(3S) 



hadrons 



i1b(2S) 



ilb(1S) 



T(1S) 




jPC = 0-+ 1— 1+- 0++ 1++ 2++ 

Figure 1.6: Spectrum and transitions of the bottomonium family (from [57]). 
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Meson 






Mass (GeV) 


r^;. (keV) 


T(l^) 




1— 


9.460 


1.26 


Xbo, Xbi , Xb2{^P) 




0++,l++,2++ 


9.860,9.893,9.913 


N/A 


T{2S) 


2 


1 — 


10.023 


0.32 


Xbo, Xbi, Xb2 (2P) 


1 ^-Po,l,2 


0++,l++,2++ 


10.232,10.255,10.269 


N/A 


T(3^) 


3351 


1— 


10.355 


0.48 



Table 1.2: Properties of bottomonia (cf. [57]) . 



1.2 Motivations 

1.2.1 Why study quarkonium physics in the high-energy regime ? 

The fact that J/ip was the first particle with a c quark to be discovered is not an accident: 
it is in fact of fundamental importance in particle physics. Let us go back to its discovery: 
as we have said, two totally different experiments were involved in it. The first was a fixed- 
target experiment, whose main advantage was to scan all the energy range below the cm. 
energy of the collision. Its main disadvantage was to produce many other particles than the 
resonance which it looked for. The second experiment was taking place at an e~^e~ collider 
with totally opposite characteristics, one given energy and a clean environment surrounding 
the expected resonance. In a sense, the J/ip discovery is as if the Higgs boson had been 
discovered simultaneously at LEP and at the Tevatron. . . 

The discovery in a fixed-target experiment shows that the J/ip leptonic decay is signifi- 
cant, making it highly detectable compared to particles only detected through their hadronic 
decays. J/ip is in fact so highly detectable that it was not a major problem to have bunches 
of other hadrons produced during the same collision. 

Let us explain why this leptonic decay is significant. Firstly, the open charm channel 
(decay to D-mesons) is kinematically blocked since rrij/^ < 2m/); secondly, colour prevents 
decays in a single gluon; thirdly, the Landau- Yang theorem [22] or the G-parity conservation 
forbids its decay into two gluons. Therefore, three gluon emissions are required and, as the 
gluon emission rate from heavy quarks is reduced^^, one is not too surprised to get a total 
width of 91 ± 3 keV, much smaller than that of usual hadronic decays. Moreover, this value 
also includes leptonic channels. On the other hand, leptonic decays are not subject to all 
these restrictions: a single photon decay is allowed, making them competitive with hadronic 
decays. More precisely, they account for one third of the total width. 

This significant branching ratio into leptons, especially into muons, is an incredible 
chance, like a metal detector to find a needle in a haystack. In high-energy hadron col- 
liders, the number of background hadrons becomes annoyingly large and complicates the 
detection of particles with constrained kinematical characteristics. However, for muons, this 
is a completely different story, as they can go through the calorimeters without interaction 
and be seen in dedicated detectors. These can measure precisely their charge, momentum 

■^'^ This comes from the fact that the scale of the process is of the order of mg and makes as relatively 
low. In fact, this results in the Zweig rules. 
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and direction, determining therefore their track. It is in turn possible to discriminate pre- 
cisely which dimuon, out of the complete set of the detected muons and anti-muons, results 
from the decay of a J/ip. The latter being a significant source of dimuons, it is therefore 
an ideal probe to understand high-energy hadron collisions where many other studies are 
prevented due to a high background. 

Concerning its discovery with an e^e~ collider, it results partly from the fact that the 
rather low J/ip mass was accessible to the available experiments whereas had it been lower, 
say 1.5-2 GeV, J/ip would perhaps have been lost among excited states of strange baryons. 

The reasonably low mass of the J/ ijj makes it highly producible, with an easily reachable 
threshold in energy. This explains why it is abundantly produced at the Tevatron, even at 
relatively large transverse momentum. This qualifies it, once more, as a good probe to study 
mechanisms associated with its production. Furthermore, its production at large transverse 
momentum, i.e. in a range where perturbative Quantum Chromodynamics (pQCD) is likely 
to be applicable, offers us a way to test the application of QCD. This statement that we 
have a sufficiently high scale to use pQCD will follow us all along this thesis. 

To complete the picture, the T's, which are heavier, will be a useful benchmark; even 
though they are less produced due to their higher mass, they set a higher scale, for which 
pQCD should apply better, even under unfavourable kinematical conditions. Their quarks 
being heavier, their velocity inside the meson is expected to be much smaller and, conse- 
quently, relativistic effects to be less important. This lowers theoretical uncertainties and 
constrains more the quantities to be tested. 

To be more definite, quarkonium production was thought, till the middle of the nineties, 
to be a good way to extract the gluon distribution function, both in electron-proton collisions 
where the modelisations are more reliable, and in proton-proton collisions where two gluon 
pdf 's enter the calculation. As we shall see, the mechanisms of production of these quarkonia 
are not yet understood and do not provide us with a powerful tool to achieve the extraction 
of these pdf 's. 

1.2.2 The eighties: the time for predictions 
The Colour-Singlet Model 

This model is^'^ the most natural application of QCD to heavy quarkonium production in the 
high-energy regime. It takes its inspiration in the factorisation theorem of QCD [23]^^ where 
the hard part is calculated by the strict application pQCD and the soft part is factorised in 
a universal wave function. This model is meant to describe the production not only of J/ip, 
ip{2S), T(15'), T(25') and T{3S), i.e. the ^5*1 states, but also the singlet S states rjc and r]b 
as well as the P (x) and D states. 

Its greatest quality resides in its predictive power as the only input, apart form the pdf, 
namely the wave function, can be determined from data on decay processes or by application 
of potential models. Nothing more is required. It was first applied to hadron colliders [24], 



before the inclusion of fragmentation contributions, 
proven for some definite cases, e.g. Drell-Yan process. 
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then to electron-proton colliders [25]. The cross sections for '^5'i states were then calculated, 
as well as also for t] and for charmonium and for bottomonium. These calculations were 
compared to ISR and FNAL data from ^/s = 27 GeV to ^/s = 63 GeV for which the data 
extended to 6 GeV for the transverse momentum. Updates [26, 27] of the model to describe 
collisions at the CERN pp collider {^/s = 630 GeV) were then presented. At that energy, 
the possibility that the charmonium be produced from the decay of a beauty hadron was 
becoming competitive. Predictions for Tevatron energies were also made [27]. 

In order to introduce the reader to several concepts and quantities that will be useful 
throughout this work, let us proceed with a detailed description of this model. 

The Model as of early 90's 

It was then based on several approximations or postulates: 

• If we decompose the quarkonium production in two steps, first the creation of two on- 
shell heavy quarks {Q and Q) and then their binding to make the meson, one postulates 
the factorisation of these two processes. 

• As the scale of the first process is approximately +p'ti one considers it as a pertur- 
bative one. One supposes that its cross section be computable with Feynman- diagram 



• As we consider only bound states of heavy quarks (charm and bottom quarks), their 
velocity in the meson must be small. One therefore supposes that the meson be created 
with its 2 constituent quarks at rest in the meson frame. This is the static approxima- 
tion. 

• One finally assumes that the colour and the spin of the QQ pair do not change during 
the binding. Besides, as physical states are colourless, one requires the pair be produced 
in a colour- singlet state. This explains the name Colour-Singlet Model (GSM). 

In high-energy hadronic collisions, the leading contribution comes from a gluon fusion 
process; as the energy of the collider increases, the initial parton momentum fraction 
needed to produce the quarkonium decreases to reach the region in x where the number of 
gluons becomes much larger than the number of quarks. One has then only six Feynman 
diagrams for the ^5*1 states production associated with a gluon^^ (see Fig. 1.7). 

One usually starts with A4{p), the perturbative amphtude to produce the heavy quark 
pair on-shell with relative momentum p and in a configuration similar to the one of the meson. 
To realize the latter constraint, one introduces a projection operator^^; the amplitude J\4{p) 
is then simply calculated with the usual Feynman rules. 

The amplitude to produce the meson is thence given by 



methods. 




(1.2) 



^ This is the dominant process when the transverse momentum of the meson is non- vanishing. 
^ In fact, this amounts to associate a 7^ matrix to pseudoscalars, 7^ to vectors, etc. 



1.2 Motivations 



23 



35 



?1 




35 




Figure 1.7: The 6 diagrams for gg 



^Sig at LO within the CSM. 



where is the usual Schrodinger wave- function. 

Fortunately, one does not have to carry out the integration thanks to the the static 
approximation which amounts to considering the first non-vanishing term of A when the 
perturbative part Ai is expanded in p. For S'-wave, this gives 



where is the wave-function in coordinate space, and ^l^l^^g ^(0) value at the origin. 
For P- waves, ^1/(0) is zero, and the second term in the Taylor serie must be considered; this 
makes appear ^^'(0). ^E'(O) (or ^I/'(0)) is the non-perturbative input, which is also present in 
the leptonic decay width from which it can be extracted. 

If the perturbative part Ai{p) is calculated at the leading order in a^, this will be referred 
to as the Leading Order CSM (LO CSM) [24, 25]. 

The Colour Evaporation Model 

This model was initially introduced in 1977 [29, 30] and was revived in 1996 by Halzen et 
al. [31], following a prescription by Buchmiiller and Hebecker [32] to explain rapidity gap 
production in deep-inelastic scattering. 

The main characteristic of this model is the total rejection of the role of colour. Con- 
trarily to the CSM, the heavy quark pair produced by the perturbative interaction is not 
assumed to be in a colour-singlet state. One simply considers that the colour and the spin 
of the asymptotic QQ state is randomised by numerous soft interactions occurring after its 
production, and that, as a consequence, it is not correlated with the quantum numbers of 
the pair right after its production. 




(1.3) 
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A first consequence of this statement is that the production of a ^5*1 state by one gluon 
is possible, whereas in the CSM it was forbidden solely by colour conservation. In addition, 
the probability that the QQ pair eventually be in a colour-singlet state is therefore ^, which 
gives the total cross section to produce a quarkonium: 

CTonium = - / dm— . (1.4) 

9 J2mQ dm 

This amounts to integrating the cross sections of production of QQ from the threshold 2mQ 
up to the threshold to produce two charm or beauty mesons (designated here by qQ) and 
one divides by 9 to get the probability of having a colour- singlet state. 

The procedure to get the cross section for a definite state, for instance a J/ip, is tanta- 
mount to "distributing" the cross sections among all states: 

Cj/i/j = Pj/tpCTonium- (1-5) 

The natural value for pj/^, as well as for the other states in that approximation, is the 
inverse of the number of quarkonium lying between 2mc and 2m£). This can be refined by 
factors arising from spin multiplicity arguments or by the consideration of the mass difference 
between the produced and the final states. These are included in the Soft-Colour- Interactions 
approach (SCI) [33]. 

By construction, this model is unable to give information about the polarisation of the 
quarkonium produced, which is a key test for the other models [34]. Furthermore, nothing 
but fits can determine the values to input for p, except perhaps the SCI approach which relies 
on Monte-Carlo simulation, but this does not help to understand the underlying mechanisms. 
Considering production ratios for charmonium states within the simplest approach for spin, 
we should have for instance cr[ric] : a[J/ip] = 1:3 and cr[Xco] '■ <^[Xci] '■ cr[Xc2] = 1:3:5, 
whereas deviations from the predicted ratio for Xci and Xc2 have been observed. Moreover, it 
is unable to describe the observed variation - from one process to another - of the production 
ratios for charmonium states. For example, the ratio of the cross sections for Xc and J/ip 
differs significantly in photoproduction and hadroproduction, whereas for the CEM these 
number are strictly constant. 

All these arguments make us think that despite its simplicity and its phenomenological 
reasonable grounds, this model is less reliable than the CSM. It is also instructive to point out 
that the invocation of reinteractions after the QQ pair production contradicts factorisation, 
which is albeit required when the pdf are used. However, as we shall see now, the CSM has 
undergone in the nineties a lashing denial from the data. 

1.2.3 The nineties: the time for disillusions 

ip' anomaly 

In 1984, Halzen et al. [26] noticed that charmonium production from a weak B decay could 
be significant at high energy - their prediction were made for ^/s = 540 GeV - and could 



1.2 Motivations 



25 



even dominate at high enough pt- This can be easily understood having in mind that the B 
meson is produced by fragmentation of a 6 quark - the latter dresses with a light quark -. 
And to produce only one b quark at high px is not so burdensome; the price to pay is only 
to put one quark propagator off-shell, instead of two for gg — > ipg. 

This idea was confirmed by the calculations of Glover et al. [27], which were used as a 
benchmark by the UAl Collaboration [28]. After the introduction oi a. K factor of 2, the 
measurements agreed with the predictions; however the pt slope was not compatible with 
6-quark fragmentation simulations. 

From 1993, the CDF collaboration undertook an analysis of tp'^'^ production. They man- 
aged to extract unambiguously the prompt component of the signal (not from B decay) using 
a Sihcon Vertex Detector (SVX) [35] 

The preliminary results showed an unexpectedly large prompt component. For the if)' , the 
prompt cross section was orders of magnitude above the predictions of the LO CSM (compare 
the data to the dashed curve on Fig. 1.8 (b)) . This problem was then referred to as the ip' 
anomaly. For J/ip, the discrepancy was smaller (Fig. 1.8 (a)), but it was conceivably blurred 
by the fact that a significant part of the production was believed to come from Xc radiative 
feed-down, but no experimental results were there to confirm this hypothesis. 




(a) J/V (b) V'' 

Figure 1.8: Cross sections obtained by Braaten et al. for LO CSM (dashed curves) of J/^ (a) 
and ip' (b), as well as for CSM fragmentation contribution (solid curves). In each 
case, the two curves depict the extremum values obtained by varying parameters such 
as rrLc and the different scales: hr, hf, l^frag (from [40]). 



In the following, ip stands for both J/tp and ip' . 

The details of the analysis are given in the following chapter. 



26 



Introduction 



Fragmentation within the CSM 

The prediction from the LO CSM for prompt tp being significantly below these preliminary 
measurements by CDF, Braaten and Yuan [36] pointed out in 1993 that gluon fragmentation 
processes, even though of higher order in as, were to prevail over the LO CSM for 5- wave 
mesons at large pr, in same spirit as the production from B was dominant as it came from a 
fragmentation process. Following this paper, with Cheung and Fleming [37], they considered 
the fragmentation of c quark into a. ip in Zq decay. From this calculation, they extracted 
the corresponding fragmentation function. In another paper [38], they considered gluon 
fragmentation into P-wave mesons. All the tools were then at hand for a full prediction of 
the prompt component of the J/ ip and ip' at the Tevatron. This was realised simultaneously 
by Cacciari and Greco [39] and by Braaten et al. [40]. We shall discuss their results after a 
digression about the formalism. 

To understand how the fragmentation process is factorised from the hard interaction and 
how it can be so enhanced that it becomes larger than the LO terms, let us consider the rjc 
production as was done in [36]. Since its charge conjugation is +1, it can be produced in 
association with two gluons. The LO CSM diagram is depicted in Fig. 1.9 (a), whereas one 
of the NLO CSM diagram is on depicted in Fig. 1.9 (b). 




"OooooO" 






(a) J/^ 



(b) V(25) 



Figure 1.9: The LO CSM diagram for rjc production in gluon fusion process (a), one of the NLO 
CSM diagram for the same process (b) . 



Normally, the two virtual gluons of the NLO process (Fig. 1.9 (b)) are off their mass shell 
by ~ Pt, as well as the virtual gluon and the virtual quark in the LO process (Fig. 1.9 (a)). 
The NLO diagram is then solely suppressed by one power of as at the chosen scale, which can 
be pt- However, in some regions of phase space, namely when the off-shell gluon attached 
to the cc is soft - with an off-shellness of the order of rric -, this diagram is enhanced by a 

2 

factor of ^ compared to the LO one. At sufficiently large values oipx, this enhancement can 
overcome the as suppression. This behaviour can in fact be traced back to the impossibility 
of large relative momentum configurations for the constituents of the bound state. 
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In this region, this NLO contribution can factorised: the first factor is the amplitude 
99 ~^ 99* where the off-shell and the on-shell gluons in the final state have a high pt and 
where the same off-shell gluon has a small off-shellness <^ p\.\ the second factor is the 
propagator of the gluon, and the third corresponds to the fragmentation function 

g* rjcQ. Taking into consideration the restriction on the relative momentum within the 
bound state, we can therefore express the differential cross section to produce the rjc with a 
momentum P as: 

da^XP) - I dzda, Dg^.^z), (1.6) 

dag (j) being the differential cross section for gg — gg* where the momentum of g* is j 
and Dg^rici^) is the gluon fragmentation function. 

It is now possible to generalise this to all quarkonia. To take into account coUinear 
production of fragmenting gluons from a higher energy parton, a fragmentation scale Hfrag 
is introduced which is meant to include logarithms of pt/^q appearing in these collinear 
emissions (ln(pT/'^Q) = ^^{pt/ fJ'frag) + ^^{fJ' frag / mq)) . Hcuce, to all orders in as, we have 
the following fragmentation cross section for a quarkonium Q: 

^ /■! p 

^QiP)-y2 dzdaii — , fx frag) D^^q{z,^ frag). (1.7) 

The scale is as usual chosen to avoid large logarithms of px/fJ'frag in ai{ — ,fifrag), that is 
fifrag — Pt- The summatiou of the corresponding large logarithms of jJifrag/'^nQ appearing in 
the fragmentation function is realised via an evolution equation [41, 42]. 




Figure 1.10: Illustration of how to obtain D{z,mc). 

The interesting point raised by Braaten and Yuan is that the fragmentation functions can 
be calculated perturbatively in as at the scale fifrag = 2mQ. For instance, in the case of gluon 
fragmentation into an r]c, the trick is to note that dzD^^z, rric) is the ratio to the rates for 
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the well-known gg gg process and gg ^cQQ (Fig. 1.10). After some manipulations, the 
fragmentation function can be obtained from this ratio by identifying the integrand in the z 
integral. This gives : 

D,^,^(z,2m,) = ia2(2^^) KWp^ _ ^^2 ^ 3(1 _ ^)in(l - z)]. (1.8) 

The other fragmentation functions of a given parton i into a given quarkonium Q, 
Di^Q{z, iJ,frag), Were obtained in the same spirit. For the Tevatron, the differential cross 
section versus px of various CSM fragmentation processes are plotted in Fig. 1.11. 

10^ 
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Figure 1.11: Differential cross section versus of various CSM fragmentation processes for J/ip 
to be compared with the LO contributions (from [40]). 

The prompt component of the J/ ip and the direct component of the tp' could in turn be 
obtained and compared with the preliminary data of CDF (see the solid curves in plots in 
Fig. 1.8 above). For the J/ip, the previous disagreement was reduced and could be accounted 
for by the theoretical and experimental uncertainties; on the other hand, for the ijj' , the 
disagreement continued to be dramatic. The situation would be clarified by the extraction 
of the direct component for J/ijj, for which theoretical uncertainties are reduced and are 
similar to those for the ip\ 

The CDF collaboration undertook the disentanglement of the direct J/ ip signal [43] . They 
searched for J/ ip associated with the photon emitted during this radiative decay: the result 
was a direct cross section 30 times above expectations from LO CSM plus fragmentation. 
This was the confirmation that the CSM was not the suitable model for heavy quarkonium 
production in high-energy hadron collisions. 
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Figure 1.12: Differential cross section versus of the CSM (fragmentation and LO) production 
to be compared with the direct production of J/^p from CDF (Adapted from [44]). 



NRQCD: successes and doubts. . . 

In 1992, Bodwin et al. considered new Fock state contributions in order to cancel exactly the 
IR divergence in the light hadron decays of Xci at LO. This decay proceeds via two gluons, 
one real and one off-shell, which gives the light-quark pair forming the hadron. When the 
first gluon becomes soft, the decay width diverges. The conventional treatment [45], which 
amounts to regulating the divergence by an infrared cut-off identified with the binding energy 
of the bound state, was not satisfactory: it supposed a logarithmic dependence of ^'(0) upon 
the binding energy. They looked at this divergence as being a clear sign that the factorisation 
was violated in the CSM. 

Their new Fock states for Xc were e.g. a gluon plus a cc pair, in a '^5'i configuration 
and in a colour-octet state. The decay of this Fock state occurred through the transition 
of the coloured pair into a gluon (plus the other gluon already present in the Fock state as 
a spectator). This involves a new phenomenological parameter, ifg, which is related to the 
probability that the cc pair of the Xc be in a colour-octet S'-wave state. The key point of 
this procedure is that a logarithmic dependence on a new scale A - a typical momentum 
scale for the light quark - appears naturally within the effective field theory Non-Relativistic 
Quantum Chromodynamic (NRQCD) [46, 47]. 

This effective theory is based on a systematic expansion in both and which is the 
quark velocity within the bound state. For charmonium, v1 ~ 0.23 and for bottomonium 
vl ~ 0.08. One of the main novel features of this theory is the introduction of dynamical 
gluons in the Fock state decomposition of the physical quarkonium states. In the case of 
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5'-wave orthoquarkonia (^^i), we schematically have, in the Coulomb gauge: 

\Qq) =0{l)\QQ[ 'S^^\) + 0{v)\QQ[ ^Pfg]) + 0{v')\QQ[ 'sf^g]) 
+ 0{v^)\QQ[ 'S^^^'^gg]) + 0{v')\QQ[ ''of'^gg]) + ... 

whereas for P-wave orthoquarkonia {^Pj), the decomposition is as follows 

\Qq,) =Oil)\QQ[ ^'Pf ]) + Oiv)\QQ[ 'S['^g]) + ... (1.10) 

In these two formulae, the superscripts (1) and (8) indicate the colour state of the QQ pair. 
The (9(f") factors indicate the order in the velocity expansion at which the corresponding 
Fock state participates to the creation or the annihilation of quarkonia. This follows from 
the velocity scaling rules of NRQCD (see e.g. [47]). 

In this formalism, it is thus possible to demonstrate, in the limit of large quark mass, 
the factorisation between the short- distance - and perturbative - contributions and the 
hadronisation of the QQ, described by non-perturbative matrix elements defined within 
NRQCD. For instance, the differential cross section for the production of a quarkonium Q 
associated with some other hadron X reads 

da{Q + X) = J2dHQQ[''^'L^/'^]+X){OQ[''^'L^j''^]), (1.11) 

where the sum stands for S, L, J and the colour. 

The long-distance matrix element (LDME) {0^[^^~^^ L^j'^^]) takes account of the transi- 
tion between the QQ pair and the final physical state Q. Its power scaling rule comes both 
from the suppression in v of the Fock state component [^'^+^Lj in the wave function of Q 
and from the scaling of the NRQCD interaction responsible for the transition. 

Usually, one defines 0^[^^~^^L^j'^'^] as the production operator that creates and annihilates 
a point-like QQ pair in the specified state. This has the following general expression 



Sr2S+i, (1,8)1 = ^t^^ (^Y,\Q + x){Q + X\^ ij^Kx 
= x'K'ilj{a^QaQ)ip^Kx, 



:i.l2) 



where^^ ip and x are Pauli spinors and the matrix K is a product of colour, spin and covariant 
derivative factors. These factors can be obtained from the NRQCD lagrangian [47]. For 
instance, 

0S[ 3^(1)] _ x^^^{aQaQ)i^^ffX- However some transitions require the presence of 
chromomagnetic (AL = and AS = ±1) and chromoelectric (AL = ±1 and AS* = 0) terms 
for which the expressions of K are more complicated. 

On the other hand, the general scaling rules relative to the LDME's [44] give: 
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The second line of Eq. (1.12) is nothing but a short way of expressing this operator. 
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where E and M are the minimum number of chromoelectric and chromomagnetic transitions 
for the QQ to go from the state [^^~^^L^J'^^] to the dominant quarkonium Q Fock state. 

The idea of combining fragmentation as the main source of production with allowed 
transitions between a Xc to a ^Si in a colour-octet state, had already been applied by Braaten 
and Yuan in the above mentioned paper [38]. Indeed, similar formulae can be written for 
fragmentation functions [48]: 

D,^Q{^,f^) = 5^rfps.i4i.B)](^,/x)(02['^+^4'''^]), (1-14) 

where d[.] {z, /i) accounts for short- distance contributions and does not depend on which Q is 
involved. But since the theoretical predictions for prompt J/ip production did not disagree 
dramatically with data and since there was no possible Xc decay to a possible enhancement 
of the Xc cross section by colour-octet mechanism (COM) was not seen as a key-point both 
for J/ip and ip' production. 

However, in the case of J/ip and ip' production, COM could still matter in a different 
manner: fragmentation of a gluon into a ^Pf is possible with solely one gluon emission and 
fragmentation into a ^S[^^ requires no further gluon emission (at least in the hard process 
described by d[.]{z,fi)). Concerning the latter process, as two chromoelectric transitions are 
required for the transition \ccgg) to |cc), the associated LDME {0'^[^S[^^]) was expected to 
scale as m^w''. In fact, d.^ (s),, the contribution to the fragmentation function of the short- 

distance process g was already known since the paper of Braaten and Yuan [38] and 

could be used here, as the hard part d[.]{z,fi) of the fragmentation process is independent of 
the quarkonium. 

In a key paper, Braaten and Fleming [38] combined everything together to calculate, for 
the Tevatron, the fragmentation rate of a gluon into an octet ^Si that subsequently evolves 
into a V'- They obtained, with^o (0^^(25) [35.(8)^^ ^ ^ 2 x 10"^ GeV^ a perfect agreement 
with the CDF preliminary data: compare these predictions in Fig. 1.13 with the previous 
ones in Fig. 1.8 (b). 

A straightforward and unavoidable consequence of this solution to the ip' anomaly was 
raised by Cho and Wise [49]: the tp', produced in the latter way, was 100% transversally 
polarised. A suggested test of this prediction was also given, i.e. the measurement of the 
lepton angular distribution in tp' which should behave as 

^I_(^'^£+r) oc (l + acos^e), (1.15) 

with a=l for 100% transversally polarised particles. Further details are given in section 2.5. 

Following these studies, a complete survey on the colour-octet mechanism was made in 
two papers by Cho and Leibovich [50, 51]. The main achievements of these papers were the 
calculation of colour-octet P-state contributions to ip, the first prediction for direct colour- 
octet J/ip production through e.g. the Xc feed-down calculation, the first predictions for T 

This corresponds to a suppression of 25 compared to the "colour-singlct matrix element" 
^0V'(2S)[3^(i)]^^ which scales as m^v^. This is thus in reasonnable agreement with a suppression. 
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Figure 1.13: Differential cross section versus of the CSM (dotted: LO; dashed: fragmentation 
and LO) production and of the COM fragmentation (solid curve) to be compared 
with CDF preliminary direct production of ^/^' (from [38]). 

and the complete set of 2 3 parton processes like ij cc + k, all this in agreement with 
the data. 

A few year later, however, the CDF collaboration carried out the study of polarisation 
of prompt ijj [52]. The outcome was a = for ip', or even a < when pj- becomes large. 
For J/ip, the same behaviour showed up, whereas comparison with theory is still difficult 
because of the Xc feed-down. In any case, no sign of strong transverse polarisation was seen, 
in total contradiction with all predictions of the COM. 

Summary 

To sum up, we have seen that the LO CSM gives smaller contributions than the fragmentation 
processes for the CSM, where the heavy quark pair produced in the hard part of the process 
has strictly the same quantum numbers as the final physical state. On the other hand, 
in order to solve problems connected with the factorisation for Xc decay, one is driven to 
introduce new contributions from higher Fock states. 

In the context of production, this allows a given physical state to come from a heavy 
quark pair with different quantum numbers, e.g. in a colour-octet state. For the case of 
fragmentation, even though these higher Fock-state contributions are suppressed, they in- 
volve lower-order contributions in and are shown to dominate for intermediate and large 
momentum transfer. 

All these statements can be easily understood by analysing the large Pt behaviour of the 
differential cross section, from typical diagrams relevant for the considered processes. 
For the LO CSM, at large Pt, the two internal lines not linked with the bound state are 
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off-shell by ~ and make the partonic differential cross section behave like ~ as 
depicted in Fig. 1.14 (a). 

To what concerns fragmentation, it behaves like a single particle production process, thus 
like -jp^ ~ -p^, if we consider that, at large Pt, the effects of the quarkonium mass become 

negligible. Even though fragmentation processes (Fig. 1.14 (b)) are of higher order in a.^ (af 
instead of a^), they scale like (g^^-yj compared to LO CSM and dominate over the latter for 
Pt > 2mc. 

The inclusion of colour-octet channels can be analysed with similar arguments. For the 
^S^^ channel, the only difference comes from the powers of as and v. No similar estimation 
for power counting of the Pt dependence can be done for the other 0{v^) channels [51], 
however the gluon t-channel dominance argument [44] provides us with the scaling Indeed 
these channels are significant at moderate Pt of the order of 2mc, as one can conclude from 
the complete analysis made in [51]. 

On the other hand, the statement that the polarisation of ip produced at large Pt be 
entirely transverse - at variance with experimental data - follows simply and only from the 
on-shellness of the fragmenting gluon and from the approximate spin symmetry of NRQCD. 

Considering that no solution of the polarisation issue has been proposed so far, we have 
decided to undertake an analysis of the non-static effect of the LO CSM, with the first 
hope to state that the static approximation was not viable for previously unknown reasons, 
but also with the goal to attain new insights about some other unknown effects linked 
with quarkonium production. As we shall see, it seems that the static approximation is 
not problematic. However, we have realised than unconstrained contributions could be 
introduced in the description of the production of the quarkonium in association with a gluon 
- as it should be at non vanishing transverse momentum. This introduction can bring the 
theoretical predictions in agreement with experimental measurements, for the cross section 
and for the polarisation. 
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(a) leading-order colour-singlet: g + g ^ cc[^Si ] + g 



gsSLSLOSiSL 




{b) colour-singlet fragmentation: g + g ^ [ccpS"} '] + gg] + g 




(c) colour-octet fragmentation: g + g ^ ccpS'j ] + g 




(d) colour-octet t-channel gluon exchange: g + g ^ cc[^Sq'\^Pj^^] + g 




Figure 1.14: Summary of the different processes arising in the direct production of ip within the 
CSM and the COM and their pr dependence (from Kramer [44]). 



Chapter 2 

Experimental results 



2.1 Foreword 

Before broaching a review of the most recent results on high-energy hadropro duct ion of J /il), 
■0' and T's, let us remind the reader of the various kinds of collisions and colliders in which 
these particles can be created. 

For most of the following discussion, we shall limit ourselves to high-energy collider 
experiments, because most of the models that we shall consider in this thesis are valid in the 
perturbative regime of QCD, i.e. when the coupling constant ag is clearly less than 1; this is 
likely to be the case for reactions with cm. energy larger than, say, 5 GeV. Furthermore, the 
factorisation hypothesis relies on the assumption of having a sufficiently inclusive reaction, 
although the quarkonium is actually detected. This is another reason to impose a large 
energy. Concerning fixed target experiments, for which factorisation is uncertain, we shall 
not include them in our survey, even though heavy quarkonium production in that regime 
is certainly of great interest. 

Limiting ourselves to high-energy collisions, the most recent published results for quarko- 
nium hadroproduction come from two accelerators: 

1. The Tevatron at Fermilab which - as stated by its name - runs at TeV energy with 
proton-antiproton collisions. For Run I ("Run lA" 1993-94 and "Run IB" 1995-96), 
the energy in cm. was 1.8 TeV. For Run II, it has been increased to 1.96 TeV. The 
experimental analyses for this Run are still being carried out. 

2. RHIC at BNL running at 200 GeV for the J/il) study with proton-proton collisions. 

There exist also data from electron-proton collision from HERA at DESY running with 
electron-proton beams with proton energy of 820 GeV and electrons energy of 27.6 GeV. 

Results for e-p collisions come for the most part from two collaborations, namely ZEUS 
and HI. Two classes of processes can be distinguished depending on the off-shellness of the 
photon emitted by the scattering electron. If the photon is on-shell, we are dealing with 
photoproduction, otherwise, we deal with electroproduction. Results on photoproduction 
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are compatible with NLO calculations within CSM [53]. Occasionally, we shall mention 
these in the following discussion. 

Finally, three e^e~ collider are yielding data on quarkonium at present time: 

1. LEP at CERN running from around Zq mass up to around 210 GeV with electron- 
positron collisions. 

2. B factories, PEP II at SLAG and KEK B at KEK running mainly at the T(45) (10.58 
GeV) resonance with electron-positron collisions. 

However, quarkonium physics at LEP is very different from that at B factories. Due to 
its large cm. energy, LEP can significantly produce ip through 77 collisions which can be 
selected via kinematical constraints. From a theoretical point of view, this is rather interest- 
ing since this way of producing ip is similar to those at work in hadropro duct ion, although 
potentially better understood. On the other hand, to what concerns B factories results, we 
think that the energy is such that the reactions are not sufficiently inclusive to be straight- 
forwardly compared to reactions that take place at HERA, RHIG or the Tevatron. However, 
it is very important to keep in mind that there exist serious disagreements between mea- 
surements and theoretical predictions of both GSM and GOM. Sometimes they are thought 
to come from missing contributions at end-point which arise when re-summation techniques 
are applied; this is indeed the sign that some kinds of "exclusive" effects become important 
at these low energies. For further information, the interested reader may have a look at [34] . 

Goncerning hadroproduction, the above-mentioned deviations between theoretical pre- 
dictions and measurements of the GDF collaboration could be actually considered as a 
confirmation of the prior UAl analysis [28]. The deviations were confirmed by D0 later 
on [54] . For the Tevatron, only GDF analysis will be presented in the following. Goncerning 
RHIG, we shall expose the analysis for J/ip in 200 GeV pp collisions from the PHENIX 
collaboration [55]. 

It is worth pointing out here that high-energy p-p and p-p collisions give similar results 
for the same kinematics, due to the small contribution of valence quarks. 

2.2 Different type of heavy quarkonium production 

As we have already explained, the detection of quarkonia proceeds via the identification of 
their leptonic decay products. We give in Tab. (2.1) the relative decay widths into muons. 



Meson 


r(/i+/i~)/r(tota/) 


J/^ 


0.0588 ±0.0010 


^(2^) 


0.0073 ± 0.0008 


T(15) 


0.0248 ± 0.0006 


T{2S) 


0.0131 ±0.0021 


T(3S) 


0.0181 ±0.0017 



Table 2.1: Table of branching ratios in dimuons [57]. 
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It is interesting to note that the CDF analysis that lead to the problem of quarkonium 
hadroproduction was in fact not devoted to production but to 6-quark production through 
its (weak) decay into ip [35]. The value of this decay was already quite well known [56] 
at that time and a measurement of the proportion of tp^s produced through this decay was 
expected to give the total number of b quarks produced. 

Their study succeeded and they extracted the percentage of ip^s produced through b 
decays and of ip^s promptly produced. The latter was unexpectedly large. Indeed, the 
production through b decays was expected to behave like a fragmentation process and to 
dominate at moderate and large pt- This was not the case and the slight disagreement 
between total production predictions [27] and UAl measurements [28] was blown away by a 
striking incompatibility between prompt i(j production predictions and CDF data. A factor 
60 of discrepancy was for instance found for ip' (see Fig. 1.8 (b)). 

In this case, no higher excited states was expected to contribute to its production (see 
Fig. 1.5), contrarily to the case of the J/ip which can be the product of a Xc (-P-state) elec- 
tromagnetic decay. In order to get a complete picture of the problem, the CDF collaboration 
decided to proceed to a study of the direct production of this meson [43]. This involved the 
detection of the photon emitted during the radiative decay, 

+ (2.1) 

They then found a direct production 30 times above theoretical predictions. 

Briefly, the problem of direct J/ tp production separation comprises three steps: 

• muon detection; 

• elimination of J/ip produced by hadrons containing b quarks ; 

• elimination of Xc radiative-decay production; 

The following table summarises the different processes to be discussed and the quantities 
linked. 



1"* step 


2"'^ step 


i""^ step 


Type 


Associated 
quantity 


pp ^ CC + X 


cc J/tP 
cc^ Xc 

cc^ V(2S') 


Xc J/V' + 7 
?/'(25) J/ip + X 


Direct prod. 
Prompt prod, by 

decay of Xc 
Prompt prod, by 
decay of ?/'(2S') 


mats, 


pp be + X 

6c — > cc + £^ + Di 

etc . . . 


cc — > J/^p 

CC^ Xc 


Xc J/ip + l 







Table 2.2: Different processes involved in J/ip production accompanied by quantities used in the 
following discussion. 



Let us explain the different fractions that appear in the latter table: 
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• F{p, ■ifj{2S)Y^^' is the ^rom^t fraction of J /il) that do not come x-, neither from tjj(2S), 
i.e. the direct fraction. 

• -F(/^)^'(25) is the prompt fraction of J/t/j that come from il){2S). 

• F{p)'ll^ is the prompt fraction of J /ip that come from x- 

• Fb (or fb) is the non-prompt fraction or equally the fraction that come for b quarks. 

Concerning %1){2S), due to the absence of stable higher excited states likely to decay 
into it, we have the following summary for the different processes to be discussed and the 
quantities linked: 



1^* step 


2'"^ step 


3'^'' step 


Type 


Associated 
quantity 


pp ~* CC + X 


cc ^ ip(2S) 




Direct/Prompt prod. 




pp bc + X 
be ~* cc + £^ + Di 
etc . . . 


cc ip(2S) 






MorFb) 



Table 2.3: Different processes involved in tp' production accompanied by quantities used in the 
following discussion. 



Let us explain the different fractions that appear in the latter table: 

• F(p)^^'^^^ is the prompt fraction of tp{2S), i.e. the direct production. 

• Fb (or fb) is the non-prompt fraction or equally the fraction that come for b quarks. 

Considering the similarities between the charmonium and bottomonium systems (see 
Fig. 1.6 for the spectrum), a substantive part of the discussion will be also devoted to it. As 
can be seen in Tab. (2.1), the leptonic branching ratio are also relatively high. The detection 
and the analysis of the bottomonia are therefore carried out in the same fashion. For the 
extraction of the direct production, the 6-quark feed-down is obviously not relevant, only 
the decays from stable higher resonances of the family are to be considered. 

All the quantities useful for the bottomonium discussion are summarised in the following 
table: 



2.3 CDF analysis for ip production cross sections 
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1"* step 


2"^^ step 


3'^'^ step 


Type 


Associated 
quantity 


pp^ bb + X 


bb~^T{nS) 
bb Xb 

bb T{n'S) 


Xb T(n5) +7 
T{nS) T{n'S)+X 


Direct prod. 

Prod, by 
decay of Xb 
Prod, by 
decay of T{n'S) 


direct 

pT(nS) 
^T(ti'S) 



Table 2.4: Different processes involved in T{nS) (n = 1,2,3) production accompanied by quan- 
tities used in the following discussion. 

Let us explain the different fractions that appear in the latter table: 

• is the direct fraction of T(?7,5'). 

• F^jI'^^'* is the fraction of T{nS) that come from Xh- 

• fraction of T{nS) that come from a higher T[n'S). 

2.3 CDF analysis for ip production cross sections^*^ 

2.3.1 Succinct description of the parts of the detector used in this 
analysis 

Without entering too much into the details, let us recall the main elements of the detector 
used for this analysis. 

• A superconducting solenoid; 

• A central muon detector (CMU: central muon system (or central muon chambers)). It 
covers the region^ |?7| < 0.6; 

• The silicon vertex detector (SVX) which among other purposes enables to detect a 
second vertex resulting from the decay of a beauty hadron in a '0. This detector 
directly surrounds the beam; 

• Two calorimeters: an electromagnetic one and a hadronic one. They cover the whole 
azimuthal angle (cf. Fig. 2.1) and spread till pseudorapidity rj of 4.2 in absolute 
value; 

^ Some plots shown in the foUowing are preUminary. These are displayed only for an illustrative purpose. 
The discussion would have been the same if the corresponding final plots had been available. The plots 
relative to final observables are however all published ones. 

^ In the following, ip stands for both J/tp and ip' . 
?7 = — ln(tan(6'/2)) where 9 is the polar angle of the particle in the cm. frame. 
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• A central tracking chamber (CTC) lying in the induction of the magnet. It enables the 
determination of the transverse momenta of the particles. It covers a pseudorapidity 
range of |?7| < 1.1. 

2.3.2 ip total production 

Muons are selected by a three-level trigger: 

• The level 1 requires two track segments in the CMU, separated by at least 5° in 
azimuth (cf. Fig. 2.1); 

• The level 2 requires that at least one of the muon track segments is matched in to a 
track found in the CTC by the central fast tracker, a hardware track-finding processor. 
These two tracks are then considered to come from a single particle; 

• The level 3 requires a pair of oppositely charged muons after full track reconstruction. 
Besides this selection, muons are to satisfy other criteria: 

• The CTC track is extrapolated to the CMU chambers and is required to lie within 
3 0" of the CMU hits, where a is the uncertainty in the extrapolated position due to 
multiple scatterings; 

• A nonzero energy deposit is required in the calorimeter tower situated in front of the 
muon-chamber segment; 

• To remain in the region of good trigger efficiency, both muons are required to have 
Pt > 2.0 GeV and one muon pt > 2.8 GeV. The two muon tracks are fitted with the 
requirement that they originate from a common point; 

• To remain in a region of good acceptance to ip decays, the ip candidate is required to 
satisfy |?7| < 0.6 and px > 5 GeV. 

• The invariant mass of the pair must lie within the region 3057 MeV < M(pair) < 
3137 MeV in order to be considered as coming from a J/ip and 3646 MeV < M{pair) < 
3726 MeV to be considered as coming from a ip'. 




Figure 2.1: General drawing defining the angles. 
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M(\|/) 



Ecm(Pair) 



Figure 2.2: Invariant mass distribution of the muon pair in which the decays. The cm. energy 



On the other hand, 60% of the muon tracks also have hits in the SVX which measures 
only r and and the resolution of their track impact parameter is (13 + 40/pT)/uni, where 
Pt (expressed in GeV) is the track momentum transverse to the beam line. 

The sample of pp collisions amounts to 17.8 ± 0.6 pb~^ at ^/s = 1.8 TeV. For J/tp, the 
considered sample consists however of 15.4 ± 0.6 pb^^ of integrated luminosity^. 

The number of ip candidate is therefore determined by fitting the mass distribution of 
the muons in the cm. frame (cf. Fig. 2.2 and Fig. 2.3 (a) and (b)) after subtraction of the 
noise, which is evaluated as follows. 

Various corrections are evaluated by Monte-Carlo simulation. They incorporate devia- 
tions due to processes as ip ^ fi~^fi~'-f. Muon momenta are smeared in order to simulate the 
finite detector resolution. As the measurement of the muon energy is a function of pt, data 
are fitted in separate px bins. 

Concerning the trigger efficiency etrigger? it is accounted for by weighting the event with 
1 /^trigger 5 which dcpcuds on Pt- 

The mass distribution is then fitted to the signal shape fixed by simulation and to a linear 
background. The fit yields the mass of the particle and a background estimate; from this, 
it is easy to find back the number of particles which has decayed and from the branching 
ratio, the number of produced particles. 

For the present study, the fits were reasonably good for each pj^-bin, the per degree of 
freedom ranging from 0.5 to 1.5. The measured width of the mass peak was from 17 MeV 
to 35 MeV for pr from 5 to 20 GeV. 

Approximatively, 22100 J/ip candidates and 800 ip' candidates above a background of 
1000 events were observed. 

The acceptance was also determined by Monte-Carlo simulation, ip^s were generated with 
uniform distribution in pt, r] and (p. A simulated detector was used and the kinematic cuts 
already discussed were applied to the generated events. This leads to an acceptance of 9% 
at Pt = 5 GeV which rises to a plateau of 28 % for px > 14 GeV. It is important to note 
that the acceptance depends on the ip polarisation. By simulation, the uncertainty on this 

This difference is due to tlie subtraction of data taken during a period of reduced level 3 tracking- 
efficiency. These data have however been taken into account for ipi'^S) after a correction derived from the 
J/ijj sample. 



amounts to ^p mass. 
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Figure 2.3: Resonance peak of J/ip (a) and 0' (b). 

polarisation is shown to produce a 15% variation of the acceptance at = 5 GeV and 5% 
at 20 GeV for prompt production, whereas the non-prompt production variation is from 6% 
to 3% for the same transverse momenta. 

The efficiency of the CMU segment reconstruction is measured to be 97.2 ± 1.2% using 
dimuon events recorded with a single muon trigger. For CTC efficiency reconstruction, it 
is determined by a Monte-Carlo technique and is 96.4 ± 2.8% for both CTC tracks. The 
efficiency of the CMU-CTC matching requirements is estimated from the number of J/ip 
before and after the matching requirements, to be 90.5 ± 1.0%. 

Systematics on the efficiency of the three triggers are evaluated by varying the functional 
form of the trigger efficiency. The J/xIj efficiency is then 97.0 ± 0.2% and the '0' one is 
92.3 ± 0.2%. 

The if) differential cross section is experimentally defined as 



dpT 



/ Cdt ■ Apr ' 



(2.2) 



where N{iIj) is the number of ip candidates in the Apx bin with efficiency corrections applied, 
and e (7^ ^trigger) is the product of the acceptances (kinematic and detector- related) , of the 
reconstruction efficiencies and of the corrections due to selection cuts. 
The integrated cross sections were then measured to be 



(T{J/ip) ■ B{J/i) /i+/i-) = 17.4 ± Q.l{stat.)tll{syst.) nb 



(2.3) 



a{'ilj{2S)) ■ B{i!{2S) /iV^) = 0.57 ± QM{stat.)tl^^{syst.) nb. (2.4) 
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where cr{ip) = a{pp — > tpX^pT^ip) > 5 GeV, \ri{ilj)\ < 0.6). The evolution oi dcr x B as a 
function of pt is shown in Fig. 2.4. 
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Figure 2.4: (nb/GeV) as a function of pt for J/^ (a) and for ip' (b). 



2.3.3 Disentangling prompt charmonia 

As aheady seen, prompt ip^s are the one that do not come from the decay of B mesons. Their 
production pattern has the distinctive feature compared to non-prompt one that there exists 
a measurable distance between the B production vertex and its decay into charmonium. 

To proceed, it is necessary to use the SVX, whose resolution is 40/im, whereas B lifetime 
is ctb ~ 450/im. Muons are constrained to come from the same point which is called the 
secondary vertex, as opposed to the primary vertex , which is the collision point of protons. 
Then the projection of the distance between these two vertices on the ip momentum, L^y, 
can be evaluated. It is converted into a proper time equivalent quantity using the formula, 
CT = prr<l'{^ — , where Fcorr is a correction factor, estimated by Monte-Carlo simulations. 



m(!/j) 



which links the ip boost factor /3t7 to this of the B [58]. 

The prompt component of the signal is parametrised by cr = (a single vertex), the 
component coming from B decay is represented by an exponential whose lifetime is r?, and 
is convoluted with the resolution function. 

The cr distribution (cf. Fig. 2.5) is fitted in each px-hm with an unbinned log-likelihood 
function. The noise is allowed to vary within the normalisation uncertainty extracted from 
the sidebands. The fraction of tp coming from fe, fbipr), is displayed as a function of in 
Fig. 2.6 (a). 
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Figure 2.5: cr distribution for J/'ip (a) and for ip' (b). The solid region is the noise contribution 
(of null lifetime), the white region (just below the peak) give the null-lifetime signal: 
the prompt production, and finally the hashed region (mainly on the right side) is 
due to B component over background. 



Uncertainties are evaluated as follows: a different choice of the fitting procedure implies 
a mean relative variation of ±0.9%, this is included in the systematics on f^. In the same 
fashion, a variation of the lifetime of -B's of 1 a modifies by 0.7%. 

The ip production cross section from B decays is thus extracted by multiplying^ ft^iPr) 
by the inclusive ip cross section. Multiplying the latter by (1 — fl^^ipr)), we obviously get 
the prompt-production cross section (cf. Fig. 2.6 (b)). 



2.3.4 Disentangling the direct production of J/iIj 

The problem here is to subtract the J/ ip coming from Xc decay, assuming that this is the only 
source of prompt J/ip besides the direct production after subtraction of FipYJ^^y the prompt 
fraction of J/ip that come from %1){2S). The latter is obtained from the ip{2S) cross section 
from the previous section and from Monte-Carlo simulation of the decays ipi^S) J/ipX 
where X = tttt, t] and vr*^. The delicate point here is the detection of the photon emitted 
during the radiative decay of the Xc- 

The events used in the analysis are collected by a dimuon trigger with same requirements 
as prior except than the lower bound on is 4 GeV instead of 5 GeV. 



^ In order to reduce statistical fluctuations, ft is fitted by a parabola weighted by the observed shape of 
the cross section [58]. 



2.3 CDF analysis for ip production cross sections 



45 




Figure 2.6: (a) Fraction of ip from B decay as a function of px', (b) -^B from the prompt 
component of [35]. The lines are the theoretical expectations based on the CSM 
(LO + fragment.)[40, 39]. 

With this selection, a sample of 34367 J/ip is then obtained from which 32642 ± 185 is 
the number of real J/ip when the estimated background is removed. 
The requirements are as follow: 

• an energy deposition of at least 1 GeV in a cell of the central electromagnetic calorime- 
ter; 

• a signal in the fiducial volume of the strip chambers (CES); 

• the absence of charged particles pointing to the photon candidate cell (the no-track 
cut). 

The direction of the photon is determined from the location of the signal in the strip 
chambers and from the event interaction point. All combinations of the J/ijj with all photon 
candidates that have passed these tests are made and the invariant-mass difference defined 
as AM = M{fi'^fi^j) — M{fj,~^fj,~) can then be evaluated and plotted (see Fig. 2.7). 

As expected, a clear peak from Xc decays is visible near AM = 400 MeV. Yet, distinct 
signals for Xc,i and Xc,2 are not resolved as the two states are separated by 45.6 MeV and as 
the mass resolution of the detector is predicted to be respectively 50 and 55 MeV. 

With regard to the background, its shape is simulated by replacing all charged particles 
(different from muons) by tt*^, ?7 and Kg with probabilities proportional to 4 : 2 : 1. The 
detector response which results from the decay of these tt", r] and Kg is obtained using 
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Figure 2.7: Number of events as a function of the mass difference after the selection procedure 
for the different p^-bins. The points are the data; the solid region is Monte-Carlo 
simulation of the background. The solid line is a fit resulting from the background 
added to the contribution expected from a gaussian distribution [43]. 



a Monte-Carlo simulation. The AM distribution for the background is consequently ob- 
tained by applying on these simulated events the above-mentioned procedure as for real Xc 
reconstruction. A simple test, which consists in trying to reproduce with this Monte-Carlo 
simulation the sideband signal, where there should not be any Xc signal, was carried out 
successfully. Indeed imposing an equal number of events, the shape of the two distribution 
agree quite well (see Fig. 2.8). 

Eventually, the AM distribution obtained from the data is fitted with a gaussian and 
with the background fixed by the procedure explained above but with a free normalisation. 
The parameter of the gaussian then lead to the number of signal events: 1230 ± 72 Xc- 

The analysis of the direct J/ip signal is done within four p^-bins: 4 < p'^J'^ < 6 , 
6 < p^^^' < 8 , 8 < p'^^^ < 10 and p^^^' > 10 GeV. The proportion of J ftp produced by Xc 
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Figure 2.8: Comparison between the AM distribution for dimuons in the J/tp sidebands, and the 
corresponding one predicted by the Monte-Carlo calculation; the two distributions are 
normalised to equal area and the vertical scale is arbitrary [43]. 



decay is given by: 

pJI'P = - (2 51 

^Jli, ^trk ^env 

where is the probabihty to reconstruct the photon once the J/ ip is found (or the photon 
acceptance), ej^^ is the efficiency of the no-track-cut and ej^^ the efficiency to reconstruct 
the photon in the presence of additional energy deposited in the calorimeter. The latter 
quantities as well as the acceptance of the photon are evaluated with the help of a Monte- 
Carlo simulation in the spirit of the previous section. For instance, this acceptance is 0.146± 
0m2{stat.) for p^!^ > 4.0 GeV, whereas el^ amounts to 97.9^1% and ej^^ to 96.5l|^%. 

The total systematic uncertainty is estimated to be ±18.9% and is correlated in the four 
bins [43]. 

For p'^^^ > 4.0 GeV and \ri'^^^\ < 0.6, the fraction of J/ifj from Xc is then 

pJ/^ = 27.4% ± 1.6%{stat.) ± 5.2%{syst.) (2.6) 

The last step now is the disentanglement of the prompt Xc production, that is the deter- 
mination of F{pyll''^ . Let us here draw the reader's attention to the fact that by selecting 
prompt J/ip (cf. section 2.3.3), J/ip produced by non-prompt Xc have also been eliminated; 
it is thus necessary to remove only the prompt production by Xc decay, and nothing else. 
This necessitates the knowledge of F(p)'^^^ , the prompt fraction of J/ip that come from x- 

The latter is calculated as follows: 



pay"* = ~ = ^-^t „ 71 
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where N^, N^'"'' are the number of reconstructed Xc and J/ip from 6's, F^, 
corresponding fractions. 

F^^"^ (or fb) has been extracted in the section 2.3.3; F^ is obtained in same way and is 
17.8% ± 0.45% for pr > 4.0 GeV. Consequently, 



F{pyj^ = 29.7% ± 1.7%{stat.) ± 5.7%{syst.) 



(2.8) 



and its evolution as a function of pt is shown in Fig. 2.9. It is also found that the fraction 
of directly produced J/ tp is 



= 64% ± 6%, 



(2.9) 



and is almost constant from 5 to 18 GeV in p^- We therefore conclude that direct production 
is the principal contribution to J/ip (see Fig. 2.10 ). 
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Figure 2.9: The fraction of J/ip from Xc as a function of px with the contribution of 6's removed. 

The error bars correspond to statistical uncertainty.The solid line is the parametrisa- 
tion of the fraction. The dashed lines show the upper and lower bounds corresponding 
to the statistical and systematic uncertainties combined [43]. 



In order to get the cross section of direct J/ ip production, it sufficient now to extract the 
contribution of obtained by Monte-Carlo simulation and of Xc obtained by multiplying 



the cross section of prompt production by the factor F(J))'^^'^ , which is a function of p^ 
The different cross sections are displayed in Fig. 2.11. 
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Figure 2.11: Differential cross section for prompt production of J/tjj — s> /i"/!"*" as a function of 
Pt- The dashed line gives the prediction of the CSM (LO+Fragmentation) [40, 39] 
and the solid lines that of the CSM + preliminary COM calculations for J/tp from 
Xc (see [43]), to be compared respectively with • and ■. 
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2.4 CDF measurement of the T production cross sections 

As T states are similar to ip states for which data are numerous and in contradiction with 
predictions, their study is just as good a point of comparison. Furthermore, many approxi- 
mations performed in the models are better for hb states. Firstly, there is no expected decay 
from t quarks; secondly, the scale set by the quark mass is higher and therefore perturbation 
methods should be more suitable; thirdly, fragmentation would be significant only for large 
Pt- Classical production mechanisms based on LO pQCD can thus be considered as the 
only source of production; finally, the velocity of the heavy quarks within the meson being 
smaller, relativistic effects should be smaller, rendering the predictions of non-relativistic 
models more reliable. 

In this section, we shall sum up the procedure and the results on T production in pp 
at ^/s = 1.8 TeV. These results were exposed in two letters [59, 60], and we shall mainly 
focus on the second one, which considered data collected in 1993-95 and corresponding to 
an integrated luminosity of 77 ± 3 pb~^. 

The parts of the detector used are the same as for the -j/^'s, with an additional muon 
detector, called central muon upgrade (CMP) outside the CMU. It covers the same rapidity 
range, i.e. \ri\ < 0.6, but it is located behind an additional steel shielding of ~ 2.4 interaction 
lengths. This increases its signal purity, although some muons from the CMU could have 
too few energy to reach it. 

2.4.1 T total production 

To what concerns muons, they are selected by a three level trigger: 

• The level 1 requires two track segments in the CMU ( with efficiencies ei{pT = 
1.5 GeV) = 40% - ei{pT > 3.0 GeV) = 93%); 

• The level 2 requires that at least one of the muon-track segments be matched in 
(within 5°) to a track found in the CTC by the central fast tracker. These two tracks are 
then considered to come from a single particle (the efficiency is here (ei(pT > 2.5 GeV)) = 
95%); 

• The level 3 performs a 3D reconstruction and requires that muons pairs have an 
invariant mass in the region from 8.5 to 11.5 GeV. 

Besides the selection, muons are to satisfy other criteria to isolate T resonances: 

• Identification by CMU of both muons is required as well as at least one identification 
by the CMP; 

• The momentum of each muons is determined using the CTC information combined 
with the constraint that they come from the beam line (no secondary vertex possible 
here) ; 

• Only muons with pt > 2.2 GeV were selected; 
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• Each muon-chamber track is required to match the extrapolation of a CTC to within 
3 0" in r — and 3.5 a in z {a is the calculated uncertainty due to multiple scatterings 
and detection uncertainties, exactly as for the ip^s.); 

• Muons are required to have opposite charge and the rapidity of the reconstructed pair 
is to be in the CMU fiducial rapidity region \y\ < 0.4. 

• The transverse momentum of the pair is finally restricted to be between and 20 GeV. 
The mass distribution obtained by CDF is shown in Fig. 2.12. 
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Figure 2.12: Resonance peaks of T(IS'), T(2S'), T(3S'). The histogram represents the data and 
the solid curve a gaussian fit to each resonance plus a quadratic background [60]. 



The differential cross section times the branching ratio in dimuon is experimentally ob- 
tained as follow: 

/^\ g(T^^V)= . (2.10) 

where Nfn is the number of T extracted from the gaussian fit on the signal for each p^-bins 
similar to Fig. 2.12, A is the geometric and kinematic acceptance determined by Monte-Carlo 
simulation, e represents the various efficiency corrections detailed in [60]. A is comparable 
for the three states and varies from 17% to 24% as a function of pr- 

Systematic uncertainties on ^^^^ arise from C (4.1%), the level 1 and 2 trigger-efficiency 
correction etrig (4%) and from the remaining efficiency corrections (6%). Uncertainties are 
also due to the fact that the T(2S') and T(3S') polarisations are unknown. They were 
evaluated by running the simulation with extreme value of polarisation for these states 
(longitudinal or transverse). For the T(IS'), as the polarisation is measured (cf. section 2.5), 
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the limits obtained are used as extrema for this variation and it is seen that the uncertainties 
are then neghgible and no systematic error from this source is then to be added. 
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Figure 2.13: Differential cross section of T(15') ^ fi /i"*" as a function of fof \y\ < 0.4. 
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Figure 2.14: Differential cross section of T(15') ^ fi /i"^ as a function of for \y\ < 0.4. 
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Figure 2.15: Differential cross section of T(15') ^ ^ yU+ as a function of for \y\ < 0.4. 



2.4.2 Disentangling direct production of T(IS') 

The analysis [61] presented here is for the most part the same as described in section 2.3.4. 
It is based on 90 pb~^ of data collected during the 1994-1995 run. The measurement has 
been constrained to the range px > 8.0 GeV because the energy of the photon emitted 
during the decay of Xb decreases with px and ends up to be too small for the photon to be 
detected properly. In the same spirit, analysis relative to T{2S) has not been carried out, 
once again because of the lower energy of the radiative decay. Concerning T{3S), except for 
the unobserved X6(3-P) (which is nevertheless supposed to be below the BB threshold), no 
P-states are supposed to be parents. 

The event used are collected by a dimuon trigger (details are given in [61]) with the 
following requirements, which slightly differ from section 2.4.1: 

• Two oppositely charged muons candidates observed in the muon chambers, 

• Both with Pt > 2.0 GeV, 

• The pair must have px > 8.0 GeV and |?7| < 0.7, 

• In order to be considered as coming from a T(IS'), the invariant mass of the pair must 
lie within the region: 



9300 MeV < M{pair) < 9600 MeV. 



(2.11) 
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With this selection, a sample of 2186 T(15') candidates is then obtained from which 
1462 ± 55 is the estimated number of T(IS') after subtraction of the background. The 
resulting mass distribution is plotted in Fig. 2.16. 

In this sample, photons likely to come from Xb decay are selected as for the J/ip case, 
except for the energy deposition in the central electromagnetic calorimeter, which is lowered 
to 0.7 GeV. 
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Figure 2.16: Muon pairs invariant mass distribution after the selection for pt > 8.0 GeV. The 
dotted line is the fit to the background, the solid line is the fit to the data; the 
region S is the signal region and the region B is the sideband region. 



The direction of the photon is determined from the location of the signal in the strip 
chambers and from the event-interaction point. All combinations of the T(15') with all 
photon candidates that have passed these tests are made and the invariant-mass difference 
still defined by AM = M(/i+/x~7)— M(/i+/i~) can then be evaluated. It is plotted in Fig. 2.17. 

The background is evaluated as in the J/ip analysis. A fit to the data then yields the 
number of signal events: 35.3 ± 9.0 for Xbi^P) and 28.5 ± 12.0 for Xfe(2P). The analysis 
of the direct T(15') signal is done within one single p^-bin, as opposed to J/tp case. The 
proportion of T(IS') produced by Xb decay is given by: 

f?J/^ - i! (2 12) 

where N^'' and A^^ are the numbers of reconstructed Xb and T(IS') respectively, is 
the probability to reconstruct the photon once the T is found, e'' is the efficiency of the 
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Figure 2.17: Number of events as a function of the mass difference after tine selection proce- 
dure. The points are the data; the solid region is a Monte-Carlo simulation of the 
background. The solid line is a fit resulting from the background added to the 
contribution expected from two gaussian distributions [61]. The inset shows the 
comparison between the background of A A/ extracted from the sidebands and with 
the Monte-Carlo simulation. 



isolation cuts. The latter quantities, as well as the acceptance of the photon, are evaluated 
with the help of simulations in the same spirit as that described in the previous section. 
For instance, this acceptance is 0.142 ± O.OOA{stat.) for Xbi^P) and 0.284 ± 0M6{stat.) for 
X6(2P), the difference between the two coming mainly from the lower energy of the photon 
for Xb(2P) decays. The efficiencies e'^ obtained are 0.627 ± 0.013(stat.) for Xbi^P) and 
0.651 ± 0.0013(stat.) for Xfe(2P). The efficiency for inclusive T(15') is supposed to be the 
same as for the sub-sample considered. 

The systematic uncertainties on F^^ associated with the variation of the shape of 
the pt spectrum and of the decay angular distribution is ±13% for Xb(lP) and ±9% for 
Xb(2P). The uncertainty from N^'' is ±7% for Xb(lP) and ±9% for Xfe(2P), this comes from 
the background of the AM distribution and, from the uncertainty on the resolution of the 
gaussian signals in this same distribution. Other uncertainties are also considered, see [61] 
for more details. Assuming the independence of all these uncertainties, they are combined 
in a total systematic uncertainty: ±16.4% for ^^(IP) and ±13.7% for Xb(2P). 

For p^^^^^ > 8.0 GeV and |?7^^^'^^| < 0.7, the fractions of T(IS') from Xb(lP) and from 
Xb(2P) are 

Pji!p\ = 27.1%±6.9%{stat.)±4A%{syst.), 
F^%% = 10.5% ± AA%{stat.) ± lA%{syst.). 
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In order to get the direct production of T(15'), one must still evaluate the feed-down 
from the S*- waves T{2S) and T{3S). This is done using Monte-Carlo simulations of these 
decays normalised to the production cross sections discussed in section 2.4. It is found that 
for p^^^^^ > 8.0 GeV the fraction of T(IS') from T{2S) and T{3S) are respectively 

pT{15) _ ^ ocr/+0.6q/ ^ ' ^ 

-^T(35) - U.«/o_o,4/o. 

Concerning the unobserved Xb(3-P), a maximal additional contribution is taken into ac- 
count by supposing than all the T{3S) are due to Xfe(3-P) and from theoretical expectation 
for the decay of this state, a relative rate of T(IS') from Xb(3P) can obtained. This rate is 
less than 6%. 

In conclusion, the fraction of directly produced T(IS') is obtained from 

pT(lS) _ _ pT(lS) _ pT{15) _ pT(lS) _ pT{15) . p.^ 

^direct Xb{lP) XbCiP) ^r{2S) -^T(35) ' y^-lO) 

for which the uncertainties come from that of cross sections and of branching ratio. It is 
eventually found that 

fJ^^I^J = 50.9% ± 8.2%{stat.) ± 9.0%{syst.). (2.16) 



2.5 Polarisation at the Tevatron 

As the considered bound states, '0 and T are massive spin-1 particles, they have three 
polarisations. In addition to measurements of their momentum after their production, the 
experimental set-up is sufficiently refined to provide us with a measurement of their spin 
alignment through an analysis of the angular distribution of the dimuon pairs from the 
decay. 

The CDF collaboration has carried out two analyses, one for the if) states [52] - for which 
the feed-down from b decay has been subtracted, but not the indirect component in the case 
of the J/i/j - and another for T(?7.5') [60]. 

In the following, we shall proceed to a brief outline of the analysis and give its main 
results. 



2.5.1 Polarisation parameters 

The polarisation state of the quarkonium can be deduced from the angular dependence of 
its decay into ^~ ■ Taking the spin quantisation axis along the quarkonium momentum 
direction in the pp cm. frame, we define 9 as the angle between the ix^ direction in the 
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quarkonium rest frame and the quarkonium direction in the lab frame (see Fig. 2.18). Then 
the normahsed angular distribution /(cos 6) is given^ by 

3 



I {cos 6) 



2 



2(a + 3) 



a COS' 



(2.17) 



where the interesting quantity is 



a 



(2.18) 

-1 corresponds to a full transverse 



a = means that the mesons are unpolarised, a - 
polarisation and a = —1 to a longitudinal one. 

As the expected behaviour is biased by muons cuts - for instance there exists a severe 
reduction of the acceptance as 9 approaches and 180 degrees, due to the pt cuts on the 
muons the method followed by CDF was to compare measurements, not with a possible 
(1 + a cos^ 9) distribution, but with distributions obtained after simulations of quarkonium 
decays taking account the geometric and kinematic acceptance of the detector as well as the 
reconstruction efficiency. 

virtual Q direction 



Q rest frame 




lab direction 



Figure 2.18: Definition of the angle 9 used in the polarisation analysis of a quarkonium Q. 

The parts of the detector used are the same as before, with the additional central muon 
upgrade (CMP) outside the CMU. 



2.5.2 Study of the il)'s 

We give here the results relative to the CDF analysis published in 2000 [52]. The data used 
correspond to an integrated luminosity of 110 pb~^ collected between 1992 and 1995. The 
three-level trigger is similar as before. Beside this on-line selection, muons are to satisfy 
other criteria to isolate if) resonances: 

^ For a derivation, see the Appendix A of [62]. 
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• Identification by CMU of botli muons is required as well as a matcliing of the corre- 
sponding tracks witliin tfie CMU witli one of tlie CTC tracks. This reduces the rate 
from sources such a.s n/K meson decay- in-flight; 

• These matches are required to pass a maximum cut of 9 and 12 (for 1 degree of 
freedom) in the and z views respectively [52]; 

• Only muons with pt > 2 GeV were selected; 
Disentangling prompt production 

The method followed here is on many aspects similar to the one exposed in section 2.3.3. 
As before, prompt events have cr consistent with zero whereas 6-quark decays have an 
exponential cr distribution, all this with a smearing due to the detector resolution. The 
CT distribution is fitted to obtain the relative fractions of prompt and fe-decay production. 
Results are that the measured fraction of J/ip mesons which come from 6-hadron decay, 
f//'^, increases from (13.0 ±0.3)% at pt = 4 GeV to (40 ±2)% at 20 GeV and for V' mesons, 
Ff is (21 ± 2)% at 5.5 GeV and (35 ± 4)% at 20 GeV. 

The following study is then based on two samples a short-lived one dominated by prompt 
production, and a long-lived one dominated by 6-decays. The definitions used to separate 
the samples are the following —0.1 < cr < 0.013 [0.01] cm for the short-lived sample, and for 
the long-lived sample 0.013[0.01] < cr < 0.3 cm, for the J/ip [ip'] analyses respectively'^. 

Within a 3-standard-deviation mass window around the J/ip peak, the data sample is 
of 180 000 J/ip events. In order to study the effect of pr, the data are divided into seven 
PT-bins from 4 to 20 GeV. Because the number of ip' events is lower, data for tp' are divided 
into three p^-bins from 5.5 to 20 GeV. 

The J/ip sample is divided into three sub-samples where the analysis is executed inde- 
pendently: the short-lived and long-lived SVX samples described above, and a third sample 
(the CTC sample) in which neither muon has SVX information and no cr measurement is 
made. 

For tp', both muons are required to be reconstructed in the SVX and the number of signal 
events is found to 1855 ±65 in a 3-standard-deviation mass window around the ip' mass. The 
sample is further divided into two sub-samples based on the cr distribution instead of 3 for 
J/ip. Once again, because the number of events is lower than in the J/ip case, 10 bins are 
used in | cos6'| instead of 20 in cos6 for J/ip. Similarly to the latter, the number of events 
in each bin is obtained by fitting the mass distribution each time. 

J/ 4' polarisation measurement 

The J/ip polarisation is determined through a fit of the measured muonic cos^ distributions 
to a set of Monte-Carlo templates. Implementing the experimental constraints with detector 

^ The boundary between the two cr regions has been determined separately for the J/ip and for the ip' . 
Depending on p^, the prompt fraction in the short-Hved sample ranges from 85%[86%] to 96% [95%], and the 
-B-decay fraction in the long-lived sample ranges from 83%[86%] to 98%[91%], for J/i{j [ip'] respectively. 
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and trigger simulations, these templates are generated by processing simulated samples of 
J/ip decays. 
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Figure 2.19: The points represent the data (with sideband signal subtracted); the dashed lines is 
the fit of the cos 6* distributions in the 12-15 GeV bin. The three plots correspond 
to the three samples defined in the text. 

The polarisation is obtained using a fit of the data to a weighted sum of transversely 
polarised and longitudinally polarised templates. The weight obtained with the fit provides 
us with the polarisation. Explanations relative to procedure used can be found in [52] and 
[62]. 

Note however that the polarisation is measured in each pt bin and that separate polari- 
sation measurements for direct J/ ip production and for production via Xc and if)' decays was 
found to be unfeasible. Let us recall here that Xc and if)' were shown to account for 36±6% of 
the prompt production (cf. Eq. (2.9)) and to be mostly constant in the considered px range. 
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For J/ip, the fit is performed simultaneously on tlie SVX sliort-lived, SVX long-lived 
and CTC samples, with two fit parameters: aprompt and afromB- As observed above, the 
long-lived and short-lived samples are not pure, this is of course taken into account. The 
-B-decay fraction in the CTC sample is assumed to be the same as in the SVX sample. This 
seems reasonable as the main difference between the two samples resides primarily in the z 
position of the primary vertex. The fit in the pt range 12-15 GeV is shown in Fig. 2.19. 

Except in the lowest pr bins, the systematic uncertainties are much smaller than the 
statistical one. The values obtained for aprompt and afromB are given in Tab. (2.5) and 
plotted with a comparison to NRQCD predictions in Fig. 2.20. 



Pt bin (GeV) 


Mean Pt (GeV) 


^Prompt 


afromB 


4- 5 

5- 6 

6- 8 
8-10 
10 - 12 
12 - 15 
15 - 20 


4.5 

5.5 

6.9 

8.8 

10.8 

13.2 

16.7 


0.30 ±0.12 ±0.12 
0.01 ±0.10 ±0.07 
0.178 ±0.072 ±0.036 
0.323 ±0.094 ±0.019 
0.26 ±0.14 ±0.02 
0.11 ±0.17±0.01 
-0.29 ±0.23 ±0.03 


-0.49 ±0.41 ±0.13 
-0.18 ±0.33 ±0.07 
0.10 ±0.20 ±0.04 
-0.06 ±0.20 ±0.02 
-0.19 ±0.23 ±0.02 
0.11±[]:i ±0.02 

-0.16 ±[5:1 ±0.05 



Table 2.5: Fit results for J/ip polarisation, with statistical and systematic uncertainties. 
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Figure 2.20: aprompt ^or J/^ fitted for \y^^^\ < 0.6. Full error bars denote statistical and 
systematic uncertainties added in quadrature; ticks denote statistical errors alone. 
The shaded band shows a NRQCD-factorisation prediction [64] which includes the 
contribution from Xc and ^' decays. 
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V'' polarisation measurement 

The procedure to obtain aprompt and ajromB is similar. Weighted simulations of the angular 
distribution are fitted to the data and the weight obtained gives the polarisation. The | cos ^| 
distributions in the two cr sub-samples are fitted simultaneously. As there is no expected 
radiative decay from higher excited charmonia, we are in fact dealing with direct production. 
The fit of the angular distribution in the three p-r-bins is shown on Fig. 2.21. 
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Figure 2.21: The points represent the ip' data for the short-lived sample only; the dashed lines is 
the fit to the | cos^| distributions. The three plots correspond to the three samples 
in pt defined in the text. Note the extension of the acceptance range in |cos^| as 
Pt increases. 



Anew, the systematic uncertainties [52] are much smaller than the statistical uncertain- 
ties. The values obtained for aprompt and afromB are given in Tab. (2.6) and plotted with a 
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comparison to NRQCD predictions in Fig. 2.22. 
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Figure 2.22: aprom.pt fof ^' fitted for \y'^'\ < 0.6. Error bars denote statistical and systematic 
uncertainties added in quadrature. Shaded bands show two NRQCD-factorisation 
predictions [63, 64]. 



Pt bin (GeV) 


Mean pt (GeV) 


(^Prompt 


^framB 


5.5 - 7.0 


6.2 


-0.08 ±0.63 ±0.02 


-0.26 ± 1.26 ±0.04 


7.0 - 9.0 


7.9 


0.50 ±0.76 ±0.04 


-1.68 ±0.55 ±0.12 


9.0-20.0 


11.6 


-0.54 ±0.48 ±0.04 


0.27 ±0.81 ±0.06 



Table 2.6: Fit results for ip{2S) polarisation, with statistical and systematic uncertainties. 



2.5.3 Study of the T(15') 

As the procedure is quite similar as of the V' study, we just present here the results. The 
measurements are made in the region of pr from to 20 GeV and with \y\ < 0.4 and the 
data are separated in four p^-bins [60]. 

Using the same technique as before, the longitudinal polarised fraction is obtained. 
Fig. 2.23 shows the cos6' distribution with the two templates used, and the resulting fit. 
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Figure 2.23: Uncorrected cos6' distribution for < 0.4 and for the highest p^-bin. The 

points are the data; the solid line, the result giving the weight of the transversal 
template (dashed line) and the longitudinal template (dotted line) [60]. 



In Tab. (2.7) are given the results for a and the same values are plotted in Fig. 2.24 in 
comparison with a theoretical prediction from NRQCD [65]. In conclusion, T(15') are shown 
to be mostly unpolarised. 




Figure 2.24: a for T(15) fitted for \y^^^^\ < 0.4. The theoretical band represents the NRQCD 
prediction [65]. 
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Pt bin (GeV) 


a 


0.0 - 3.0 
3.0 - 5.0 
5.0-8.0 
8.0-20.0 


+0.05 ±0.14 
+0.01 ±0.14 
+0.10 + 0.17 
-0.12 + 0.22 



Table 2.7: Fit results for T(IS') polarisation. 



2.6 PHENIX analysis for J/ijj production cross sections 

In this section, we present the first measurements of — > J /ip + X at RHIC obtained by 
the PHENIX experiment [55] at a cm. energy of ^/s = 200 GeV. The analysis was carried 
out by detecting either dielectron or dimuon pairs. The data were taken during the run at 
the end of 2001 and at the beginning of 2002. The data amounted to 67 nb~^ for ^'^ ^~ 
and^ 82 nb~^ for e+e~. The 5-decay feed-down is estimated to contribute less than 4% at 

= 200 GeV and is not studied separately. 

Concerning the dielectron study, it relied on: 

• Two central spectrometer arms covering |?7| < 0.35 and A0 = 90°; 

• An electromagnetic calorimeter (EMC); 

• A ring imaging Cerenkov detector (RICH) with a rejection threshold for pions of 
4.7 GeV. 

whereas for muons, were used: 

• Two muon-spectrometer arms covering A<^ = 360°, of which only one was used (1.2 < 
V < 2.2); 

• A muon-identifier system (MuID); 

• A muon-tracker system (MuTr). 

Concerning dielectrons, the trigger requirement was either a deposition of 0.75 GeV in a 
2x2 pile of EMC towers or 2.1 GeV in a 4 x 4 pile. Electron candidates had to be associated 
with a RICH ring (with at least 2 photo-tubes hit) and with an EMC hit (within +4(T 
for the position). Furthermore, a correlation between the energy deposition in the EMC, 
E, and the momentum of the reconstructed track, p, was required through the condition 
0.5 < E/p < 1.5. Finally, to avoid pion signals from the RICH, a 5 GeV cut-off for the 
electrons was imposed. 

On the other hand, the muon trigger required two deeply penetrating track seeds ( "roads" ) 
in two distinct azimuthal quadrants of the MuID [66]. Then these seeds were matched to the 

® The difference results from different cuts on the extrapolated vertex position. 
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cluster hits in the MuTr stations. A fit with energy-loss correction of the MuID and MuTr 
hit positions and of the collision vertex finally gave the momentum. 

The net yield of J/ip within the region 2.80 GeV < M{pair) < 3.40 GeV was found 
to be 46.0 ± 7.4 for electrons, whereas for muons, it was 65.0 ± 9.5 J/tp within the region 
2.71 GeV < M{pair) < 3.67 GeV. The mass distributions obtained are shown in Fig. 2.25 
for both dimuon and dielectron selections. As an indication of the background, a like-sign 
pair distribution is shown by dashed lines. 
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Figure 2.25: The invariant mass distribution for dielectron (left) and dimuon (right) pairs. The 
solid lines are for unlike-sign pairs and the dashed lines for like-sign pairs (from [55]). 



The cross sections were in turn extracted using 



B 



(fa 



N 



dydpT (J Cdt)AyApTebias(^trig Aer 



(2.19) 



where Nj/^ is the measured J/ip yield, / Cdt the integrated luminosity, etias is the minimum 
bias trigger which is used to measure the luminosity, etrig is the trigger efficiency, A the 
acceptance for J/ip reconstruction and finally erec its efficiency. 

We shall not describe here the method to evaluate the systematics, we simply reproduce 
verbatim the table given in [55] . Note however that the analysis of errors did not take into 
account possible polarisation effects. 
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Quantity 






Yield 


± 5% 


± 5% 




0.026-0.010 ± 13% 


0.038 - 0.017 ± 13% 




(2x2) 0.87-0.90 ± 5% 






(4x4) 0.30-0.74 ± 36% 




^bias 


0.75 ± 3% 


0.75 ± 3% 


J Cdt 


82 nb-i ± 9.6% 


67 nb-i ± 9.6% 


Total 


± 15%(sys) ± 10%(abs) 


± 14%(sys) ±10%(abs) 



Table 2.8: Table of quantities and their systematic-error estimates. Ranges are given for px 
dependent quantities. For the yU^/i" case, the values of Aerec and ^trig ai'e combined. 
The absolute cross-section normalisation uncertainty from e^ias and J Cdt is kept 
separate and is labeled "(abs)" (taken from [55]). 



From Eq. (2.19), the cross section as a function of px is obtained for the two analyses. 
For pt larger than two 2 GeV, it is compared in Fig. 2.26 with LO GSM and COM without 
fragmentation but with Xc feed-down. These theoretical predictions are from Nayak et 
al. [67]. 
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Figure 2.26: The differential cross section measured to be compared with theoretical predic- 
tions from CSM and COM (without fragmentation but with Xc feed-down) of [67] 
(adapted from [55]). 



Chapter 3 

Building a new model 



3.1 Foreword 

We first present a study [68] of tlie tlieoretical uncertainties of the LO CSM as an indication 
of tfie uncertainties of all theoretical model. As we shall see, these are far from being sufficient 
to explain the discrepancies with experimental measurements. 

These theoretical uncertainties can arise from errors on the leptonic decay width deter- 
mination. The latter is used to determine the wave function at the origin, which is the 
non-perturbative input of the model. Similar inputs exist for the COM and they are asso- 
ciated with the probability for the QQ pair to eventually evolve into a given physical state. 
These probabilities are obtained from fits and are subject to similar uncertainties. 

Another main source of uncertainty is the parton distribution function. In the case at 
hand, only gluons are actually involved, as the energy of the colliding hadrons is very large, 
and the required momentum fraction of partons consequently low. 

After the presentation of this study, we shall expose the guidelines that we shall follow 
in order to build a new model; in a sense, the model serves as a non-static extension of the 
CSM, but it could also be seen as a generalisation of the COM for QQg states with the gluon 
having a large momentum, whereas, in NRQCD, the momentum of the gluon associated with 
the heavy quark pair in the considered Fock state is necessarily lower than Mv where v is 
the typical velocity of the quarks within the bound state of mass M. 



3.2 Theoretical uncertainties of the LO CSM 
3.2.1 Wave function 

At LO in a^, the leptonic decay width reads: 

r(^^i - = ^ '-j^. (3.1) 

where a is the fine structure constant and the quark charge in units of e. This is one of 
the consequences of the factorisation postulate. 
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We thus find that the error on r35j^^^ introduces an error of at least 10% on the cross 
section. 



Meson 


|^(0)|2±a|^(0)|2 


Relative error 


J/i, 


0.041 ± 0.0042 GeV^ 


10% 




0.024 ± 0.0025 GeV^ 


10% 


TilS) 


0.397 ±0.015 GeV^ 


4% 


T(2S) 


0.192 ±0.030 GeV^ 


16% 


T(35) 


0.173 ± 0.032 GeV^ 


18% 



Besides, one may consider the Colour Singlet Model as a level of approximation of the Non 
Relativistic Quantum Chromodynamics (NRQCD) (already exposed in the first chapter). 
Indeed this formalism describes the production of quarkonia by two main classes of processes, 
the first when the quark pair is perturbatively produced within a colour- singlet configuration 
- it can then have exactly the same quantum number as the meson, as in the GSM, or be 
in a different spin configuration-, the second class is when the quark pair is perturbatively 
produced within a colour-octet configuration. 

In fact, as already mentioned, NRQGD predicts that there is an infinite number of Fock- 
state contributions to the production of a quarkonium Q. Practically, one is driven to 
truncate the series; this is quite natural in fact since most of the contributions are known to 

2 

be suppressed by factor of C(^), where v is the quark velocity in the bound state. 

For definiteness, the latest studies accommodating the production rate of ^5*1 states retain 
only the colour singlet state with the same quantum numbers as the bound state and colour 
octet P-wave states and singlet S'-wave. In this context, the GSM can be thought as a 
further approximation to the NRQGD formalism. 

However, in this formalism, factorisation tells us that each contribution is the product 
of a perturbative part and a non-perturbative matrix element, giving, roughly speaking, 
the probability that the quark pair perturbatively produced will evolve into the considered 
physical bound state. If one transposes this to the GSM, this means the wave function at 
the origin corresponds to this non-perturbative element. This seems reasonable since the 
wave function squared is also a probability. 

Yet, one has to be cautious because one actually links production processes with decay 
processes. In NRQGD, two different matrix elements are defined for the "colour singlet" 
production and decay, and they are likely to be different and independent. 

The only path left to recover the GSM is the use of a further approximation, the vacuum 
saturation approximation. The latter tells us that the matrix element for the decay is linked 
to the one for the production. This enables us to relate the wave function at the origin 
appearing in Eq. (3.1) to the colour-singlet NRQGD matrix element for production. This 
gives: 

(0t{^S,)) = lS\m\^ + 0{v'). (3.2) 

The conclusion that could be drawn within NRQGD is that the extraction of non- 
perturbative input for production from the one for decay is polluted by factors of O(w^). 



3.2 Theoretical uncertainties of the LO CSM 
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3.2.2 Parton distribution functions 

Let us now analyse the uncertainties due to the parton distribution functions (pdf). As 
there are two gluons coming from the initial hadronic state, these can induce non-negligible 
uncertainties on the cross section, either directly or via the associated tts(Q^). 

For coherence, we have limited ourselves to LO pdf's. As a consequence, we have not 
felt the necessity to refine the study to the point of using the very last pdf on the market. 
This choice is further justified because we want to have an estimate of the effects due to 
extreme-case choices. Here are different plots obtained for a set of 12 pdf's (see the legend 
of Fig. 3.1) coded in PDFLIB [69] which illustrate these effects: 
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Figure 3.1: Variation of the cross sections obtained with various pdf's [69] due to the pdf's only. 
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Figure 3.2: Variation of the cross sections due to as only. 
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Figure 3.3: Variation of the cross sections obtained with various pdf's [69] due to «s and the 
pdf 's. 




Figure 3.4: Ratio of the cross sections obtained with various pdf's [69] to the cross section 
obtained with MRS(G) 2-94. 



3.3 Assumptions for a new model 
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Figure 3.5: Variation of the cross sections obtained with various pdf's [69] for T{3S). 

We thus found that the overall factor resulting from these 3 sources of uncertainties is 
at most 2-3 and roughly constant for different values of pt- Hence, the discrepancy of the 
CSM with data is not due to conventional theoretical errors. Furthermore, these could not 
change polarisation results either, for which data seem to contradict NRQCD. 

3.3 Assumptions for a new model 

As we have already stated, no solution of the polarisation issue has been proposed so far, 
while we just have seen that conventional theoretical uncertainties are under control. We 
therefore propose to go beyond the traditional treatment which is to consider the heavy 
quarkonia as purely non-relativistic states with the consequence of neglecting a priori any 
effects that may come from the internal dynamics of quarks in the bound state. This approxi- 
mation is sensible and supported by arguments such as the success of potential models, with 
instantaneous interactions. Nevertheless, it has never been demonstrated that non-static 
effects are systematically negligible and that the internal dynamics as well as the quark 
off-shellness in the meson do not provide any sizable modification of their description. 

In order to tackle internal dynamics, we shall describe the heavy quarkonia with a phe- 
nomenological 3-particle vertex, made of a Dirac structure inferred from the quantum num- 
bers of the considered quarkonium and a vertex function inspired by the Bethe-Salpeter 
formalism. The description and the justification of the different steps leading to this 3- 
particle vertex is made in the next chapter. 

Nevertheless, as we shall see in chapter 5, the introduction of vertex functions breaks 
gauge invariance in production processes associated with a photon or a gluon. Our solution 
to this problem is the introduction of new 4-particles vertices restoring the invariance. A 
complete derivation of their characteristics will be undertaken. Gauge-invariance restoration 
leaves some freedom to add autonomous 4-particle vertices, which are gauge invariant alone. 
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In chapter 6, we shall then apply our approach to the determination of the imaginary 
part of the amplitude for gg ^Sig. This amplitude is the expected main contribution for if) 
and T production at the Tevatron and RHIC at low Pt- This will enable us to calculate the 
various cross sections for these. In turn, we shall introduce the above-mentioned autonomous 
contributions that will allow us to fit the data rather well. A further and complete agreement 
with the data, namely with polarisation, will then be accomplished by a combination of 
our prediction with the colour octet ^Si fragmentation. Indeed, we eventually obtain a 
polarisation parameter roughly constant and compatible with an unpolarised production. 



Chapter 4 

Two-particle description of heavy 
quarkonia 

In this chapter, we shall introduce the quantities required to study processes involving heavy 
quarkonia, such as decay and production mechanisms. In field theory, all the information 
needed can be parameterised by a vertex function, which describes the coupling of the 
bound state to the quarks and contains the information about the size of the bound state, 
the amplitude of probability for given quark configurations and the normalisation of the 
bound-state wave functions. 

The vertex function has close links with the classical non-perturbative wave function 
introduced in quantum mechanics. We shall first consider those, and give their characteristics 
if one thinks that the quarkonia are purely non-relativist ic bound-states. This will show us 
also that the mass of the quark is actually a poorly constrained parameter, which can be 
chosen with some freedom. 

Then, we shall start the description of our approach to build the quarkonia vertex func- 
tions. Along the way, we shall explain why we do not rely on the Bethe-Salpeter Equation [70] 
(BSE), usually used to constrain the properties of vertex functions, for example, to relate 
the mass of the bound state to the mass of its constituent quarks. Beside problems with 
gauge invariance, all predictions coming from BSE are realised within Euclidean space (see 
e.g. [71, 72, 73]) and we shall show that the continuation to Minkowsky space is in general 
problematic. We shall make the choice to remain in Minkowsky space, at variance with 
other phenomenological models (see e.g. [74, 75]). The price to pay for not using BSE will 
be additional uncertainties due to the functional dependence of the vertex. Fortunately, 
most of these can be cancelled if one fixes the normalisation from the study of the leptonic 
decay width. 

We shall thus explain only in few words how in principle from a given choice for the 
vertex functions, we may be able to fix some of the parameters by the use of the BSE. 
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4.1 The charmonium model 

As soon as the J/ip and the ip' were discovered and the community of particle physicists 
started to convince itself that they were composed of a c and a c quark, Applequist and 
Politzer [76] proposed that these states should share some properties with positronium states 
made of and e~. Their statement was mainly motivated by the discovery of asymptotic 
freedom in QCD [77, 78, 79, 80]. 

The identification of the J/ip and '0' as "ort ho- charmonium" ^ was then straightforward. 
Concerning the decay, as the c-quark mass sets a sufficiently high scale compared to Aqcd, 
a formula similar to that for positronium - derived from the QED perturbative expansion 
- was therefore believed to be sufficient, with the natural replacements of a by as and rrie 
by rric. In the latter formula, the use of a coulombic wave function gave the right order of 
magnitude for the leptonic decay width r« [81]. 

As the Bohr radius ((la^mc)"^) of these states was of the order of 1 fm (a typical hadron 
size), they could not be purely coulombic. A strong binding force operating at large distance 
was then added. It resulted in a coulombic potential with a linear rising term preventing 
the freedom of quarks. 

This gave birth to what is referred to as the charmonium model [82] whose potential was 
initially of the form: 



the wave functions "^nimi^) and the binding energies En£ were then calculated, the corre- 
sponding mass spectrum being then built using M„£(cc) = 2mc + Ene- 

4.1.1 Refinements 

Later on, in order to solve the entire spectroscopy of the charmonium and bottomonium 
family, other forms for the potential were proposed: 



V(r) = + kr, 

S r 



(4.1) 



where k is the string tension. 

From the Schrodinger equation 



(4.2) 



• The "Cornell" (COR) potential [84] 



V{r) 



k 



r 



(4.3) 



+ 



2' 



r 



a 



with2 k = 0.52, a = 2.34 GeV"^ ^ ^nd = 1.84 GeV or mb = 5.18 GeV. 



1 S'-wave (n = 0,£ = 0), J=l. 

^ Comparing this form to the one of Eq. (4.1), we see that as ~ 0.38 makes k ~ be 0.52. 
^ ^ here plays the role of the string tension k of Eq. (4.1). 



4.1 The charmonium model 
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The "Buchmiiller & Tye" (BT) potential [85] such that 



V{r) ~ kr {rcrit < r), 
V{r) ^ 



rln(l/A^^^r2) 



l + O 



(4.4) 



where k is again the string tension, Aqcd is the QCD running scale, Vcru is fixed by 
^ - 100 (see [85]) and = 1.48 GeV or = 4.88 GeV. 



QCD, MS 



The "Logarithmic" (LOG) potential [86] such that 



V{r) = -0.6635 GeV + (0.733 GeV) log(r ■ 1 GeV 



(4.5) 



with rric = 1.5 GeV or rub = 4.906 GeV. 



the "Power-law" (POW) potential [87] 



V{r) = -8.064 GeV + (6.898 GeV)(r ■ 1 GeV ) 



0.1 



(4.6) 



with rric = 1.8 GeV or rrib = 5.174 GeV. 



For illustrative purpose, we show in Fig. 4.1 some wave functions obtained by Kope- 
liovich et al. [83] from Eqs. (4.2-4.6). As one may see, even though the values of the charm 
quark mass, rric-, differ substantially from one potential to another, the wave functions do 
not differ much for large r. On the other hand, they differ by 30 % at the origin, value which 
is a key input in the static approximation. Usually, it is extracted from the leptonic decay 
width but one is free to rely on these determinations from the resolution of the Schrodinger 
equation. 

Please note also the presence of a node in the wave function for the first radial excitation 
[ip') as should be from general theorems of quantum mechanics. This node at around 0.4 fm 
produces the suppression of photoproduction of iIj' compared to J / ip [88] . 
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Figure 4.1: The radial wave function \&n^(r) for the J/ip (IS) and ^p' {2S) states obtained by 
Kopeliovich et al. [83] for four potential: BT, COR, LOG and POW. 



4.2 Phenomenological approach to the vertex function 
4.2.1 Wave function vs. vertex function 

We have seen in the previous section how one can calculate the wave function for the quarko- 
nia. One thing is however lacking, that is how to include this information in a calculation 
involving Feynman diagrams. The usual procedure can be summarised as follows: 

• the quarks are to be on-shell within the quarkonia and therefore have to be produced 
on-shell; 

• the Schrodinger wave function provides with the binding-probability amplitude ; 

• the quantum numbers of the bound state come uniquely from the quarks configuration. 
The parity P and charge conjugation C are therefore given by P = {—)^~^^ and C = 
{—)^~^^ where L and S are given by the corresponding eigenvalues of the Schrodinger 
equation. 

From this, it follows that the amplitude for producing a bound state of quantum numbers 
J, L, S is 

M,ot = ^^^l-^^^) E / 7^-^A,A(p1^tf (P)*, (4.7) 



4.2 Phenomenological approach to the vertex function 
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with p, the relative momentum between the two quarks, A, A, the spin projections of the 
quarks along the 2;-axis (quarkonium direction), Mx\{p), the amplitude for producing the 
heavy-quark pair and '^^"^''{p}, the wave function of the bound state. 
The important point here is that the latter reads: 



LzSz 



* (p) — v{ 



(4.8) 



where "^iL^ip) is the non-relativistic wave function, v and u are the spinors for the two 
on-shell heavy quarks - as stated in the assumptions - and Vss, is a spin-projection operator 
whose explicit expression is [89] 



2 ^ 1 1 

'Pss,{PuP2) = ^ Xl^o^i' 9^2|^^2)f(p2, A2)n(pi,Ai) 



where n^^^ 



2m ^'2 '2 

Ai,A2 

1 '][)2 — m - 7° + 1 + m 



2m ^p2fi + m"^^" 2\p2 ^P\f) + vn 
-75 {S = 0) 

re;{s^) {S = l) 



(4.9) 
(4.10) 



The expressions in Eq. (4.10) are what we shall use for our vertex, for respectively a 
pseudoscalar particle and a vector particle. But we are going further by lifting the constraint 
that the quark be on their mass shell. This generalises structureless vertices of the usual 
pQCD to bound states. 



4.2.2 Our choice for the vertex function 

It has been shown elsewhere (see e.g. [72]) that other structures than 7^ Dirac matrix for 
vector bound state are suppressed by one order of magnitude in the case of light mesons, 
and that this suppression is increasing for the meson </> (ss). We have therefore reasons to 
believe that their effect are even more negligible in heavy quarkonia. 

As a consequence, the following phenomenological Ansatz for the vector-meson vertex, 
inspired by Eq. (4.10), is likely to be sufficient for our purposes. It reads: 



V,{p,P) = r(p,P)7^, (4.11) 

with P the total momentum of the bound state, P = pi — p2, and p the relative one, 
P = {Pi +P2)/2 as drawn in Fig. 4.2. This Ansatz amounts to multiplying the point vertex 
(corresponding to a structureless particle) by a function, T{p,P). 
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Figure 4.2: Phenomenological vertex obtained by multiplying a pointvertex, representing a struc- 
tureless particle, by a vertex function (or form factor). 

The function F, which we shall call vertex function, can be chosen in different ways. Two 
simple Ansatze are commonly used in the context of Bethe-Salpeter Equation. We shall 
consider both. They correspond to two extreme choices at large distances: a dipolar form 
which decreases gently with its argument, and a gaussian form: 

A^ 

ro{p,P) = - -T-^, (4.12) 

and 

ro{p,P) = Ne^, (4.13) 

both with a free size parameter A. As T{p, P) in principle depends on p and P, we may shift 
the variable and use — ^^^2 instead of p"^ which has the advantage of reducing to — 
in the rest frame of the bound state. We shall refer to it as the vertex function with shifted 
argument. In the latter cases, we have - in the rest frame -, 

A^ 

and 

T{p,P) = Ne^. (4.15) 

Excited states 

As one knows, the number of nodes in the wave function, in whatever space, increases with 
the principal quantum number n. This simple feature can be used to differentiate between 
IS" and 2S states. 

We thus simply have to determine the position of the node of the wave function in 
momentum space. To what concerns the vertex function, working in the meson rest frame, 

the node comes through a prefactor, 1 which multiplies the vertex function for the 

IS state. Explicitly, this gives, 

P) =N'(l- M_] (4.16) 



anode J (1 + ^) 



4.3 Normalising: the leptonic decay width 
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and 

T2s{p,P) = N'(l- e^. (4.17) 

In order to determine the node position in momentum space, we can use two methods. 
The first is to fix anode from its known value in position space, e.g. from potential studies, 
and to Fourier-transform the vertex function. In the case of a gaussian form, this can be 
carried out analytically. 

The second method is to impose the following relation^ between the IS* and 25* vertex 
functions: 

/" Ipj'dlple^ - = 0. (4.18) 

The two methods give compatible results. The introduction of a second node to describe the 
35* states is similar. 



4.3 Normalising: the leptonic decay width 

To fix the normalisation parameter A^, we shall use the leptonic decay width. Besides, as 
the calculation is similar to the one for the production processes, we shall give a detailed 
explanation of this calculation as a pedagogical introduction to the production calculations 
which are the core of this work (see chapter 6). 

The width in term of the decay amplitude M., is given by 

where d2{PS) is the two-particle phase space [90]. 

The amplitude is obtained as usual via Feynman rules, for which we shall use our vertex 
function at the meson-quark-antiquark vertex. At leading order, the square of the amplitude 
is obtained from the cut-diagram drawn in Fig. 4.3. 

In terms of the sub-amplitudes A'^'^ , B^'^ and C'^'^ defined in Fig. 4.3, we have^ : 

I \M\'d,{PS) = 1a,,,A- B^^'{^)C^'^\ (4.20) 

where the factor A^i, = {g^^y — ^j^) = '^i£i,fj.£t,fi' results from the sum over polarisations 
of the meson and the factor | accounts for the averaging on these initial polarisations. 

^ inspired from the orthogonality between the IS* and 25' wave functions. 

^ In the Feynman gauge. The calculation can be shown to be gauge-invariant, though. 
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Figure 4.3: Feynman diagram for '^S'l 



4.3.1 First sub-amplitude 

From the Feynman rules and using the vertex functions discussed above, we have (see Fig. 4.4) 



zA^^ = -3e 



Q 



-3e 



Q 



(2^ 



Tr (zr(p,P)7'^ 



i{]f: - \lf + m) 



-(ze7 



(k — — m? + is {k + — m? + ie 



r(p,p) 



{{k - 1)2 - m2 + ie)){{k + 1)2 - m2 + ie)) ' 



(4.21) 



eg is the heavy-quark charge, —1 comes for the fermionic loop, 3 is the colour factor. 




Figure 4.4: Feynman diagram for '^Si 7. 

Since the relative momentum is k, we shall have in the bound-state rest frame, for the 
dipolar form of the vertex function as a function of k , 

nP,P) = - (4.22) 



and for the gaussian form. 



r(p, P) = Ne^ 



(4.23) 



4.3 Normalising: the leptonic decay width 
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The expressions for the vertex functions of excited states immediately follow from the above 
discussion and are not written down. 

We can easily motivate our choice of the argument of the vertex function, beside its 
simple relation with that of wave functions. Its main virtue is to regularise the integration 
in all spatial directions. In the meson-rest frame, the integral can be done by standard 
residue techniques, and the remaining k integral is guaranteed to converge. If one uses /c^ 
instead of as an argument, it is possible to show that along the light cone, one obtains 
logarithmic divergences, which would presumably need to be renormalised. Besides, the k^ 
integral then also becomes dependent on the singularity structure of the vertex function. 
These two reasons make us prefer the shifted vertex in our calculation. 

A notable simplification can be obtained by guessing the tensorial form of A'^'^. Indeed, 
current conservation (gauge invariance) for the photon can be expressed as: 



A'^'P^ = 0. (4.24) 

It is equivalent to 

zA'^'' = zF{g^'' - = iFA^\ (4.25) 

the coefficient F being since A^^ = A^^g^^ = 3F. This assumption can be easily verified 
in the bound-state rest frame P = (M, 0, 0, 0) for which Eq. (4.24) reduces to A^^M = =^ 
j^fio _ g shall check this at the end of this calculation. 

It is therefore sufficient to compute the following quantity, where we set |A;|2 = K^: 



POO 

iA''^ = -3eQ A7rK^dKT{-K^)x 
Jo 



dkp ^-2{k^-^)+Am^) 

(27r)4 ((A; - f )2 - m2 + ie){ik + f )2 - + ie)) ' 



Integration on ko 

Let us first integrate on ko by residues. Defining E = ^K'^ + m2, we determine easily the 
position of the pole in /cq (still in the bound-state rest frame) 

P M2 

{k + — )^ -m^ ^%e = kl^ k^M + ^ - n? ^ ie = 

(ko + y)2 -E' + te = {{ko + Y^+E- le'm, + ^)-E + le') = (4.27) 

-M -M 
(ko - (^ + E + te'))iko - (^ + E- te% 
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and 

P 

{k - — f -m'^ + ie = kl- k^M + ^ K'^ - + ie = 

{h - -E' + ie = {{ko -f) + E- le'm, - f) - E + is') = (4.28) 

{ko-{f-E + ze')){ko-{f + E-ze')). 

To calculate the integral on k^, we choose the contour as drawn in Fig. 4.5. Two poles 
—E — ^ and —E + ^ are located in the upper half-plane and the two others in the lower 
half-plane. We shall therefore have two residues to consider. The contribution of the contour 
Cr vanishes as R tends to oo. 



-3ec 



K^dKTi-K^W 



dko{-2k^ + 2K^ + ^ + Am 



{ko- 



E + te)){ko - - ^ - te))iko - {f - E + ie)){k^ - (f + E - te)) 



K^dKT(-K^)x 



2m 



-2(-i?-f)2 + M! + 2E2 



2m2 



-E + ff + ^f 



2E^ + 2m2 



[-2E){-M){-2{E + f )) M{-2{E - f ))(-2E) 



-3e, 



Q 



K'^dKT{-K^ 



2 

2m 

Ume 



2EM + 2m2 2EM + 2m^ 



E + 



M 



E 



M 
2 



(4.29) 



Setting the constituent quark mass 

This way of proceeding is correct only if the poles of the upper-plane are on the left and 
those of the lower-plane are on the right^. 

This crossing would first occur between the points —E + ^ and E — Solving for E^ 
we have: 

M „ M _ M 



~E + — = E ^ E = — = 

2 2 2 

M2 



2 _L ^2 



— = K^ + m^ (4.30) 
4 

M2 „ „ 

= >Q. 

4 



If they cross, the integral acquires a discontinuity from the pinch and becomes complex. 
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Figure 4.5: Illustration of the contour chosen to integrate on ko 



The crossing is therefore impossible when M < 2m. Therefore, to get a coherent description 
of all charmonia (resp. bottomonia) below the open charm (resp. beauty) threshold, we shall 
set the quark mass high enough to avoid crossing for all of these. This sets rric to 1.87 GeV 
{mi)) and mt to 5.28 GeV {ms)- Considering the variety of results obtained from potential 
model, this seems to be a sensible choice. 

Integration on K 

Putting it all together, we get 

_ r ,,.K^n-K^) (2K^ + 3m^) 

One is left with the integration on K for which we define the integral / depending on the 
vertex function whose normalisation is pulled off, 

/ is a function of A through the vertex function T{—K'^) and is not in general computable 
analytically. In the following we shall leave it as is and express A^^ as: 

A^"^ = ^^iV/(A, M, m) ^ A^"" = ^NI{A, M, m)A^"'. (4.33) 
Check of the assumption on A^"^ 

As said before, we can test our assumption of gauge invariance by checking that A'^^ = 0. 
For y4*°, this is straightforward as the integrand is odd. This is a bit less trivial for A^^. 
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Let us go back to Eq. (4.21) and write down A^^{K) omitting the integration on K and 
neglecting the constant factors: 



dkc] 



oc 



oc 



{{k - f )2 - m2 + ie){{k + f )2 - m2 + le) 
dk,{kl + E^-^) 



oo{ko-i^-E + te))iko ~{^-E- te)){ko - {f - E + ie)){k, - {f + E - te)) 
E^ + ^ + EM-^ + E^ E^ + ^-EM-^ + E^' 



-2E)(-M)(-2(E + f )) ' M(-2(E-f))(-2E) 



+ 



OC 



-1 1 



0. 



(4.34) 



4.3.2 Second sub-amplitude 

From the Feynman rules, we have 



[le] 



(4.35) 



where a and h are the momenta of the lepton and antilepton (a + 6 = P), and mi their mass. 
As TTii is much smaller than M and m, we shall neglect it. 

To perform the integration on the two-particle phase space, we shall use the following 
well-known relations [90] 

,2 ^2 ia'-h'^Y' 



7rA/f2 ^ / papP 



1 + — + 2- 

M2 M2 



^7rM2 



pa p/3 



(4.36) 



d2{PS)a.h = j d2iPS)a''b^gaf3 = ^ (4 + 2) ""^^ 



24 ' ' 4 

where A = A(l, a^/M^, b^/M^) = A(l, mj/M^, mj/M^) = 1 - 2m'j/M^ ~ 1 if we neglect the 
lepton mass. 

Therefore, we have 



BP''' =Ue? 



24 



M2 



-.pp 



pppp 



(4.37) 



4.3.3 Third sub-amplitude 

We simply have 



(4.38) 



4.3 Normalising: the leptonic decay width 
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4.3.4 Results 

Now that all quantities in Eq. (4.20) are determined, we can combine them. Thence, we 
obtain^: 

(4.39) 

The leptonic decay width eventually reads from Eq. (4.19) : 



Numerical results 

The value of is in practice obtained by replacing Vu by its measured value, eg by y for c 
quarks and ^ for h quarks and, finally, by introducing the value obtained for / for the chosen 
value of A. The formulae to be used for the different quarkonia are shown in Tab. (4.1). 



Meson 


r,,(keV) 


A^ 




5.40 


/(5.4xl0-S)277r5M3^^ 
Y e4/(A,Mj/^,mc)^ 




2.12 


/(2.12xlO-f^)277r^A/^, 
Y e^I{KM^,,m,)'^ 


T{IS) 


1.31 


1 (1.31xlO-6)1087r5M|j^gj 
Y e-iI{A,Mr(lS),mty 


T{2S) 


0.57 


1 (0.57xl0-6)1087r5M3^2<j^ 
Y e-iI(A,Mr(2S),mb)^ 


T(35) 


0.48 


1 (0.48xl0-6)1087r5M3j3gj 
Y e'i/{A,Af-f{3S),m6)2 



Table 4.1: Formulas to be used to determine A^. 



We sketch some plots of A^ for the J/tp and for the T(IS') to show its dependence on 
A and rriQ and the value it actually takes. The values obtained are also given in tables in 
Appendix B.l and will be used in the production calculations presented in the following. 

We have chosen the range of size of the mesons, i.e. A, following BSE studies [91] and 
other phenomenological models [74, 75], where the commonly accepted values for the J/ip 
are from 1 to 2.4 GeV, and for the T(IS') from 3.0 to 6.5 GeV. 

For illustration, we show in Tab. (4.2) several values of A obtained with BSE [91], ac- 
companied with the vertex-function normalisation. 

Recall that the projector A^^ satisfies A^^A^'' = 3 and A^^.A'^"' = A/'. 
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rric (GeV) 


Vertex functions 




Gaussian 


Dipole 






A (GeV) 




A (GeV) 


1.6 


1.59 


2.34 


1.33 


2.87 


1.7 


2.28 


2.08 


1.87 


2.54 


1.87 


3.31 


1.97 


2.76 


2.31 



Table 4.2: Set of values for nic, A and A^ obtained for the J/ip within BSE approach [91]. 

To what concerns the ip', we give the results for A=1.8 GeV, rric = 1.87 GeV and a 
gaussian vertex function (Tab. (4.3)). The constraint for the node position (Eq. (4.18)) gives 
o-node = 1-46 GeV for these values of A and m^. The normahsation A^' = 16.36 will then be 
the value that we shall retain for the numerical applications. Note that the normalisation is 
strongly dependent on the node position; it is even possible to make the decay width vanish, 
with the consequence of obtaining an infinite value for A^. 



a-node (GeV) 


N' 


0.8 

1.46 

2.0 


2.02 

16.36 

4.12 



Table 4.3: Normalisation for the ijj' as a function of a^ode- 




Figure 4.6: Normalisation for a dipolar form for J/i/,' as a function of A. 
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J/\|/ 

11 I , , ^ 

10 - 
9 

8 - 




1 1.2 1.4 1.6 1.8 2 2.2 



A (GeV) 

m.= 1.60GeV m.=1.78GeV 

171^=1.69 GeV mc=1.87GeV 

Figure 4.7: Normalisation for a gaussian form for J/^ as a function of A. 




Figure 4.8: Normalisation for a dipolar form for J/ip as a function of the c-quark mass. 
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Figure 4.9: Normalisation for a gaussian form for J/'ip as a function of the c-quark mass. 




Figure 4.10: Normalisation for a dipolar form for T(15') as a function of A. 



4.3 Normalising: the leptonic decay width 
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T(1S) 

5.5 r , , , 

5 -\„^ 
4.5 
4 - 




3 3.5 4 4.5 5 5.5 6 



A (GeV) 

171(3= 4.80 GeV 171^=5.28 GeV 

171^=5.00 GeV mb=5.50GeV 



Figure 4.11: Normalisation for a gaussian form for T(IS') as a function of A. 




Figure 4.12: Normalisation for a dipolar form for T(IS') as a function of the 6-quark mass. 
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T(1S) 

6 I , , , 




1 -: 



I ' ' ' ' ' ' 1 

4.8 4.9 5 5.1 5.2 5.3 5.4 5.5 

nib (GeV) 

A= 3.00 GeV A=5.20 GeV 

A=4.10GeV A=6.25GeV 

Figure 4.13: Normalisation for a gaussian form for T(15') as a function of the 6-quark mass. 



4.4 Another point of view 

In appendix B.2, we give some details about another method for determining the parameters 
of the vertex function. The latter relies on a Wick rotation of the relative momentum in order 
to work in Euclidean space. This procedure is absolutely common among the community of 
physicists which employ the Bethe-Salpeter equation and the Dyson-Schwinger equation to 
obtain the mass, the normalisation and the decay width of mesons. 

We have been unable to justify this Wick rotation, especially in the case of the gaussian 
vertex function which reaches infinite values when I/cqP oo and |Xm /cqI > ^o|- How- 
ever, as one can see in Tab. (4.2), the results are similar to the ones obtained by our method 
and considering that this more classical method also gives A, m and A^, we have thought 
useful for the reader to sketch some of its results. This can be seen as an a postenon justifi- 
cation for the values of the parameters that we used, such as A = 1.8 GeV for m.c= 1.87 GeV. 



Chapter 5 



Bound states and gauge invariance 



As already mentioned in the previous chapter, in order to tackle processes involving bound- 
states, we shall phenomenologically generalise the Feynman rules for point-like particles to 
bound states by introducing vertex functions or form factors, keeping by hypothesis the same 
Dirac structure as for the corresponding point-like particle. 

However, before using this formalism, one needs to deal with gauge invariance. Indeed, 
if we consider the production of a bound state associated with a photon or a gluon^, this 
may couple with the constituent before the binding but also after, even though the bound 
state is neutral for the considered interaction. In the framework of Bethe-Salpeter equations 
(BSE) [70], one has to go back to the details of the interactions occurring in the bound 
state. We have to take into consideration which interactions occur and which type of possible 
diagram we shall have. This renders this study cumbersome [93, 94]. 

In the approach we follow, which we refer to as "phenomenological" , gauge invariance 
is of course also broken. The vertex functions get different values for different kinematical 
configurations and the different diagrams, contributing at a given order in the expansion 
parameter, have different weights. As one knows, only the sum of the whole set of diagrams 
at a given order provides gauge-invariant results. Therefore gauge invariance is broken^. 

The aim of this chapter is to introduce the reader to our method to restore gauge invari- 
ance. We shall also choose a phenomenological point of view. Indeed our method is based on 
the introduction of effective new 4-particle vertices as a supplement to the 3-particle ones. 
We shall also show what constraints these effective new vertices have to satisfy. 

First we shall present our derivation for a pseudoscalar particle, produced in association 
with a photon. Compared to our final application to heavy quarkonium production {e.g. J/ijj 
and T), this brings two main simplifications to the discussion: the absence of colour factors 
and a single vector particle, namely the photon. Thence we shall give a similar derivation 
for vector- meson production with a gluon. Finally, we shall state that, as there is an uncon- 
strained degree of freedom in this restoration, we are free to add any suitable contribution to 
the amplitude. This will be accompanied with the derivation of the general form for these. 

As said above, we consider the production of a bound state Q associated with a photon 

^ This depends of course on the bound-state content: quark, lepton,. . . 
^ To the exception of the decay process of the previous chapter. 
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or a gluon, the simplest processes are therefore qq Q7 or qq Qg. We shall consider 
these here and this will be used in the calculation of the imaginary part of the amplitude for 
gg Vg as done in the next chapter. 

5.1 The simpler case of qq X07 

We start from the usual amplitude at leading order (LO) for qq ^o7; where Xq stands 
for a pseudoscalar structureless (point) particle and 7 for a photon: 

M = -i9eu{p^Vp^^,^^eliq)u{p) 

'^-"'^ (5.1) 

where the momenta are defined as in Fig. 5.1, e is the QED charge and g is the qqX^ coupling. 
To emphasise the different Dirac structures appearing in it, we can rewrite it as: 

u{p')v^'o^nt^l{q)u{p) = u{p') {a^\i> - i)r + Pi'r + srif + ih' + >^ri') 

(5.2) 

The factors a, . . . , A are functions of scalar products of momenta of the problem and take 
the following values in this case: 



1 

^point — 7 v3 9 

[p — qY — 

point _|_ _ jj^2 



As announced, should we generalise this amplitude in order to describe the production 
of a bound state, the introduction of vertex functions would break gauge invariance. To 
restore the latter, we may add to these two diagrams a new contribution, constructed from 
a new vertex. This new vertex would account then for missing contributions, akin to direct 
interactions between the photon and the bound state, Xq. We shall call it V^. 

In turn, this novel vertex has to satisfy constraints, e.g. to have the same symmetries as 
the sum of the previous contributions (which we shall refer to as the point ones). To enforce 
this symmetry constraint, we shall merely impose it on the partial amplitude built only from 
this effective vertex. 

A first constraint comes from charge conjugation which can be assimilated to the crossing 
of the quark and the antiquark and reads (denoting the transformation T): 

Vip^p', g, P, m) = T{V^{p,p', g, P, m)) = -f'^V^^-p, -p, g, P, -m)7°. (5.4) 



" ; Ppoint — / \n n TnP^ , 

{p - qY - 
' ' - (p' + qY - m2 " ""^2 • 



5.1 The simpler case of qq —>■ Xq'j 
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p-q 



P / \ P 




(a) Direct 



(b) Crossed 



Figure 5.1: Feynman diagrams for the leading order contributions for qq Xq'j. 

Indeed, returning to the usual form for the amphtude (where we omitted the spinors u, the 
polarisation vector e and the other factors), 



point 



{p — qY — Tn^ 



+ 



ip' + qY " T^"^ 



(5.5) 



and using the following well-known relations, (7^)^ = 7^, 7'^7^ = —7^7^^, (7^)^ = 7°7'^7°, we 
get: 



T(\/^) = -f^V^-p', -p, q, P, -m)7 



{—p' — qY — 
V^ip,p\q,P,m). 



(-75)((-J^ + ^) -772)7^^ 



-p + qY — 



(5.6) 



To sum up, the operation T(-) = 7°(-)^(— p', —p-, q, P, —m)'-/'^ should have no effect on the 
sought-after gauge-invariance restoring sub-amplitude: 



V^{p,p\ g, P, m) = -l^V^\-p\ -p, q, P, -m)7". 



(5.7) 



In the following, instead of deriving directly the effective vertex restoring gauge invariance 
in cases where the qqX^ vertex factors are simply the point ones multiplied by a function 
(this procedure is illustrated in Fig. 5.2), we prefer to derive explicit and general constraints 
expressing how gauge invariance would be broken and consequently how to restore it taking 
into account the above symmetry. 
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Figure 5.2: Phenomenological vertex obtained by multiplying a point \/ertex, representing a struc- 
tureless particle, by a vertex function (or form factor). 

These explicit and general constraints will actually relate the parameters «,..., A. Having 
these constraints at hand, it will then be easier to construct our effective contribution from 
the general expression of Eq. (5.2). We shall simply input in Eq. (5.2) a, . . . , A fulfilling the 
obtained relations between them. 

5.1.1 Symmetry: general constraints 

In order to extract such relations for the symmetry of Eq. (5.4), we shall first rewrite the 
amplitude Eq. (5.2) by ordering the Dirac structures to insure immediate identifications of 
the factors of the latter structures. 

=a-f^{^ - + /37^7'^ + 5-i^{i)' + ^)7^ + A7'^7^ 

= - a-i^-i^'i) + 2a7 V - ^7^7''^ - 2a^^q^(3j^Y - SfYl^ 

+ 2 + ^T^'^t' + A7''7' (5.8) 

^^=^ma7'^7^ + 2^7 V - «7^7^^ - 2a7^g'' - (3^^^ - mS^^^ + 2S'j^p^ 

- 57''7^^ + A7^7l 

The equality (*) is obtained by means of Dirac equations for the on-shell quark and antiquark: 

^{p) = mu{p) and u{p')f = mu{p') (or v{—p'){—f) = —mv{—p')). 

By inspection of Eq. (5.8), we see that we only have three independent Dirac structures: 
^Ai^5 7^*7^^. We have however to keep in mind that p, p' and q are independent, and 

that the identification must be true for all their values. 
Schematically, we shall write: 



Structures 


Factors 


7''7^^ 


2ap>' - 2aq^ + 2 V 
ma — [i — m6 + A 
—a — 5 



Now, we shall apply the transformation T(-) = 7°(-)^(— p', —p, q, V, — m)7° = 7°(-)i^27' 
to the different structures and their coefficients. For the structures, we have 



5.1 The simpler case of qq —>■ Xq'j 
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7^'^ - I'iYlWl' = 7%°^7°(7')W°7^°7° = ^7'^7' = 7^7^^ + 2g^7', 
and for the coefficients^: 

niai — Pi — m5i + Ai — —ma2 — P2 + '^^2 + A2 

— «! — ^1 ^ — Q!2 — (^2 

Finally, we get, 



Structures 


Factors 1 


Factors 2 


7^ 


2aip^ - 2aig^ + 28lp'^' 
mai — (3i — m6i + Ai 
—ai — 61 


2a2p''' + 2a2g^ + 2<52P'' - (02 + S2)2q'' 
-ma2 - P2 + m62 + A2 
—a2 — S2 



which gives 4 different identities: 



«! = 62 

Si = a2 

mai — Pi — mSi + Ai = —ma2 — P2 + '^<^2 + A2 
—ai — Si = —a2 — S2 



(5.9) 



(5.10) 



(5.11) 



In short, we have 



ai = S2 




Si = a2 


(5.12) 


-Pi + Ai = -p2 + A2 





These three last equations contain all the constraints in terms of a, . . . , A, coming from 

Eq. (5.4). The ffist thing to notice is that the point case obviously fulfils these conditions. 

1^2 

Indeed, noticing that Pi — > P2, 

1) Spoint,! — P2 ^ Spoint,2 ~ Pi ~ Olpoint,ly 

2) Oipoint,! ~ Pi ^ C>Lpoint,2 ~ P2 ~ Spoint,!-, ^ 

3) — Ppoint,! + Spoint,! = —TTlP^ ^ + mPg ^ 

^ ~Ppoint,2 + ^point,2 = —{~^)P2 ^ + {—Tn)Pi ^ = Xpoint,! ~ Ppoint,!- 

^ The indices 1 and 2 stand for the two sets of variables {p,p', q, V, m) and (— p', ~p, q, V, —m) for which 
these factors are to be evaluated, e.g. ai — a{p,p' , q, V, m) and a2 — a{—p', —p, q, V, — m). 
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The most important result is that we now have a means to deduce in a straightforward way 
whether or not given contributions, for instance the generahsation of these point contributions 
to bound state ones, have the symmetry required. 

5.1.2 Gauge invariance: general constraints 

A4, the amphtude under consideration in Eq. (5.1) involves e*^{q), the polarisation vector of 
the photon. Being a four- vector, the latter possesses four degrees of freedom, one is removed 
by the orthogonality condition: 

e;{q)q>^ = 0, (5.14) 

three remain, whereas we know that photons have only two. Supposing the photon travels 
in the z direction, one possible choice for these two polarisations could be 

= (0, 1, 0, 0) and 5^ = (0, 0, 1, 0) (5.15) 

The third degree of freedom can be expressed by noticing that any new vector obtained 
by a shift in the photon momentum direction q is also valid: 

^''^ = <, + xg^ (5.16) 

A unique definition of the polarisation vectors is then only achieved by imposing a fur- 
ther condition upon x- This condition is equivalent to fixing the gauge (see e.g. [92]), or 
equivalently, two different values of x correspond to two different gauges ^ . 

Gauge invariance of an amplitude Ai involving n external photons is therefore realised 
by the absence of x in its expression. Hence, as x appears solely in association with^ q^, a 
necessary condition for a gauge-invariant result resides in the current-conservation condition: 



V,,...,X' = 0, V^, (5.17) 

defining by M = Vf,,,„^„Ue''^{qi). 

To what concerns A4 of Eq. (5.1), the requested condition reads 

Vq,, = 0. (5.18) 

We shall perform the same ordering as before but for the amplitude contracted with 
g^, the momentum of the photon. As the contraction vanishes, it will be in turn expressed 
in terms of identities for the coefficients a, f3, 6, A. 

^ Let A be the gauge-transformation parameter such that A'^^ = + d'^A{x). The wave function 
transforms hke ■(/;' = e'^-^^^'fip and the polarisation vector hke et^e^''^ = e^'e''^^ + 9^ / A.{x)e^''^dx. We may 
then identify x with / A{x)e^''^dx. 

^ /i is the index of the photon for which we look at gaugc-invariancc or current conservation. 
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= v>^q, =a^\i> - m + f^iH + m' 



- Q!7 # + 2aYP-q - ai q + PTi - 

— ma'j^^ + 2ap.q'j^ — aq'^'j^ + Pj^^ + mS'j^^ 



26p'.qY + 5g 7 - )^ri 



Finally we get, 



Structures 


Factors 


7^ 
7^^ 


2ap.q — aq^ + bq^ + 2(5p'.g 

—ma + /3 + — A 



(5.19) 



And we can write -the two factors (right column of the last table) being identically zero- 



-aPi + = 
/? — A = m(a — S) 



(5.20) 
(5.21) 



We may check that the point case follows these relations. This is straightforward since 

(5.22) 

5.1.3 Restoring gauge invariance for the "phenomenological approach" 



1) — OipointPl + SpointP2 — —Pi ^ Pi + P2 ^ P2 — 0, 

2) P point ^point TTlP^ 171 P2 TTT'{,Oipoint ^ point) 



Having at our disposal the latter constraints, we are now ready to construct an effective 
4-particle vertex that will restore gauge invariance for the cases where a bound-state vertex 
function enters the calculation. 

The values of a,/?, 5, A for our "phenomenological approach", where we introduce a 3- 
particle vertex function F (V3 = F7^), are easily shown to be: 



Fi mFi F2 , mV2 

Pi ^ Pi P2 P2 



(5.23) 



where Fi = F(2p — q) and F2 = T(2p' + q). This amounts to the following modifications 
compared to the point case: 



a = a. 



point 



^ Ti, f3 — jSpoint X Fi, 5 — 5point X A — Xpoint X 



(5.24) 
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We directly see that Eq. (5.21) is automatically satisfied by this choice contrarily to 
Eq. (5.20) for Fi ^ Y2- This is actually true as long as /5 = ma and A = mb. We have thus 
expressed through one single relation (Eq. (5.20)) involving one single Dirac structure (7^) 
the breaking of gauge invariance. 

The fact that Eq. (5.21) is satisfied shows us that either the effective 4-particle vertex 
should not contribute to this equation - no term in 7^7^ at all (before contraction with q'^) 
or that its contribution should be transverse, i.e. should vanish when contracted with q^. 

Concerning the crossing symmetry, let us first analyse the constraint ai = 62'. the sym- 
metry of the vertex function, Fi — > F2, and Opomt,! = ^point,2 imply that the latter also holds 
in the "phenomenological approach" . The same is true for the other symmetry relations as 
can easily be seen. 

Let us turn back to the derivation of the new vertex. For simplicity, we shall first introduce 
terms that explicitly do not contribute to Eq. (5.21). On the other hand, Eq. (5.20) is not 
satisfied by our approach and corrections from the 4-particle vertex should then appear in 
this equation, i.e. with a 7^ Dirac matrix for sole Dirac structure. This implies that it 
should be of the general form k'^'j^. The vertex should as well satisfy the symmetry relation 
Vt{p,p', g, V, m) = Tiy:) = 7°^i(-p', -p, g, V, -771)7°. 



First solutions 

Setting k'^ to {p+p')^ in k^'-j^ , defining 7/ as a multiplicative factor, we can take 27/(p+p')'*7^ 
which can be rewritten as: 

2ri{p + p'Y-i^ = 7/(2p''^7^ - 7717^^7^ + 2p^7^ - 7^7^777) = 7/(7^/7^ + 7^;^7''). (5.25) 

As we have seen, under the transformation T, 7^ changes to —7^ and rii{p + p') will give 
—V2{p' +p), therefore, we have to impose that 

7/1 = 7/2. (5.26) 

Concerning gauge invariance, a similar equation to Eq. (5.20), with this new term added 
and with the values of a and 6 for the "phenomenological approach" , reads 

= -aPi + 5P2 + viP + p')-qi'' = -^Pi + §A + 2ri{p + p').qj\ (5.27) 

-r 1 -r 2 

Solving for rj, this determines the effective vertex: 

Vr = 2r/(p + pTi' = I'ip + P'Y (5.28) 

[p + p').q 

We also see that 7/ satisfies the symmetry relation Eq. (5.26) as Fi transforms into F2. 
The similar vertex 

yA'= f'~l' iHp-pT (5.29) 
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will suit also. Contracted with g^, it indeed gives (Fi — r2)7^ as required to guarantee gauge 

invariance. The symmetry constraint requires t] to be odd in T. This is also verified as 

„ ^ ri-r2 

'/ 2{p~p').q- 

If we consider a virtual photon, ^.e.with ^ 0, we can also choose 

= ^^^^^'^^ (5-30) 
nevertheless V^e{q)fj, = due to Sfj,{q)q^ = 0. 
Singularity structure 

Now that we have seen that different choices are possible, let us turn to another issue, i.e. the 
singularity structure of the 4-particle vertex. Indeed the above mentioned vertices have 
all different denominators, with a different singularity structure, whereas we shall demand 
that the gauge-invariance restoring contribution does not introduce new singularities in the 
problem. The easiest way to insure this is to limit our choice to vertices with the same 
denominators as the ones already in the problem. 

Of course any combination of the previous solutions is also correct. Thus any structure 
that would not contribute to Eq. (5.20) and Eq. (5.21) and with the correct symmetry can 
be added to them. This freedom will be discussed thoroughly for the qq — > Vog case (see 
section 5.3). 

Since all structures similar to will satisfy the gauge-invariance constraint, it is 

sufficient to concentrate on its symmetry. Since the odd character of (Fi — r2) is compensated 
by 7^, it is obvious that any combination of the vectors of the problem will yield a correct 
k (since k is both in the numerator and in the denominator). 

If we choose k^ = p^{p'.q) +p'^{p.q), we get 

Except for the term g^, the denominators of this solution are equivalent to the propagators 
of the problem. For g^ = 0, they do not bring any new singularities in the problem. But we 
look for a more general case. . . 

A slight different choice for k^ gives: 



_ {2p - qY{2p'.q + g^) + {2p' + g)^(2p.g - g^) 5 
^ ^ ' ^' 2(2p'.g + g2)(2p.g-g2) ^ 

,^ ^ , ( {-2p + qY ^ {2p' + qY \ , 
= (^^-^^) i 2P, ^^IfT)^ 



(5.32) 



The great advantage of this vertex, contrarily to Eq. (5.31), is that it does not introduce 
new singularities in the problem for any value of g^. It can therefore be used in model 
calculation for pion structure functions where the probing photon is off-shell [95, 96]. 
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5.2 Derivation of an effective vertex for qq V^g 



We shall here deal with a vector particle and with a gluon field. We start from the amplitude 
at LO which we rewrite in a simplified form comparable to Eq. (5.1): 



M = -ijj.r»ii_jfi(p2V""'e*(«K(P)t.(p) 



{p — qY — m? 



(5.33) 



where the momenta are defined as in Fig. 5.3, gs is the QCD coupling, g the qqV'^ coupling, 
T" the colour matrices and the colour factor for a singlet state. To emphasise the different 
structures^ appearing in it, we can rewrite it as : 



(5.34) 



Note that we have added a new structure g'^'^ to get a more general discussion in the following. 




p- q 




(a) Direct 




(b) Crossed 



Figure 5.3: Feynman diagrams for the leading order contribution of qq — > Vog. 



The factors a, . . . , k, then take the following values in the QCD perturbative case for 
point-like particles: 



^ Note that these are not necessarily independent. 
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1 



a. 



Pi 1 Ppoint — / NO 2 — 'rnPi , 



point I \n 9 — -'15 y point i no o 

\P — qY — \j) — qY — m'^ 

^point 



with the short notations for the propagator denominators (p^ = m^,p'^ = m?^ q^ = 0), 

1—1 11 

Pi^ = 1 9 = i\ — and P^^ = — ^ = -— . (5.36) 

[p — qY — ni'^ ^P-q [p + q) — fn 2p .q 

Similarly to the pseudoscalar case, we work out a first constraint from the quark crossing 
symmetry, which reads: 

V^''{p,p, g, P, m) = -7V^^"^(-p', -p, g, P, -m)7°. (5.37) 

Indeed, returning to the usual form for the amplitude (where we omitted the spinors u, the 
polarisation vectors e and the other factors), 

{p — qY — m? [p' + qY — rn? ' 

and using the following well-known relations , '^^'^'^ = —'^'^'^^ + 2g^'^, (7^)^ = 7°7'^7'', we get: 

- - nYjY + ^) - nT'h^ 

{—p' — qY — Tn^ {—p + qy — Tn^ 

-V^'{p,p',q,P,m). 

(5.39) 



riy^"") = -f^V^'^^i-p', -p, g, P, -m)7° 



Note the change of sign compared to the pseudoscalar case (Eq. (5.4)). 
5.2.1 Symmetry: general constraints 

As for the pseudoscalar case, we shall now extract from Eq. (5.37) relations between the 
parameters a, . . . ,k. Using the relations 'y'^^'y^ = 2'y'^q'^ — 7^7^^^ = 27^g^ — 2g^^^^ + 7''7^^, 

V^"" = (x-i\i> - + /37"7^ + b-i^{il + ^)7" + A7^7'^ + Kg^"" 

= 2aY{p - qY + oi-i^Yiji -i)- '^o^g^'^iji -i)- h^Y + '^Pg"" + A7''7^ 

+ 2 ^^7^^ - 5j(^'7'^7" + 257'^g" - 57'^7'^^ + K^?''" (5.40) 

= 2^7'' (p - qY + ma7^7'^ - 07^7'^^ - 2ma^'^^ + 2a^'^^^ - (S-i^'Y + 2/3^'''' 
+ A7''7'' + 25p'^7'' - m57^7'' + 25-i^q'' - b-^^-i^ + Kg^"" . 
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The equality (*) is obtained by means of Dirac equations for the on-shell quark and antiquark: 
'f)u{j)) = mu{p) and u{p')'^' = mu{p'). 

By inspection of Eq. (5.40), we see that we have to deal with six Dirac structures: 7^, •y", 
1^1^ 1 1^1'^ ii 9^'^ ^-iid g^'^^'i we also have to keep in mind that p, p' and q are independent, 
and that the identification is expected to be true for all their values. 

Schematically, to help the identification, we shall write: 



Structures 


Factors 




2ap>' - 2aq^ + 2bp'^ 




25q'' 




ma — [i — m6 + A 




—a — 5 




-2ma + 2(3 + K 


9^^ 


2a 



Now, we shall apply the transformation — T(-) : 
different structures and their coefficients. We obtain 



{.~'P' ^ ~'PA^ to the 



7 



-7°(7^)V = -7^ 

-I'wh' = -7^ 



9' 



-9'\ 



(5.41) 



and for the coefficients^: 



2aip^ - 2aig^ + 2(5ip'^ - 


-2a2p'^ - 2a2q^ - 


252P 


25ig^ - 


- 252g^ 




mai — Pi — m6i + Ai - 


-ma2 -132 + ni52 




— tti — 61 - 


^ -a2 - 62, 




—2mai + 2[3i + ki - 


2ma2 + 2/32 + '^2, 




2ai - 


2a2. 





(5.42) 



Finally, we get. 



Structures 


Factors 1 


Factors 2 




2aip^ - 2axqi' + 25ip'^ 


2a2p'^ - 2S2q^ + 2^2^'' 




2biq- 


2a2q'' 


YY 


ma\ — Pi — m5i + Ai 


-ma2 - P2 + 1^62 + A2 




—ai — 61 


—a2 — 62 




—2mai + 2(3i + ki 


—2m52 — 2X2 — K2 




2ai 


252 



^ The minus sign of the transformation — 7'^(')i^27'^ already taken into account in the transformation 
laws for the structures. 
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which gives the following 8 identities relating out of which 5 are independent: 

2ai = 2(52, 
-2ai = -2,52, 



25i = 2^2, 
25i = 2^2, 

mai — Pi — m6i + Ai = —ma2 — P2 + ^^2 + A2, 
—ai — 5i = —a2 — 62, 
—2mai + 2(3i + /ti = —2m52 — 2A2 — 1^2 
2ai = 252. 



(5.43) 



These relations reduce to 4 conditions: 




Let us check that the point case fulfils these conditions. 

1) ^point,! ~ P2 ^ ^point,2 ~ ^1 ~ Olpoint,!: 

2) Olpoint,! ~ Pi ^ Olpoint,2 ~ P2 ~ ^point,!i 

3) - Ppoint,! + Apomt,i = -mPf ^ + UlP^^ (5.45) 

=^ ~Ppoint,2 + ^point,2 = ~{~'^)P2 + {—T^)Pi = ^point,! ~ Ppoint,!, 

4) — {Upoint,! + f^point,2) = = 2(Ppoint,! + ^point,2) = mP^ ^ + (—m)P^ ^. 

We have as in the pseudoscalar case a means to deduce in a straightforward way whether or 
not given contributions, for instance the generalisation of these point contributions to bound 
state ones, have the symmetry required. 



5.2.2 Gauge invariance: general constraints 

If we now turn to gauge invariance, we shall perform the same ordering but for the amplitude 
V'^'^ contracted with g^, the momentum of the gluon. Current conservation implies that the 
contraction should then vanish. This annulment is in turn expressed in terms of identities 
for the coefficients a, /?, 6, A, k. 
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= v^'^q, =a^\i> - m + ^''i + m' + + Hi" + ^^q" 
= - cyiH^ + 2aYp-q - aYq^ + PiH - ¥iY 

+ 2^47' + Sq'^Y + 'iXq" - XY^ + f^q" (5.46) 
= — ma'-y'^^ + 2ap.q'y'^ — aq^'-y" + /57'^^ — 2m6q'' + ■mS'j''^ 
+ 26p'.qY + Sq^Y + 2Ag^ - A7''^ + nq" 

We are left with 3 independent structures, 7^^, 7^^ and g"^ with the following coefficients: 



Structures 


Factors 


Y 


2ap.q - aq'^ + 26p' .q + 6q^ 


Yi 


—ma + (3 + m5 — A 




-2m5 + 2X + K 



This gives 



-aPi + 5P2 = 


(5.47) 


P — \ = m{a — 6) 


(5.48) 


K = 2{m5 — A) 


(5.49) 



The point case follows these relations. Indeed, 



1) — apointPl + SpointP2 — ^-^1 ""^-^1 + -^2 ^-^2 — 

2) Ppoint ~ Xpoint = ^ — ^ = TTl^apoint — ^point) (5.50) 

3) 2{m6point - \point) = 2{mP2^ - mP^^) = = KpoM 

5.2.3 Restoring gauge invariance for the "phenomenological approach" 

Let us construct an effective 4-particle vertex restoring gauge invariance, broken by the 
"phenomenological approach". If we introduce a 3-particle vertex function F {V^ = T'-f"), 
the values of a, P,6, X are : 

Ti n "^ri F2 mF2 . - 

a = f3 = -—,6 = —,X = -—,K = 0. (5.51) 

-r 1 ^1 -'2 -'2 

where Fi = F((2p - g)^) and F2 = F((2p' + g)^). 

a = CXpoint X Tl; /3 = Ppoint ^ Fi, 5 = ^point ^ 1^2; A = Xpoint ^ F2. (5.52) 
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Prom the point of view of gauge invariance, we directly see that Eq. (5.48) and Eq. (5.49) 
are automatically satisfied by our phenomenological 3-point vertex contrarily to Eq. (5.47) 
at least as long as Fi ^ T2- Concerning the symmetry under — T, it can easily be seen that 
it is also satisfied. 

The single relation Eq. (5.47) derived from coefficient identification of one single Dirac 
structure (7*^) will be also our starting point to restore gauge invariance in this case. However 
we shall proceed bearing in mind not to add any contribution to Eq. (5.48) and Eq. (5.49), 
i.e. no term in 7'^7'^ at all (before contraction with g^), or only a transverse one (vanishing 
due to the contraction). 



First solution 

Starting from contribution whose Dirac structure is 7*^, we can choose the following one 
(with the multiplicative factor 77): 

27]{p + p'Y-i'' = r]{2p'^'Y - m-f^'Y + 2p''Y - 7^7'^m) = 7]{YfY + Yff^)- (5.53) 

As we have seen, under the transformation — T(-), 7"^ changes to — 7"^ and ?7i(p + p') 
changes* to —ri2{p' + p), therefore, we have to impose that 

Vi = V2- (5.54) 

Concerning gauge invariance, the introduction of this new contribution adds a new term 
to Eq. (5.47): 

= -aPi + 5P2 + 2r/(p + p').q = -^Pi + + 2r]{p + p').q. (5.55) 

Pi P2 

Solving for rj, we have: 

vr = 2v{p + p'YY = Yip + V'Y. (5.56) 

[p + p).q 

T] being i^^^^i^- We in turn see that r] satisfies the symmetry relation Eq. (5.54) as Fi 
transforms into F2. 



The similar vertex 



vr= f' ^' I'ip-p'r (5.57) 

also satisfies all the criteria. Contracted with g'^, it indeed gives (Fi — F2)7'^ as required to 
obtain gauge invariance. The symmetry constraint would then require rj to be odd. This is 
also verified as 77 = tJ-^^^tt-. 

As in the pseudoscalar case, if we consider a virtual gluon, i.e. with 7^ 0, we might 
have also chosen 

Fi^Fs 
q' 



vr = -^^Yq^, (5.5J 



which does not contribute due to e^{q)q^ = 0. 

® The minus sign of —T{-) is considered only the transformation of the structure. 
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Generalising. . . 

Similarly to the pseudoscalar case, where ^^^q^ l^k^ was a possible gauge-invariance restoring 
vertex for any /c, we have here ^^^^^ I'^k^ ■ 
If we choose = p^{p'.q) + p'^{p.q), we get 



=^^'-^'^ 2ip'.q)ip.q) ^ =^"^-"n^'"^'" 



(5.59) 



In this case, this choice contains the same propagators as the point case: it cannot 
introduce new singularities in the problem. As before, there exist also choices of k to make 
appear the denominator of the problem when 7^ 0, i.e. for an off-shell gluon. If we choose 
/c^ = (2p - qY{2p'.q + q^) + {2p' + qY{2p.q - q"^), we get 

^ .p _ P . (2p - g)^(2p^g + q') + (2p^ + g)^(2p.g - g^) , 
^ ^ ' 2(2j9'.g + g2)(2p.g-g2) ^ 

(Fi-r^) /(2p-g)^ (2p- + gr \ p./ (2p-g)^ (2p- + gr \ , 



2p.g - g2 2p'.q + q^ J ' ' ^ V 2Pi 2Pi 



(5.60) 



suitable for an off-shell gluon. 

The choice of Eq. (5.59) associated with the contributions derived from our phenomeno- 
logical 3-point vertex will be what we shall refer to as the Gauge-Invariant Phenomenological 
Approach (GIPA). The latter will be considered as our benchmark for the study of non-static 
effects and will be studied in details when applied to the calculation of the quarkonium cross 
section at the Tevatron and at RHIC in chapter 6. 



5.3 Autonomous contributions for qq V^g 

Following the discussion for qq V^g, we may generalise the vertex by adding new contri- 
butions that are gauge-invariant. This is in fact equivalent to what is implicitly done when 
passing from one gauge-invariance restoring contribution to another . 

For this reason we shall derive here a generic form for a gauge-invariant contribution that 
can be added to the amplitude for the bound state. We shall nevertheless restrain ourselves 
to certain class of contributions that satisfy reasonable constraints, namely 

• They should vanish for a structureless bound state, that is when Fi = 

• The colour factors, as well as the coupling constant from QCD, are to be the same as 
in the corresponding perturbative contribution with the same external particles. 

• Once again, no new singularities can be brought into the problem. 
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• The new contributions should have the same symmetries under quark crossing as the 
structureless ones. 

As we are dealing with two vectors particles and two of spin 1/2, we have have to have 
the following generic form for the amplitude 



u{p')V^''el{q)el{P)u{p\ (5.61) 

(g) is the gluon helicity vector and e*{P) that of the vector meson. 
The only independent tensorial structures that we can introduce are made of Dirac ma- 
trices and of the momenta of the problem. From momentum conservation, we might choose: 
P, q and p + p' (recall that p — p' = q + P). Therefore, V^^ could be constructed from a 
product like the following 



1^/^^ = (Fi - Ts 



\{p + pTJ 



r 



(5.62) 



\{p + pT/ 



and with g^"^, but this simple combination does not insure gauge- invariance. 



The latter can be insured by use of the projector 



9' 



. For (yf^^, this gives 



yi'^ = (Fi - T2)ig^'' - ^). As the index of the bound-state u is not affected by the gauge 
constraints, let us consider the factors in Eq. (5.62) with index fi first: 























r') 




\{p 


+ P' 





\{p + p' 



Em \ 



/ q.P 



q.P 



(5.63) 



We can drop q^ as the index is the one of the gluon and e{q).q = 0. We do the same for 
P" {e{P).P = 0) and finally for P^. Indeed because the gluon is real and thus transverse, 
we have e{q).k = for all time-like vectors, k, like q + P -which reduce to (v/i, 0, 0, 0) in the 
CM frame-. e{q).{q + P) = and e{q).q = obviously give e{q).P = 0. We therefore have 







/ r 


\ 





















\ip+p'] 




\ip+p'] 


7 



This drives us to the generic contribution : 

= (r, _ F2)(ai7'' + a2q'' + a^{p + p'^l^il^ + Hp + P'Y) + io^" (5-65) 
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Please remember that the latter expression is without terms that explicitly do not con- 
tribute due to contraction with helicity vectors of the bound-state and of the gluon. V^'^ in 
Eq. (5.65) is not meant to verify current conservation as is. We would rather write to be 
complete: 



=(ri -r2)((ai7" + «2P" + a3g'^ + a4(p + p' 



q.P 



q.P 



^9 



)) 



(5.66) 



As long as the coefficients a, /3, 7, k, A, ^ are not singular anywhere, the second condition 
that we required is satisfied. To what concerns the crossing symmetry, first let us recall that 
Fi — r2 changes signs since Fi — > r2, the rest of the expression must then also change sign. 
Under — T(-), we have: 



+p'y ^ +i^{,p+p'Y {p+p'Yip+p') 



Yip + p"" 
pt^qf^ 



-{p + p'Y{p + p') 



^9 



9 



+p^(p 



p. 



(5.67) 



the other combinations are trivially obtained from these. As we required them to change 
sign, to cancel that of Fi — F2, we have to impose that '-/'^{p + p')" as well as q'^{p + p'Y have 
coefficients changing sign contrary to {p + p')'^{p + p'Y , whose coefficient can be a constant. 

On the other hand, we have to rule out 7'^7^, the term —2g^'^ cannot be cancelled by 
the terms proportional to g^'^ . Only more complicated combination starting from a more 
general product: 



= (Fi - F2 



pfj. 
q^ 



\ 



( 



\{p + p'YJ 



\ 



( 



r 

pv 

q^ 



\ 



(5.68) 



\(p + p')7 



would allow to have terms including 7'^7''. 

For a first approach, we shall retain here only vertices with a constant multiplicative 
factor. This enables us to drop some terms in Eq. (5.65) which eventually becomes: 



V^- = (Fi - F2) (a ^q-" + (3{p + p'Yip + p'Y + e 9n- 



(5.69) 



5.3.1 Our Feynman rules 

Considering the process qq Vog in the GIPA, where the qqVo vertex is V{p1^i, P)'~^^ {prei 
is the relative momentum of the quark and the antiquark and P the momentum of Vq), the 
choice given by Eq. (5.59) is written explicitly in the first line of Tab. (5.1). 
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To what concerns autonomous contributions, we limit ourselves to the ones where the 
free multiplicative factor is a constant. The Feynman rules for these three possible choices 
derived in Eq. (5.69) are also given in Tab. (5.1). 



p Q 







-igsT^,{V{2p+q) - r(2p'-,)) 



+ 



r 



■ rpa {r{2p+q)~T{2p'-q)) 



p.q p'.q^ 



r 



= -% gs a (r(2p+,) - Ti2p'-q))Y'q' 



-igsT^,C{r{2p+q)-r{2p'-q))g 



cf.Eq. (5.59) 



First term 
in Eq. (5.69) 



Second term 
in Eq. (5.69) 



Third term 
in Eq. (5.69) 



Table 5.1: Feynman rules for one choice of gauge-invariance restoring vertex (first line) and for 
the three possible choices of autonomous contributions with a constant multiplicative 
factor. 



Chapter 6 



J/^ and T production at RHIC and at 
the Tevatron 

We shall now calculate the inclusive polarised cross sections for pp — *• i/^X (or pp —>■ ipX) and 
pp — i> TX (or pp TX) in the experimental conditions of the Tevatron collider based at 
Fermilab and RHIC at BNL. As already mentioned, the Tevatron is a collider with two beams 
of equal energy, one of protons and one of antiprotons. The energy in center of momentum 
frame amounts^ to 1.8 TeV. There exist two main detectors operating on this collider, CDF 
and D0. We shall here solely consider the results coming from the CDF collaboration as their 
analysis reckon, oppositely to that of D0, on a central magnet which enables to identify the 
charges of the detected particles. As a consequence, their results are slightly more reliable. 
Concerning RHIC, we shall focus on the PHENIX analysis for collisions at -/s = 200 GeV. 

The main study of the CDF analysis, as seen in chapter 2, has been centered on J/tp, 
if)' and the T{nS). For the latter, the direct production has been extracted only for T(IS'). 
The Pt range was from 5 to 20 GeV and the pseudorapidity range from -0.6 to 0.6 for if) 
and from -0.4 to 0.4 T. PHENIX studied only total J/ip production in the central region 
\y\ < 0.35 and in the forward one —2.2 <y< —1.2. The Pt range was much smaller, from 
0.4 GeV to 5.5 GeV. 

In those conditions (very high energy and hadron beams), we have already explained that 
that the main partonic contribution will be gluon fusion. Therefore, one can easily convince 
oneself that a priori the LO sub-process to produce J/ if) and T (generically denoted Q) at 
medium transverse momentum should be gg — > Qg, where the final state gluon is required 
to provide with the desired transverse momentum. 

We first explain the kinematics, underlining the link between the hadronic cross section 
and the partonic one. The latter is then expressed in terms of helicity amplitudes, which will 
enable us to calculate the polarisation parameter apoi measured by CDF (see section 2.5). 



In run I where the data were taken. 
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6.1 Generalities on the kinematics 

The link between the partonic cross sections and the hadronic ones rehes on the following 
general formula: 

E^= [\x,dx2Gi{x,) G2{x2) -^6{s + i + u-M^), (6.1) 
di Jq tt at 

where P = {E, P) is the momentum of the meson in the cm. frame of the colliding hadrons. 
Note on the way that these variables are experimentally accessible. The functions Gi are the 
parton distribution functions, which give the probability that a parton i take a momentum 
fraction Xj off the parent hadron. 

In the cm. frame of the colliding hadrons, introducing the rapidity the transverse 
momentum Pt = {Px, Py) and the transverse energy^ we may decompose 

P as follows (satisfying the on-shell constraint P^ = M^): 



'M2 + P2 + p2 cosh y, P^,Py, y'M2 + P2 + p2 gj^h y) = {Et cosh y, Pt, Et sinh y) , 

(6.2) 

with 

Now, let us express the measure d^P in terms of the variables Pt and y, i.e. d^P = 
dPzdPT = dPzdcj) Pt dPT. As we want to keep Pt, we shall only express the z-component 
in terms of in order to get '^^^ ■ 



dy 

dP. 



Pz = Et sinh y 



ETCOshy = E, (6.4) 



dy 

which enables us to rewrite the integration element as: 

rfV ^ d^a 

dP^ dyPTdPTdf ^ ■ ^ 

As the integrand is symmetric in </>, we obtain the double-differential cross section upon 
Pt and y from Eq. (6.1): 



da 

dydPT Jo 



/ dxidx2Gi{xi)G2{x2)TsPT^5{s + i + u- M"^). (6.6) 
Jo dt 



When there might be an ambiguity, the Mandelstam variables relative to the partonic pro- 
cess are discriminated from the hadronic ones by circumflexes. To simplify the notations. 
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especially when dealing uniquely with partonic processes, we shall sometimes drop the cir- 
cumflexes, as in the next section. 

At this stage, we can perform the integration upon X2 (or xi) using the delta function, 
which gives us the jacobian: 

1 

d(s+i+u-M^) ' (6.7) 

dx2 

To evaluate the latter, we have to express s, t and u in terms of xi, X2, Pt and y. If the 
momenta of the initial gluons are ki and /c2 and that of the flnal ones is (see Fig. 6.1), we 
have: 

s = {k, + k2f = {q + Pf, 



i={kr-Pf = {k2-q)\ 
u = {kr-qf = {k2-Pf. 



(6.8) 



ki 



k, = P 




Figure 6.1: Definition of ki, k2, k-s (P) and k^ (g). 

Deflning the momenta of the two hadrons as pi and p2, we may always choose the axis 
such as to suppress their transverse components (along the x and y axes for deflniteness) . 
For massless initial hadrons, we then have in their cm. frame 



PI = (^,0,0.^), 

p. = (^.0,0.-^). 
Now, by deflnition of the momentum fractions Xi and X2'- 



(6.9) 



{ki + k2 



(xi +a;2,0,0,xi - X2) 



SX1X2, 



(6.10) 



and 



i = (A;i - py 



s 





p^ 

hp 



-2 ki -P 

X'lPl 



-.M^-2x,^{E-P^). 



(6.11; 
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From Eq. (6.2), we have [E — Pz) = -Er(coshy — sinhy) = Etc ^ and 

i= M"^- xie~y^sET. (6.12) 

Similarly, we get 

u = M"^ -X2ey^sET. (6.13) 

Writing down s + t + u — = sxiX2 + — xic""^ ^/sET — X2e'^^/sET and solving for X2, 
we obtain: 

For completeness, solving for xi we have 

^ X2ET^sey-M^ 

^s{^sx2- Ere-y)' ^ ' ' 
We also easily compute the jacobian of the 5 integration: 

1 1 

V~s{^xi - Erev) ' (6-16) 



(s+t+-u-M2) 



dx2 

As X2 is now a function of xi, we are left with a one-dimensional integration which is restricted 
by 

0<a;2<l. (6.17) 
Supposing that'^ Xi > ^x^^ ggt the following condition: 

- y^iEre-y - v^) ' ^ ' 

Therefore, the double differential cross section upon Pt and y takes the following form: 

da 2sPTGi{xi)G2ix2ixi)) da 

dXi r -y:. (6.19) 



dydPx Jx^"^ \fs{\/sxi — Exey) di 

The last step is now to relate the partonic differential cross section ^ to the amplitude 
calculated from our model. To this end, we use the well-known formula: 

§ = T7^^l^l^ (6-20) 
dt IQns'^ 

where is the squared amplitude averaged on the initial spin and colour degrees of 

freedom and summed over the final ones for an unpolarised cross section, or averaged only 
for the colour for polarised cross sections. 

The calculation of the latter amplitude within the framework of our model is the subject 
the following section. 

The conditions t <Q and m < give xi> and X2 > ^^^'-^ ^'^'^ always less restrictive than 

xi > "^Jf- and X2 > " for Pt ^ 0. Furthermore, as long as \y\ is not of the order of In(^) (|j/| = 3.8 
at Pt — 40 GeV), xi > is always less restrictive than Eq. (6.18). 
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6.2 Amplitude for gg ^Sig 

In this section, we shall give our procedure to calculate the imaginary of the amplitude for 
gg — s> ^Sig in our formalism through Cutkovsky rules. We shall consider this amplitude at 
leading order in as- 

First, we further develop the kinematics to avoid any ambiguity and we derive expressions 
for the polarisation vectors of all particles. Then we write down the expressions for the 
helicity amplitudes. The cutting rules impose further constraints on the internal kinematics, 
which we give. We then organise the calculation by classifying the integrals that enter the 
amplitude. 

ki h = P 




k2 ^4 = q 

Figure 6.2: Kinematics of the process gg ^Sig. 



6.2.1 Further kinematics and polarisation vectors 

The momenta of the particles are as shown in Fig. 6.2. The on-shellness conditions are then: 
k\ = ki = ki = 0, fcg = M^. The Mandelstam variables are defined as usual (dropping the 
circumflexes) 



s = {ki + k2f = 2ki.k2 = (^3 + kif = M^ + 2^3.^4 , (6.21) 

t= {ki- ksf = -2ki.k3 + = (^2 - k^f = -2k2.ki , (6.22) 

u={ki- kif = -2ki.ki = {k2 - k^f = -2k2.hi + , (6.23) 

s + t + u = . (6.24) 

and all scalar products, /c^./cj, are functions of the latter: 

ki.k2 = -,ki.k3 = — - — ,ki.k^ = --, (6.25) 

- M , , t , , s- , , 

k2-kz = ,k2.ki = --,k3.ki = — . (6.26) 



In order to compute polarised cross sections, we have to introduce explicitly the polari- 
sation vectors for the vector meson and for the gluons. It is however sufficient to determine 
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for all polarisations the scalar product between these vectors and the other vectors of the 
problem, i.e. ki, . . . , k^. 

Let Ei = e{ki) be the polarisation vector of the different particles considered in the present 
case, we shall perform a Sudakov decomposition of these along the momenta of the problem 
and along a vector er (ct-Ct = — 1) orthogonal to the plane of the scattering. 

A generic helicity vector 5^, where i is the particle label {i = 1, ..4), h the helicity state 
(Ti or T2 for the gluons and Ti, T2 or L for the vector meson), then becomes: 



= a'lki + b'lk2 + dlki + rf^T- (6.27) 

We first impose that each vector be orthogonal to the momentum of the particle it 
describes: 

h.e'l = 0, (6.28) 

for all polarisation h. 

Then we require that vectors of different helicities {h 7^ h') and of one given particle be 
orthogonal to each other: 

e14' = 0. (6.29) 
We in turn impose the correct normalisations, 

{e^f = -1. (6.30) 

The previous conditions are valid for all polarisations. In the case of the transverse 
ones, we shall introduce this transverse property by requiring that only spatial components 
remain. In other words, the scalar product between transverse polarisations and a vector of 
the generic form v = {vq, 0, 0, 0) should vanish. 

In the cm. frame (of the gluons), ki + k2 or k^ + k^ have the latter property since 
/ci + ^2 = ^3 + ^4 = {\/s, 0, 0, 0), and thus the condition for transverse polarisations becomes: 

{ki + k2).ej = {k3 + ki).ej = 0, (6.31) 

Taking into account the conditions of Eq. (6.28), we then derive 

ej.k2 = 0, 
e^.ki = 0, 
£3 .A;4 = 0, 
eth = 0. 

Here, we are free to choose ej'^ to be a vector orthogonal to all the others, namely 

efi = 41 = = = er, (6.33) 
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since this satisfies all conditions cited above. 

Having fixed ej^, we shall determine ej^ imposing sj'^.ej^ = 0. We have: 



Ut S 

s u 



.T,_T,_ h (6.34) 



^ ^ ^' ut t + u ' 
dl^ = df = df = = 0. 

Finally, requiring e^Ke!^ = = e^^ef^, we fully determine : 

Thus, the complete expressions of the polarisation vectors are the following: 

^1 — ^2 — ^3 — — ^T, 



= a/ (tki + uA;2 + sk4) = e?. 



1 / , fs + AP 



Compatibility with the standard choice 

We shall show here that our result can give in a non-covariant form the standard choice for 
instance taken in [50, 51]. Let us do that for the longitudinal meson polarisation vector. We 
have with the notations explained above: 

whereas the standard choice used in [50, 51] is in the meson rest frame {k^ = (M, 0, 0, 0, )): 

4 = (0,0,0,1), (6.38) 
First we rewrite in terms of ks, this gives 
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Choosing the z-axis direction to be the same as — /C4, we have = (/c4,o, 0, 0, — |A;4,2|), which, 
for an on-shell gluon, gives = (/c4,o, 0, 0, — A;4,o)- We are thus left to determine k^^o, which 
is given by 



s = {k3 + kif = + + 2k^.ki = + 2Mkifl k^^ 
Finally, we have 



2M 



1 

M 



((M, 0, 0, 0, ) - (M, 0, 0, -M)) = (0, 0, 0, 1). 



(6.40) 



(6.41) 



6.2.2 Expressions for the amplitude 

As we shall consider only its imaginary part, the amplitude of the process can be diagram- 
matically presented by a product of two scattering amplitudes gg — > qq and qq — > Qg where 
we put the quark lines with momenta i + ki and i — k2 (Fig. 6.3) on their mass shell, so that 



they satisfy the relations: {£ + ki 
propagators by S functions: 

1 1 

{i + kiy-m^ {l-k2Y-m 



m 



-k. 



m , or equivalently we replace their 



2 > -(27r)25+((£ + k,Y - m2)r ((£ - k^y - m^) 



{2'Kf5{2i.{ki + k2))5{i^ + 2i.h - m^) x 



(6.42) 



X e{t + A;^)^(-r + kl) 



where the factor ^ comes from ^mA = |Disc(^), and the discontinuity is obtained with the 
introduction of the above-mentioned 6 functions. 

At leading order in ids), we have three diagrams for gg — > qq and also three for 
qq Qg, in our framework. The six diagrams'^ are drawn in Fig. 6.3. 




(QQ) 



(QQ) 




[III 



fllll 



(QQ) 




[Vi 



Figure 6.3: Different contributions to the imaginary part of the amplitude for qq Qg. 



^ Compared to notations used for the derivation of the gauge-invariance restoring vertices (chapter 5), 
we have p £ + ki, p' ^ £ — k2, q ^ and P ^ k^. 
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These scattering amplitudes have the following explicit forms, where p, g, r and s are 
the helicity states of the particles, 

[O] = {-gs)r'M^ + h) [{h - hYg^'' + + ihYg''^ - {2h + hYg^^] x 

X (XTfc2'2 i-^9sTi)Y'u{i + (6.43) 



X e^^ihYlih), (6.44) 

[I] = {-^9s?T^iTiu{l - h) ' ^^'^1^2^ u{i + k^Yl{k^Yl{k^\ (6.45) 

[II] = {~^g.fTlTt,u{^ - k,) ' ""^^[f^jl'S^lylZ^ u{l + k,Yl{k,Yl{k,l (6.46) 

[III] ^ ^.g;rt^^mi^ - ..) - "^^ ^(^ + x 

X ^'^(^3)^(^4), (6.47) 

2 [7 ^((/+^i-^4) + m) r] 

(i + ki- k^Y - 

X £^'^(^3)^(^4), (6.48) 

[V] = [Vi] + [V2] (6.49) 



[IV] ^ i-^r^s^^Y-WsT^^uii - ..) - J 'U i^ k.) X 



[Vi] = -^ {T, - r,) (-^^.)T,=,^(£ - k2)Yjl^jfz^k^^2^i^ + ki) X 

X £-(^3)5^(^4), (6.50) 

[V,] = -z (Fi - r^) (-z^7.)T^.n(^ - k2)Yjf^^^^^^^2^(^ + ^1) X 

X e^;{ksY':{k,). (6.51) 

Before going further, let us draw the reader's attention to some points concerning these 
amplitudes: 

• Gauge invariance imposes that ([0] + [I] + [II])£:''(A;2)^i = and not simply ([0] + [I] + 
[II])A;f = 0. The same remark holds also for k2. Indeed, using the Dirac equation for 
the on-shell quarks, we have ([0] + [I] + [II])A;5' = ig^^[T'',T^]u{i + ki)f^iu{i - ki) 
which is orthogonal to ev{k2). Their product therefore vanishes. 

• Gauge invariance for the right part is insured by the introduction of the diagrams [Vi] 
and [V2] that come from the choice of the following gauge-invariance restoring vertex 
derived in section 5.2.3, Eq. (5.59) (see also Tab. (5.1)) 

^'\Pi ^\{l + k^-kY)^-m^ {l-k2 + k^Y-'^^) 

(6.52) 
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5. 



• We have split the contribution [V] in two parts as they contain two different propaga- 
tors. 

• We have chosen to omit the ^ factor in the colour factor for the quark-antiquark- 
meson vertex. Of course, the same choice has been made during the normalisation 
procedure (cf. Eq. (4.21)). 

• The contributions of the triple-gluon vertex [0] will vanish when contracted with the 
right part of the diagram due to Furry's theorem. Indeed, a loop with an odd number 
of vertices of vectors particles vanishes. We have nevertheless checked numerically that 
it was the case, also for the gauge-invariance restoring vertex [V]. 

• This set of contribution is what we shall refer to as the Gauge-Invariant Phenomeno- 
logical Approach (GIPA). 

The vertex functions depend on the relative momentum of the quarks at the vertices as 
follows 

T, = T(^i+]^{h-k2 + kS^ , T2 = T(^i + k,-h,}j . (6.53) 

Multiplying the amplitudes from the right with the ones from the left and taking the 
trace in Dirac matrix space and in colour space, we define the following quantities' 

tr^m = - ^?fTr [T'^t'T^] Tr [7^^ (/ + m) 7^((/ - h) + m) 

xr ((/ -h + h) + ^) Yii^ + h) + m)] el{k{)el{k,)e^;{k,)el\k,) 

trf:,V = - 9> [T^T'T'^] Tr [7^^ (/ + m) - h) + ^) 

yiY ((/ + h - h) + m) + h) + el{k,)el{k2)e^;{k,)el\k,) 

^Af-iii = - Tr [T^T»T^] Tr [Y ((/ + ^1 - h) + m) ^ii^ " h) + m) 

((/ -h + h) + m) 7^((/ + h) + ^)] el{k,)el{k,)e^;{k,)el\k,) 

tr^rV = - ^?'Tr [T^T'^T^] Tr ((/ + ^1 - + m) ^{{^ - %) + m) 

X7^ ((/ + - U + m) 7'^((/ + h) + m)\ el{k,)el{k2)e';{k,)el\k,) 

i^Tk = - 9> [T^T'T^] Tr [7^^ (/ + m) ^{{^ - U + m) 

X7''((/ + ^1) + m)] (r + kl)el{h)el{k2)e^;{k,)el\k,) 

trr4 = - glTv [T-T'T'^] Tr [7^^ (/ + m) - h) + ^) 

xY{{^ + h) + m)] (-)(r - k^,)el{k,)el{k,)e^;{k,)el\k,) 

tr?r_V. = - 9> [T'T^T'^] Tr ((/ + ^1 - ^2) + m) ^{{^ - + m) 
xYii^ + ^1) + m)] (r + k'[)elik^)elik2)e;\ks)et\k,) 



(6.54) 

(6.55) 
(6.56) 
(6.57) 
(6.58) 
(6.59) 
(6.60) 



^ Note that the following quantities are defined so that they contain neither the denominator of the 
quarks propagators nor the meson vertex function. 
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tr?r_V. = - 9!Tr [t'^T^T^] Tr [Y ((/ + - h) + m) 7^((/ - h) + m) 
X7''((/ + ^i) +"^)] (-)r - A;2>^(A;i)4(fc2)er(A;3)C(fc4) 

Using the latter quantities, it is now straightforward to rewrite the different contributions 
to the imaginary of gg ^Sig by including the integration over the internal quarks momenta 
constrained by the two 6 functions coming from the cutting rules which is designated by 
dps = ^S{2e.{h + k2))5{P + 2£.ki - m^)e[e + + k^) in the following, 

[I][III](p,g,r,.) = /^P^r.trr;, _ _ 

[I][lV]{p,q,r,s) = / dpsT2trf_!'jy- TT^ 2 

J £^ — [i + ki — k4^)^ — m'^ 

f 11 

|iil|in|(p.,. r,.) = J dp.r.tr;r,„ ^ _ ..^^ _ _ ^ ^^y _ 

|IIl[IV](p,,,r,.) = /<ip.r,tr?r, ^ _\ _ „ _\ _ 

1 ^ (6-62) 
[I] [Vi] (p, g, r, s) = / dpsiTi-T2) trf ^ — — — 

f 11 
[I] [V.] (p, q,r,s) = Jdps{r,-T2) tr?- _ ^ _ 

[II] [V,] ip,q,r,s) = JdpsiT,-T2) iAr\ ^^^k^X?-m- + fc, - fc,)^ - 
[II] [V,] (p, g,r,s) = ldps{r,-T2) tr^f.V, + _ _ (f - fc, A)^ - 

The full polarised amplitude of gg — > is given by summing all these expressions. 
6.2.3 Integration over internal kinematics 

To extract explicit forms for the conditions constraining the integration phase space, we shall 
decompose the vector i along three independent light-cone vectors ki, k2, k^ and 

£ = xki + yk2 + zk^^ + Ae^- (6.63) 

In terms of these new variables (x, |/, 2, A), we can rewrite the element of integration d'^i as 

1 X 1 

d'^i = -Vstu dx dydzdX = -\/ stu dxdydz —= [dX^ (■)Iva=a2 + ^^"^ (■)Iva=-a2] • (6-64) 
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From Eq. (6.63) and Eq. (6.21)-Eq. (6.23), we find 



{xy)s — {xz)u — {yz)t + A e^., 
^{ys - zu), 



{xs — zt), 

{x{s + u) + y{s + t)- z{u + t)), 



■-{xu + yt), 



and 



S{2e.{k, + k2)) 



{u + t) 



6 -z- 



u + t 



(x + y) 



(6.65) 
(6.66) 

(6.67) 

(6.68) 

(6.69) 

(6.70) 



The integration over z can thus be done. The second 5-function in Eq. (6.42) allows to 
integrate over A^: 



6{{i + ki'^ 
We have for d'^t. 



1 V stu 
Au + t 
1 \/ stu 
lu + t 

1 V stu 



SZL /_2 — \ / — 9 — \ 9 

[x + X) -\ [y^ — y) + m 



u + t 



{dxdy)dz5 ^z -j— ^ 

(dxdy) 



u + t 

(x + y) 
1 



(6.71) 



d{X' 



-^A^ + [••])(■) 



^ + X + ^(f -y) + ^m2l 

u+t I U^"^ ^ ' su J 



Au + t 



(dxdy) 



(6.72) 



We can rewrite the latter square root {i.e. A), 



-su 



u + t 



x'^ + X ^ — {y"^ — y) + 



u + t 



u 



su 



-su 



u + t 



{x - xi{y)){x - X2iy)) (6.73) 



so that the condition X2{y) < x < Xi{y) giving the limits of integration over x appears 
clearly: 



xi,2iy) = --±7;\h--^-[y^-y + 



2 2 



u 



u + t 



St 



(6.74) 



From this expression we also have a constraint upon the y-integration. Similarly, the latter 
square root can be rewritten: 



,t / „ _ u + t , 

1 — 4— { y — y -\ 

u \ st 



-4:-{y -yi){y -y2)- 



(6.75) 
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We see that iit/u > 0, then the expression is real for y2 < y < yi, where 



1 1 lu + t s-4m^ 1 1 J(s - M2)(s -4m2) 

= 2 ± 2 V ^ ^ = 2 ± 2 7^^! ^'-''^ 

From the boundary values yi and ?/2, we also infer that s > Af^ and s > Am'^. The 
integral from Eq. (6.72) then becomes: 

dH{.) = / dy / dx^ 



Au + t 



m Jx2{y) y/-^(a: - Xi{y)){x - X2{y)) 



-^A^ - xi{y)){x - X2{y)) 



Au + t 



6.2.4 Organising the integration 

In the same fashion, we rewrite the propagators in terms of the new variables of integration 
{x,y, z, A): 

(6.78) 

EI' (6-79) 

t 

(6.80) 





1 


u + 1 1 


u + t 1 




— m^) 


sn X — -y 


st y -jx' 


1 




U + t 1 


u + t 1 




— m^) 


X — -y -\ 


- St w - 


1 




1 1 


1 1 


((£ - A;2 + A;4)2 


— m^) 


u x+{{y-l\ 


) t y - 1 + f x ' 


1 




1 1 


1 1 


({£ + ki-hy 


— m^) 


U X + + 1 


t y + f(x+l)' 



(6.81) 
(6.82) 

Concerning the vertex factors Fj, as discussed previously, we shall employ either the 
dipolar form or the gaussian form both with the shifted argument pj = pf— ''^''^p instead 
of pf: The relative momentum squared pj is calculated from the momenta at the meson 
vertex 

. -ux-tiy-l) , s-4m2 + (2x + l)u+(2v-l)t 
2 ^ 1^ L ^ V i \^ J_ 

2 4 4 .gggX 

2_u{x+l)+ty 2 _ s - 4m2 - (2x + l)n - (2?/ - l)t ^' ^ 

P2 = ^ + - — - I 

and we therefore obtain the following expressions for the argument of the vertex functions 
at the considered relative momenta: 

2 f s-4m^ + {2x + l)u+{2y-l)t , (ux + t{y - 
Phsh = ( + Jjr2 ) (6.84) 

P2,sh = 



4 4M2 
s - 4m2 - {2x + l)u - {2y - l)t {u{x + I) + ty - M'^f 
4 ^ 4M2 
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and 



N 



(1 + 



A2 



or Fj = A^exp(- 



Pi,sh - 



All integrals that we have to calculate can be recast in the following forms: 
1 1 „ , , 1 \fstu 



d 



PS 



{e^ - m"^) {{i - k2 + hy - m^) 



r, X [. 



iTT^ 4(m + t) 



u + t 



[dxc 



1 



(6.85) 



(6.86) 



1 



ips 



{P - m"^) {{£ + ki ~ - m'^) 
~^ [dxdy] 



r,; X\... 



' stu 



X^x-^y x + ^y + 1 



87r2 4(m + 1) 
r,; X [. . .1 



(6.87) 



d 



PS 



{{i + ki- k2Y - m?) ((£ -k2 + k^f - m2) 



r,, X [. 



u + t 



[dxc 



VX'^-iv-'^x+iiy-l) 



'stu 



87r2 4(m + 1) 



r. X [. 



(6.88) 



1 



1 



ips 



((£ + A;i - k^Y - m?) ((£ + k^ - k^f - m^) 

1 1 



~'~ [dxdy] 



su^ 



Jy^x -^y-^x + -y + l 



r. X [. . .] = 
r,, X \. . 



1 \/ stu 

hx^4:{u + t) 



(6.89) 



with A = ±y'-^(x - xi{y)){x - X2{y)) in [. . .]. 

Let us introduce a short-hand notation for the denominators, the vertex functions and 
the jacobian -j=^: 



Ds{di,ej,Tk) = 



1 1 



1 



X"^ X + di X + Co- 



wit h 



di = — y., 
u 



do 



t u + t 

-y + 



u u 

ei = +-{y-l), 
u 

t _ 

62 = H — y + l- 

U 



(6.90) 

(6.91) 
(6.92) 
(6.93) 
(6.94) 
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defined so that 



1 



u + t 1 



(£2 - m2 



su X + di 
u + t 1 







— /c2)2 — m2) 






1 












1 



1 1 



((^+fci-fc4) 



u X + ei 
1 1 
M a; + 62 ' 



The amphtudes can then be written as^: 



[I 
[I 
[11 

[11 

[I] 
[I 
[11 
[II 



[III] {p^q.r, s) 

[IV] (p,g,r, s) 
[III] iP,q,r,s) 



-C J [dxdy]D3{di,ei,ri)tTf_!jjj 
+C j [dxdy]Dsidi,e2,T2M 
+C j [dxdy]D3{d2-, ei,T i)iT 



[IV](p,g,r, s) = -C J[dxdy]D3{d2,e2,T2)tT^j 

[Vi](p,g,r, s) = +1 

[V2](p,g,r, s) = -I 

[Vi](p,g,r, s) = -I 

[V2](p,g,r,s) = +1 



with 



C 



' stu 



327r2 sm2 



(6.95) 



pqrs 
'l-III 


(6.96) 


jpqrs 
I-IV 


(6.97) 


pqrs 
II-III 


(6.98) 


pqrs 
'l-IV 


(6.99) 


-D,{d,,e2j2)MZv, 


(6.100) 


- 1^3(^/1, ei, 12)] trr4 


(6.101) 


-Z}3(rf2,e2,r2)]tr?f_V, 


(6.102) 


- 1)3(^2, ei,r2)]trff1V2 


(6.103) 



(6.104) 



Apart from the denominators and the vertex function in D3, the only dependence on 
left is contained in the trx-v ^ through scalar products involving Indeed, as we have seen, 
the decomposition of £ will introduce x, z and A factors. The z factors are dropped in 
favour of x and y. On the other hand, as A has non-linear dependence on x and y, we shall 
keep it explicitly. 

It is also easy to convince oneself that we can decompose any trx-v as: 



Lm,n 



(6.105) 



the + or — signs before C simply come from Eq. (6.95). 
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where l,m,n = (0, 1, 2, 3, 4) and I + m + n < 4. This follows simply from the number of i in 
each trx-Y- Nothing here will be stated about Kf^^'^ ^ except that they do not contain x 
and i.e. the integration variables. 

Therefore, we are simply to perform one type of integrals 

/(rfi,e,-,rfc;/,m,n) = j [dxdy]\'x^rD^{due,,Tk) (6.106) 

[dxdy\— — _ — Tfc 



^J-^A^-Xl{y)){x-X2{y)) 



X + di X + Cj 



[dxdy]^^-—- —— r, . (6.107) 

^/X2 X + di X + Cj 

We already know that the integrals will vanish for / odd, as the two branches of the integrals 
over cancel each other. 

From this feature, we can also draw some conclusions on some helicity amplitudes. In- 
deed, the only way to make A appear and not A^ is through the scalar product i.eT = A; all 
other possible scalar products with i (see Eq. (6.65)) do not involve an odd power of A. Since 
eT appears solely in ej^, we can check that all amplitudes with an odd number of particles 
polarised along er (vector orthogonal to the collision plane) vanish, as should be. 



6.3 Polarised cross sections 

Combining Eq. (6.105) with the expressions for the polarised amplitudes as in Eq. (6.96) 
and the generic expression for the integrals on the internal phase space Eq. (6.106), we have 

M{p,q,r,s) = ^[X][Y](p,g,r,s) = ^ ^ A'J^;^/(cii, e^-, T^; /, m, n), (6.108) 

X,Y X,Y l,m,n 

the coefficients Ki are determined by FORM [97] and are functions^ of s, t and u whereas 
the integrals I{di,ej,rk',l,m,n) are performed numerically for each sets of s, t and u after 
having chosen the parameter defining F. 

Squaring the sum of the expressions [X][Y](p, g, r, s), multiplying them by ^ for the 
averaging on initial colour, by ^ for the averaging on initial spin, we get |A^(p, q, r, s)p to be 
plugged into Eq. (6.20). This in turn gives the polarised differential cross section ^^fiEi^Ijfl^ 

As we are concerned with the polarisation effects only for the vector meson, we shall 
sum the various helicity contributions for the gluons and define the following polarised cross 



^ As now, we are dealing back with hadronic and partonic processes, we put the circumflexes back for 
the partonic Mandelstam variables. 
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sections: 



dar^ _ ^ da{p,q,Ti,s) 



dt ^ rr dt 

p,q,s=Ti,T2 



1^ TS (6.109) 



dt ^ dt 

p,q,s=Ti,T2 

doL _ da{p,q,L,s) 
di ~ ^ ^ di 

p,q,s=Ti,T2 

The hadronic cross sections are therefore obtained with 



dcTh , 2sPTGi{xi)G2ix2ixi)) dah ,n -,-,r.\ 

— dXi r -pr, (6.110) 



dydPr jx^^" \/s{\/sxi — Exey) di 
for h = Ti,T2 or L. Specifying s, y and Pt, we get the final resuh. 

6.4 Numerical results and comparison with the CSM 

The most natural comparison of our calculation is with the LO CSM predictions. Indeed, 
the structure of that calculation is similar to that of our model with solely the imaginary 
part of the amplitude. We recall here that the two quarks forming the meson are put on 
their mass shell in the CSM. 

The comparison will be done with the unpolarised partonic cross section for gg — ^ ^Sig 
obtained in [24]: 



2\2l 



da _ 207r^Maf|V;(0)p [s\s - M^] + [lP{i - M^] + [u\u - M'^ 
'dt ~ 9P {s - M^fii - M^y{u - M2)2 ■ 

We shall use thegluonpdf's from CTEQ6L0 fit [98] and MRST2001LO [99] evaluated 
at the scale Q = + P^. The coupling constant is computed as 



with 5 active flavours and I^qcd = 0.165 GeV [98]. Indeed, Q is above as soon as Pt 
reaches 3-3.5 GeV for the J/ip . 



6.4.1 J/i/j at the Tevatron 

Setting y/s to 1800 GeV and considering the cross section in the range |?7| < 0.6, we get 
the following results for J/ip production at the Tevatron in the experimental conditions of 
the CDF collaboration whose measurements are detailed in chapter 2. The flrst plot (see 
Fig. 6.4) shows our prediction {axox, c"t and a^) for nic = 1.87 GeV and A = 1.8 GeV to be 
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compared with the LO CSM (which does not depend on rric) and the CSM fragmentation. 
We see that our prediction for the imaginary part of the amphtude is lower than the LO 
CSM, except at large Pt and is mostly longitudinal. 

On the other hand, these results include the gauge-invariance restoring vertex whose 
choice is arbitrary. The conclusion drawn here might be only relevant for the GIPA, i.e. the 
standard pQCD set of diagrams with a phenomenological vertex for qq ^ Q and the 
gauge-invariance restoring term of Eq. (5.59). Differences between possible choices of gauge- 
invariance restoring vertices are described by what we call autonomous contributions, their 
effects are analysed in the next section. 




Figure 6.4: Polarised (cr^ and cr^) and total {aTor) cross sections obtained with a gaussian vertex 
functions, nic = 1.87 GeV, A = 1.8 GeV and the MRST gluon distribution, to be 
compared with LO CSM and fragmentation CSM contributions. 



Then, we illustrate (see Fig. 6.5 (a)) the effect due to a change in the form of the vertex 
functions, i.e. from a dipole to a gaussian. The effects due to a variation in A (see Fig. 6.5 (b)), 
rUc and A (see Fig. 6.5 (c)), due to a change in the scale Q = -|- (see Fig. 6.5 (d)) 

are also shown. 
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Figure 6.5: (a) Comparison between the cross section obtained with the dipole and the gaussian 
vertex functions (MRST); (b) Variation of the cross section due to a change in A for 
a fixed value of the quark mass (MRST); (c) Variation of the cross section due to a 
change in nic and A (CTEQ); (d) Variation of the cross section due to a change of 



a factor two in the scale Q 



P| (CTEQ). 



As the changes in the vertex function illustrated above are of little influence on the final 
results, all the following plots for the J/ip, unless specified differently, will be made for 
rric = 1.87 GeV, A = 1.8 GeV, Q = + P| and a gaussian vertex function. 
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6.4.2 T(IS') at the Tevatron 

Here are similar plots for the T(IS') at the Tevatron (77 = 0). The deviation of the GIPA 
with the GSM is slightly more marked, perhaps indicating that the deviation behaves like 
a power of In the following, all plots for the T(IS'), unless specified differently, will be 

made for m^, = 5.28 GeV, A = 6.0 GeV, Q = a/ + and a gaussian vertex function. 

1 i 1 1 ' 1 1 ' 1 i 
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Figure 6.6: Polarised (cr^ and cr^) and total {(Jtot) cross sections obtained with a gaussian vertex 
functions, nib = 5.28 GeV, A = 6.0 GeV and the CTEQ gluon distribution, to be 
compared with LO CSM and fragmentation GSM contributions 
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Figure 6.7: Variation of the cross section due to a change in A for a fixed value of the quark 
mass (GTEQ) 
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6.4.3 J/^ at RHIC 

Here are plots for the J/^ at RHIC (v^ = 200 GeV) (for -0.5 <y< 0.5). 
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Figure 6.8: Polarised [ax and cr^) and total (cttot) cross sections obtained at ^/s = 200 GeV 
with a gaussian vertex functions, rric = 1.87 GeV, A = 1.8 GeV and the CTEQ gluon 
distribution, to be compared with LO CSM. 
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Figure 6.9: Polarised [ar and ai) and total [aTor) cross sections obtained at ^/s = 200 GeV 
with a gaussian vertex functions, rric = 1.87 GeV, A = 1.8 GeV and the CTEQ gluon 
distribution, to be compared with LO CSM. 



132 



J/ip and T production at RHIC and at the Tevatron 



6.4.4 ip' at the Tevatron 

Here we plot the cross section for the ip' at \/s = 1800 GeV with rric = 1.87 GeV, A = 1.8 GeV 
and for three values of the node parameter anode- 0.8, 1.46, 2.0 GeV (see Eq. (4.17) and 
Tab. (4.3)) We see that this parameter can have a significant infiuence. This comes from the 
fact that the normalisation strongly depends on the node position, whereas the production 
process is not very sensitive to this parameter. 

In the following, we shall use anode = 1-46 GeV, nic = 1.87 GeV, A = 1.8 GeV and a 
gaussian vertex function. 
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Figure 6.10: Polarised (ctt and ai) cross sections obtained at ^/s = 1800 GeV with a gaussian 
vertex functions, uic = 1.87 GeV, A = 1.8 GeV and the CTEQ gluon distribution 
for various values of anode- 



6.4.5 Taking stock of the situation 

Following the discussion concerning the J/ip at the Tevatron and RHIC, the T at the Teva- 
tron and the ip' at the Tevatron, we may draw some conclusions concerning the effects of 
changing the energy involved in the collision, the mass of the quarks and the internal struc- 
ture of the quarkonium. We can first say that the main features of our predictions are not 
drastically affected by these changes. 

To what concerns the deviation between the GIPA and the LO GSM, it seems that it is 
a permanent feature at low Pt- As we have said before, it indicates that the real part can 
be large. One important result of our study of non-static effects - or relativistic effects - is 
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the existence of a large uncertainty in ip' production linked with the position of the node in 
the vertex function. 

We have shown that it is possible to go beyond the static approximation and that our 
procedure does not show any instabilities and uncontrollable behaviour due to changes in 
and A. Our procedure to normalise the vertex function also seems to work and to provide 
sensible values. 

We have the hope that tuning the internal dynamics of the quarkonium in production 
processes could be applied to their study in media; for instance, to study their suppression 
due to an evolution in their size or due to an alteration of the binding of their constituting 
quarks. 

However, numerous choices of gauge-invariance restoring vertices are possible. No firm 
conclusions should be drawn before we have an idea of the effects due to a change in them. 
The easiest way to get insights on this issue comes from the analysis of contributions that 
were derived in section 5.3. This is the subject of the following study. 



6.5 Analysis of autonomous vertices 

We now consider the contributions from autonomous vertices as in Eq. (5.65). We have 
simply to introduce sub-amplitude [VI], [VII], . . . , whose expressions are obtained from the 
vertex chosen. 

We limit ourselves to choices for the vertex coming from Eq. (5.69) with constant factors, 
we have therefore three possibilities, that we rewrite in the conventions of this chapter: 

1. (El - E2)a7^g^ (Ei - E2)a7'"A;f the "a-term" 

2. (El - Es)/? {p + p'Y{p + p'Y ^ (El - T2)(3 {2£ + k, - k2y{2£ + ki - k^f the "/3-term" 

3. (El - Es)^ g^"" (El - Ea)^ gP" the "^term" 

Let us consider the first one. If we include it in our calculation, it will induce the following 
sub-amplitude to be added to the set [III] + [IV] + [V], 

[VI] = -la (El - E2) {-^gs)n^'^{^ - k2)iYklu{l + k^)e'-; {k^)el\k^) . (6.112) 

Let us now have a closer look at the product k1e'^p{k^). The latter is zero for r = T 
(Eq. (6.32)), whereas for r = L it is (Eq. (6.36)) 

^'iks).h = '-^. (6.113) 



This contribution will thus produce solely longitudinal vector mesons. It is the ingredient 
needed for a possible solution of the polarisation issue. We shall come back to this fact when 
we analyse the numerical results. 
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In practice, the introduction of these contributions does not modify much the procedure 
to get the polarised cross section as explained in the previous sections. The main difference 
is the introduction of a new quantity 



which appears in the following new integrals, 



[I][VI](p,g,r,s) = +C' j[dxdy][D2{duTr)-D2{drj2)Wi%i (6.115) 
[II][VI](p,g,r,s) = -C j[dxdy][D2{d2ji)-D2{d2j2)]tAT-vi^ (6-116) 



where C is -^^^^j^ and the quantities ti^'^^yj and ti^jY^yj are 



trf'-^^ = - Tr [T'^T'T^] Tr [7^^ (/ + m) 7^((/ - ^2) + m) 



tr?f_V, = - aglTi [t'T^^] Tr ((/ + ^1 - h) + 7^((/ - h) + 
X7'^((/ + h) + ^)] fc4'^^(fci)4(fc2)er(^3)C(^4) 



(6.117) 



(6.118) 



6.5.1 Numerical results 

To get a first idea, we have introduced these vertices in the amplitude one by one with a, (3 
or ^ equal^ to 1. In Fig. 6.11, we show the individual contributions of these terms. Note that 
some of these produce only one helicity. One example already foreseen concerns the a-term 
whose transverse cross section is zero. However, if we want to fit the data, the autonomous 
vertices will need to be dominant, hence the GIPA and any interferences terms will not be 
important. 



Note that in this first appHcation the coefScicnts a, /3 and ^ do not have the same dimension. 
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Figure 6.11: Polarised (cry and ai) cross sections for the J/'ip at the Tevatron obtained with 



three different autonomous vertices alone. Note that a 



0"} 



0. 



The most striking result is that this latter contribution (a-term) has a slope compatible 
with the data. The contribution from the ,^-term has no longitudinal contribution^ and its 
transverse component is slightly bigger than the GIPA. Yet, one might expect, since its 
dimension is lower, that the actual value of ^ would enhance it. The contribution of the 
/3-term has no remarkable characteristic and considering its dimension, we do not expect /3 
to enhance its contribution. 



6.5.2 Fitting the data 

We are then going to make use of these two unconstrained structures {a and ,^-terms). One 
is purely transverse whereas the other is purely longitudinal. This is a motivation to try 
to reproduce the behaviour of the polarisation parameter Opoi as measured by CDF (see 
section 2.5). 

First we shall rescale the multiplicative factors a and ^ to get dimensionless parameters: 



a 



a 



f"^' (6.119) 

Note that this rescaling is purely arbitrary. It is far from obvious which quantity among mg 
(mass of the quarkonium), mq , A, !s.qcdi V^, ... is the most suitable. However, except for 
^QCDi they are all of the same order. Besides, the uncertainties coming from the missing 
real part may modify the final values obtained in the fit. 

^ This can be also foreseen since g^" e'^* {kz)e^* {k^) = e^* {k^) .e'^* {ki) is zero ior r — L (see section 6.2.1). 
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J 14) at the Tevatron 

Setting a to 8 and ^ to 37.5, we reproduce the cross section for J/ip at -y/i = 1800 GeV [43] 
as displayed in Fig. 6.12. However since ai is much larger than ctt for Pt > 10 GeV, we 
predict that Upoi will be negative for -Pr ^ 6 GeV at variance with COM fragmentation 
predictions, for which apoi is significantly positive. This fact can easily traced back to the 
steeper slope of the ,^-term compared to the a-term (see Fig. 6.11 (b)). 



T(15') at the Tevatron 

Setting a to 8 and ^ to 10, we reproduce the cross section for T(IS') at y/s = 1800 GeV [60] 
as shown in Fig. 6.13. We observe the same behaviour as in the J/ip case to what concerns 
(Tl relative to ar and we predict that apoi will be negative for -Pr ^ 8 GeV. 



J/ij at RHIC 

Since the Pt range of the data [55] is not very wide and the error bars are rather large, we 
use the same values for a and ^ as for the Tevatron energy. Furthermore, our calculation 
(Fig. 6.14) agrees with the data, even though slightly improved agreement could be obtained 
with a different set of values for a and ^. To what concerns polarisation, we predict that 
apoi will be positive for Pt ^ 4 GeV and negative for larger values of Pt- 

Note that our LO GSM curve includes solely gluon fusion at variance with curves of 
Nayak et al. [67] which include all parton-fusion processes. 



ip' at the Tevatron 

Setting anode to 1.46 GeV and a to 27.5, we get the cross section for ip' at ^/s = 1800 GeV 
as shown in Fig. 6.15. This fits perfectly the data obtained by GDF [43]. There is no 
need to introduce any contribution from the (^-term and we predict that the cross section is 
longitudinally polarised. 
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Figure 6.12: Polarised {ar and (Tl) and total {aTor) cross sections obtained for J/ip at = 
1800 GeV with a = 8 and ^ = 37.5 to be compared with LO CSM, the fragmentation 
CSM and to the data of CDF [43]. 
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Figure 6.13: 



Polarised {ax and (Tl) and total {aTor) cross sections obtained for T(15') at ^/s = 
1800 GeV with a = 8 and 10 to be compared with LO CSM and to the data 
of CDF [60]. 
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Figure 6.14: Polarised (cr^ and ai) and total {(Jtot) cross sections obtained for J/i}) at = 
200 GeV with a = 8 and ^ = 37.5 to be compared with LO CSM and to the data 
of PHENIX [55]. 
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Figure 6.15: Polarised (ctt and ai) and total (error) cross sections obtained for tp' at ^/s = 
1800 GeV with a = 27.5 to be compared with LO CSM, the fragmentation CSM 
and to the data of CDF [43]. 
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6.5.3 Refining for the polarisation 

Following the discussion of Tab. (1-14), we know that fragment at ion- like processes should 
dominate for sufficiently large Pt- We have not yet included these in our calculation. Since 
the COM contribution is mainly transverse^" and can account with appropriate LDME 
(Long-Distance Matrix Element) value for the CDF data [50], we are going to combine 
it with our contribution from the a-term which is purely longitudinal and which can also 
account for the data. 

The amplitude to produce a ^Si quarkonium Q by fragmentation of a colour octet state 
^5"^^^ is given by 

\M{gg - 'Sf'g - 'S,g)\' = Y,\A{gg - ' Sf^ g)\' {O'^f sf\) (6.120) 

T 

with for a transverse polarisation state [50] 



Y.\^{99^'Sfg) 



T 



X 



144M3 (s - M2)2 

27{si + iu + us) - 19M^ 
[{s- M'^){i- M^){u- NP)f 



M^f + + u^ + 2M^ 



X 



(6.121) 



and for a longitudinal polarisation state [50] 

VUfoo ^ 'S^'^a)\' - J±^^M!^(P , 27(gt + fe + ng)-19M^ 

(6.122) 

and where the LDME {0^['^S^^]) depends on Q. The longitudinal cross section is suppressed 
by and therefore always small compared to its transverse counterpart. 

Combining our contribution and COM fragmentation with smaller coefficients (roughly 
decreasing the cross section of each contribution of a factor of 2), we should recover a good 
description of the measured cross section as a function of Pt but also a reasonable description 
of the evolution of the polarisation parameter function of Pt ■ 



J/ip at the Tevatron 

Setting a to 4.2, ^ to 26.5 and {O^/'^'i^ Sf^)) to 4 x lO'^ GeV^ we got the cross section 
for the J/ip at = 1800 GeV shown in Fig. 6.16. This fits rather well the cross section 
measured by CDF [60]. To what concerns polarisation data [52], they are only available for 
prompt J/ip and not for the direct component of the signal for which our prediction are 
relevant. However, we can say that we obtain the same trend (see Fig. 6.17). 

We shall not consider its longitudinal part in plots, even though we included it in our calculation for 
the sake of completeness, plotting them would have diminished the readability of the plots. 
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Figure 6.16: Various cross sections (see legend) obtained for the J/^ at = 1800 GeV with a = 
4.2 and ^ = 26.5 and with a COM contribution with {0-^/^{^s[^^)) = ixlQ-^ GeV^ 
to be compared with LO CSM and to the data of CDF [43]. 



+ 



CM 



CDF DATA (PROMPT) 




a=4.2 ^=26.5 + C0Mld^;|e=4 10"' 



10 



Pj (GeV) 



12 14 16 18 



20 



Figure 6.17: The polarisation parameter apoi = "J^j^2Jl obtained for the J/if) at = 1800 GeV 



with a = 4.2 and ^ = 26.5 and with a COM contribution with (C-^/^(3^f ^)) 
4 X 10-3 GeV^ 
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at the Tevatron 

With anode = 1-46 GeV, setting a to 17.5 and {0^'{^S[^^)) to 4.5 x lO^^ GeV^ we get the 
cross section for the tp' at ^/s = 1800 GeV shown in Fig. 6.18 and the polarisation parameter 
apoi, in Fig. 6.19. The agreement is quite good both for cross section and for the polarisation. 
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Figure 6.18: Various cross sections (see legend) obtained for tp' at y/s = 1800 GeV with a = 17.5 
and with a COM contribution with {O"^' {^S[^'^)) = 4.5 x 10^^ Q^y^ compared 
with LO CSM and to the data of CDF [43]. 
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Figure 6.19: The polarisation parameter apoi = "J^j^2Ji obtained for the ip' at ^/s = 1800 GeV 

with a = 17.5 and with a COM contribution with {0^'{^s{^^)) = 4.5 x 10"^ GeV^ 
to be compared with the CDF measurements [52]. 
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In brief 

Our fitted values for a and ^ are: 

1. Witliout COM contribution: 

• for tlie J I'll) at tfie Tevatron: a=8 - .^=37.5 

• for the J/ip at RHIC: cf. Tevatron 

• for the T(IS') at the Tevatron: a=8 - .^=10 

• for the ip' at the Tevatron: a=27 - ^: < ^ < 30 (not constrained) 

2. With COM contribution: 

• for the J/%Ij at the Tevatron: a=4.2 - ^=26.5 - {0'^''^'(^Sf'^)) = 4.0 x lO'^ GeV^ 

• for the ip' at the Tevatron: a=17.5 - ^: < ^ < 30 (not constrained) - 

(C^'(35f^)) =4.5 X 10^3 Q^y3 



Chapter 7 

Conclusions and outlook 



"Use NRQCD !". I shall never forget this wise and humoristic reply of Eric Braaten to 
me when I asked him what to do to go beyond the CSM - and its static approximation - 
without having problems with gauge-invariance. Indeed, that approximation implying that 
whatever the hard part could be, the non-relativistic nature of the quarkonium should reduce 
the problem to a static one, was seriously disturbing me. Of course, NRQCD was a way to 
probe the internal dynamics of the quarkonium, but still, the quark in the bound state were 
always considered at rest. 

Now, are we in position to reply that we could go beyond this static approximation 
and that we included relativistic effects in the calculation ? On the one hand, yes. We 
have introduced a relative momentum in the problem, we have tackled the issue of gauge 
invariance, we have constrained our model and we have applied it to processes for which the 
CSM was failing. 

On the other hand, no. We have not calculated the real part of the amplitude and the 
contribution from the imaginary part is lower than the LO CSM at low Pt- However, we 
have shown that the choice of the gauge-invariance restoring vertex is not unique: this is the 
price to pay to introduce non-static effects. And, then, things are getting very interesting 
since we can reproduce the data, both from the Tevatron and from RHIC. This is the first 
time that a model can reproduce both the rate and the polarisation of the data. 

What we did 

Let us recall one by one the achievements of this thesis. As we have just said, our aim was 
to introduce non-static relativistic effects in the problem of hadroproduction of quarkonium, 
especially for the J/ip, ip' and the T for which data are numerous, reliable but in evident 
contradiction with the models derived from pQCD. 

In chapter 3, we studied the theoretical uncertainties affecting the LO CSM, as an indi- 
cation of those affecting all models. These uncertainties being a factor 3, our conclusion was 
that something new was to be introduced, like the colour-octet contributions brought in by 
NRQCD. 

The way that we have chosen to introduce non-static effects was via a phenomenologi- 
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cal 3-particle vertex, made of a Dirac structure inferred from the quantum numbers of the 
considered quarkonium and a vertex function inspired by phenomenological arguments. Of 
course, the vertex function was to be constrained. Firstly, its size parameter - teUing us 
how rapidly it falls when the relative momentum of the quarks inside the quarkonium in- 
creases - was chosen to be compatible with existing theoretical studies. Secondly, the factor 
multiplying the vertex - the normalisation - was fixed to reproduce the leptonic decay 
width. 

In the case of the ifj', we showed that the normalisation, or equivalently the decay width, 
was strongly affected by a change in the position of the node in the vertex function. 

Whereas the choice of the vertex function was driven by phenomenological considerations, 
the method, which we have proposed to restore gauge invariance was based on a rigorous 
and ab initio derivation exposed in chapter 5. First, we stated how gauge invariance was 
broken by the introduction of a 3-particle vertex function in the case of a pseudoscalar bound 
state. This was followed by the description of how to restore gauge invariance, keeping in 
mind that this procedure, which introduces new 4-particle vertices in the calculation, was 
not unique. 

We considered the same problem for vector particles associated with a gluon. The same 
conclusions were drawn, and we also studied in details how to parametrise the freedom left 
in the 4-particle vertex. Indeed, we have given simple and easily manageable Feynman rules 
for what we have called "autonomous" contributions, which correspond to different choices 
in the gauge-invariance restoring vertices. Note that the new contributions disappear in the 
case of point-like particles, for which there is no problem with gauge invariance. 

In chapter 6, we have applied our model to the hadropro duct ion of quarkonia associated 
with a gluon, with the sole consideration of the imaginary part of the amplitude. This 
constraint comes from the necessity of having internal quark on-shell to derive the new ver- 
tices. We made the calculations, first, for one choice of a gauge-invariance restoring vertex, 
then, for the same contribution with autonomous vertices whose weights were fitted to repro- 
duce the cross-sections and, finally, with the further addition of colour-octet fragmentation 
contributions to reproduce the polarisation data. 

What we have obtained 

The main result that we have obtained is that one autonomous contribution is naturally 
dominant and produces only longitudinally polarised quarkonia. Moreover, the Pt slope of 
this contribution is compatible with the data. Indeed, it is expected that any non vanishing 
scalar product with a helicity vector of a longitudinally polarised quarkonium should scale 
like s, which, at threshold, is proportional to P^. This explains why longitudinal components 
have a less steep slope. 

If we combine this contribution with colour-octet fragmentation, which is also naturally 
large and produces only transversally polarised quarkonia, we can reach a perfect agreement 
with all existing data and with a minimal number of free parameters. 

This combination with colour-octet fragmentation was actually motivated by the fact 
that we were missing fragmentation contributions, which are expected to be significant, due 
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to their Pt slope (see Fig. 1.14). 

Our new contribution is the only one, which is able to produce only longitudinally po- 
larised quarkonium. 

Of course, we obtained other results. We have managed to introduce relativistic effects 
in the sense that we took account of contributions with nonzero relative momentum - non- 
static configurations - between the two constituent quarks. Our estimate of the imaginary 
part is lower than the LO CSM at low P^, but of the same order at large Pt- We also got 
some interesting insights on the influence of the node of the 25" states. 

More generally we have provided a formalism to study effects that come from the internal 
structure of the quarkonia and that were never considered before. One very interesting appli- 
cation would be to study the effect of a variation of the binding potential on the production 
cross-section of quarkonia. This variation of potential could for instance be introduced to 
study the influence of the quark-gluon plasma. 

What we still miss 

Even though these results are new, there are still some drawbacks. The first is the missing 
real part. If we take the LO CSM as a benchmark, the real part seems to be dominant over 
the imaginary part at low Pt- One way to calculate it would be to use dispersion relations. 
However, due to the expected asymptotic logs behaviour of the partonic amplitude, we have 
to introduce a subtraction constant in the dispersion relation. Furthermore, this constant is 
a function of Pt- We could fix it by fitting the data at a given energy and make predictions 
for different energies. In the present situation, the data from RHIC have too small a Pt 
range to proceed through the analysis. The second missing point it that our description of 
the quarkonia does not take into account cuts in the vertices. 

These missing contributions make us believe that the fits which we have made are reliable 
only up to a factor of order 3. 

Another thing that we miss is an understanding of how the autonomous contributions 
arise. One of the three autonomous contributions that we have considered, namely the 
^-term {g^'^)-, could come from off-shell contributions. This can be seen easily by setting 
a different quark mass in the off- shell quark propagator in Eq. (5.33) and performing the 
classification of the different Dirac structures with their factors. However, for the other ones, 
we have currently no simple argument to motivate their presence, except that they exist in 
the standard contribution and that their slope in Pt is expected to be quite low, explaining 
why they dominate at large Pt- 

What is expedient to do 

A first necessary test of our formalism would be to apply it to production regimes other than 
hadroproduction. The application to photoproduction at HERA is the first that comes to 
mind. Moreover, results on photoproduction are compatible with NLO calculations within 
the CSM [53] and no contribution from colour-octet mechanisms are needed. In that context, 
a smaller value of the LDME's fixed by the fits on hadroproduction data would be welcome. 



146 



Conclusions and outlook 



while our model provides with smaller ones. For all these reasons, we do plan to consider 
photoproduction in a close future. 

Besides photoproduction, electroproduction is also expected to provide us with interest- 
ing tests and constraints. Compared to photoproduction, the kinematics would be slightly 
modified. There exist data from LEP and a comparison with those would also be useful. 

Finally, in the case of S-factories, the situation is even more appealing, since both 
the CSM and the COM diagree with Belle data. For some observables, the discrepancy 
amounts to more than one order of magnitude. Let us mention the double-charm produc- 
tion, e.g. e^e~ — >■ J/ijj + T]c, and the inclusive production rate. However, a confirmation from 
BaBar of theses discrepancies would be greatly welcome. 

In a totally different context, another attractive application is the modelisation of Gen- 
eralised Parton Distributions (GPD). In previous works, we have worked out these for the 
pion in a simple but rough model where the restoration of gauge invariance was achieved by 
symmetrising the vertex function 2) — > Pi + P2. The result were encouraging, however 
we have now the tools to greatly refine this model. An important asset of our approach is 
the possibility to cast off the obligation to work on the light cone. 

Moreover, as mentioned a couple of times, the application to production in media is 
opportune. For some time now, J/ip has been thought to be an ideal probe for quark-gluon 
plasma through its suppression. Besides the fact that the absolute rates are not completely 
understood - as we have seen throughout this work -, the way to introduce modification 
in the quark binding induced by the hadronic medium are not completely settled. Our 
formalism may be used to study these processes. 

If we now turn back to the formalism rather than to its possible applications and tests, 
we are thinking to try to justify the presence of these autonomous 4-particle vertices from 
first principles, e.g. from standard perturbative Feynman rules applied to some typical sets 
of diagrams. We would then obtain an estimate of the free parameters a, /3 and ^, possibly 
in a similar fashion than the velocity scaling rules in NRQCD. 

Another direct improvement of our formalism would be of course the consideration of the 
real part of the amplitude. We have already addressed the possibility of making use of dis- 
persion relation with the drawback which is the presence of arbitrary subtraction constants. 
Another way may be to generalise our gauge-invariance restoring vertices to off-shell quarks. 
Then, we would be in position to consider the real part at any s or to determine it at a 
precise value of s, providing us with the value of the subtraction constants of the dispersion 
relations. 

By showing that it was possible to build a model reproducing all available data on high- 
energy hadroproduction of quarkonia, we do believe that we have provided a new subject 
of discussion and research and, hopefully, shed some light on what may be considered to be 
the path to a resolution of this issue. 



Appendix A 

Acronyms and abbreviations 



AGS Alternating Gradient Synchrotron 

BNL Brookhaven National Laboratory 

BSE Bethe-Salpeter Equation 

CDF Collider Detector at Fermilab 

CERN Conseil Europeen pour la Recherche Nucleaire 

CES Central Electromagnetic Strip chambers 

cm. center-of-momentum 

CMP Central Muon Upgrade 

CEM Colour Evaporation Model 

CLEO CLEO(patra)^ or Cornell's Largest Experimental Object (or Operation or Organization) 

COM Colour Octet Mechanism 

CSM Colour Singlet Model 

CTC Central Tracking System 

DESY Deutsches Elektronen-Synchrotron 

DSE Dyson-Schwinger Equation 

EMC ElectroMagnetic Calorimeter 

FNAL Fermi National Accelerator Laboratory (Fermilab) 

GIPA Gauge-Invariant Phenomenological Approach 

GPD Generalised Parton Distribution 

HERA Hadron Elektron Ring Anlage (Hadron Electron Ring Accelerator) 

ISR Intersecting Storage Rings 

KEK Kou Enerugii Butsurigaku Kenkyusho 

LDME Long-Distance Matrix Element 

LEP Large Electron - Positron (storage ring) 

LO Leading order 



Named sugcsted by Chris Day: short for Cleopatra 
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Acronyms and abbreviations 



MuID 


Muon IDentifier 


MuTr 


Muon Tracker 


NLO 


Next-to-Leading Order 


NRQCD 


Non Relativistic Quantum ChromoDynamics 


pdf 


parton distribution function 


pQCD 


perturbative Quantum ChromoDynamics 


QCD 


Quantum ChromoDynamics 


QED 


Quantum ElectroDynamics 


RHIC 


Relativistic Heavy Ion Collider 


RICH 


Ring Imaging C(h)erenkov 


SCI 


Soft Colour Interactions 


SLAC 


Stanford Linear Accelerator 


SPEAR 


Stanford Positron Electron Accelerating Ring 


svx 


Silicon Vertex Detector 



Appendix B 
Normalisation 

B.l Results relative to the normalisation 



A 


JV 


A 


JV 


1.00 


5.95 


1.63 


2.64 


l.OC 


5.36 


1.69 


2.48 


1.13 


4.86 


1.76 


2.34 


1.19 


4.44 


1.82 


2.21 


1.25 


4.07 


1.88 


2.09 


1.32 


3.76 


1.95 


1.98 


1.38 


3.48 


2.01 


1.88 


1.44 


3.23 


2.07 


1.79 


1.51 


3.01 


2.14 


1.70 


1.57 


2.82 


2.20 


1.62 



Table B.2: Normalisation for a dipolar form for J/Tp as a function of A for rric = 1.60 GeV. 



1.00 


7.93 


1.63 


3.17 


1.06 


7.04 


1.69 


2.96 


1.13 


6.31 


1.76 


2.77 


1.19 


5.69 


1.82 


2.60 


1.25 


5.16 


1.88 


2.44 


1.32 


4.71 


1.95 


2.30 


1.38 


4.31 


2.01 


2.17 


1.44 


3.97 


2.07 


2.06 


1.51 


3.67 


2.14 


1.95 


1.57 


3.40 


2.20 


1.85 



Table B.4: Normalisation for a dipolar form for J/-ip as a function of A for nic = 1.69 GeV. 
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Normalisation 



1.00 


9.39 


1.63 


3.54 


1.06 


8.28 


1.69 


3.29 


1.13 


7.3C 


1.76 


3.07 


1.19 


6.59 


1.82 


2.87 


1.25 


5.94 


1.88 


2.69 


1.32 


5.39 


1.95 


2.52 


1.38 


4.91 


2.01 


2.38 


1.44 


4.50 


2.07 


2.24 


1.51 


4.14 


2.14 


2.12 


1.57 


3.82 


2.20 


2.01 



Table B.6: Normalisation for a dipolar form for J /ip as a function of A for rric = 1.78 GeV. 



1.00 


10.61 


1.63 


3.84 


1.06 


9.31 


1.69 


3.56 


1.13 


8.23 


1.76 


3.31 


1.19 


7.34 


1.82 


3.09 


1.25 


6.59 


1.88 


2.88 


1.32 


5.95 


1.95 


2.70 


1.38 


5.40 


2.01 


2.54 


1.44 


4.93 


2.07 


2.39 


1.51 


4.52 


2.14 


2.25 


1.57 


4.16 


2.20 


2.13 



Table B.8: Normalisation for a dipolar form for J/ip as a function of A for nic = 1.87 GeV. 





N 


rrtc 


N 


1.60 


5.95 


1.74 


8.81 


1.61 


6.33 


1.76 


9.03 


1.63 


6.67 


1.77 


9.25 


1.64 


6.99 


1.78 


9.46 


1.66 


7.29 


1.80 


9.66 


1.67 


7.57 


1.81 


9.86 


1.69 


7.84 


1.83 


10.05 


1.70 


8.10 


1.84 


10.24 


1.71 


8.34 


1.86 


10.43 


1.73 


8.58 


1.87 


10.61 



Table B.IO: Normalisation for a dipolar form for J/tp as a function of the c quark mass for 
A = 1.0 GeV. 



1.60 


3.39 


1.74 


4.54 


1.61 


3.55 


1.76 


4.63 


1.63 


3.69 


1.77 


4.71 


1.64 


3.82 


1.78 


4.79 


1.66 


3.94 


1.80 


4.87 


1.67 


4.05 


1.81 


4.95 


1.69 


4.16 


1.83 


5.02 


1.70 


4.26 


1.84 


5.09 


1.71 


4.36 


1.86 


5.17 


1.73 


4.45 


1.87 


5.24 



Table 8.12: Normalisation for a dipolar form for J/ip as a function of the c quark mass for 
A = 1.4 GeV. 



B.l Results relative to the normalisation 
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1.60 


2.25 


1.74 


2.82 


1.61 


2.33 


1.76 


2.86 


1.63 


2.40 


1.77 


2.90 


1.64 


2.47 


1.78 


2.94 


1.66 


2.53 


1.80 


2.98 


1.67 


2.58 


1.81 


3.02 


1.69 


2.63 


1.83 


3.05 


1.70 


2.68 


1.84 


3.09 


1.71 


2.73 


1.86 


3.12 


1.73 


2.78 


1.87 


3.16 



Table B.14: Normalisation for a dipolar form for J/^ as a function of the c quark mass for 
A = 1.8 GeV. 



1.60 


1.62 


1.74 


1.94 


1.61 


1.67 


1.70 


1.97 


1.63 


1.71 


1.77 


1.99 


1.64 


1.75 


1.78 


2.01 


1.66 


1.78 


1.80 


2.03 


1.67 


1.81 


1.81 


2.05 


1.69 


1.84 


1.83 


2.07 


1.70 


1.87 


1.84 


2.09 


1.71 


1.90 


1.80 


2.11 


1.73 


1.92 


1.87 


2.13 



Table B.16: Normalisation for a dipolar form for J/ip as a function of the c quark mass for 
A = 2.2 GeV. 



A 


JV 


A 


JV 


1.00 


5.36 


1.63 


2.43 


1.06 


4.84 


1.69 


2.29 


1.13 


4.40 


1.76 


2.16 


1.19 


4.02 


1.82 


2.04 


1.25 


3.70 


1.88 


1.94 


1.32 


3.42 


1.95 


1.84 


1.38 


3.17 


2.01 


1.75 


1.44 


2.95 


2.07 


1.67 


1.51 


2.76 


2.14 


1.59 


1.57 


2.58 


2.20 


1.52 



Table B.18: Normalisation for a gaussian form for J /ip as a function of A for rric = 1.60 GeV. 



1.00 


7.38 


1.63 


2.95 


1.06 


6.54 


1.69 


2.76 


1.13 


5.86 


1.76 


2.58 


1.19 


5.28 


1.82 


2.43 


1.25 


4.79 


1.88 


2.29 


1.32 


4.37 


1.95 


2.16 


1.38 


4.01 


2.01 


2.04 


1.44 


3.69 


2.07 


1.93 


1.51 


3.41 


2.14 


1.83 


1.57 


3.17 


2.20 


1.74 



Table B.20: Normalisation for a gaussian form for J/ip as a function of A for nic = 1.69 GeV. 
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Normalisation 



1.00 


9.00 


1.63 


3.34 


1.06 


7.90 


1.69 


3.11 


1.13 


7.00 


1.76 


2.90 


1.19 


6.26 


1.82 


2.71 


1.25 


5.63 


1.88 


2.54 


1.32 


5.10 


1.95 


2.39 


1.38 


4.65 


2.01 


2.25 


1.44 


4.25 


2.07 


2.13 


1.51 


3.91 


2.14 


2.01 


1.57 


3.61 


2.20 


1.91 



Table B.22: Normalisation for a gaussian form for J/ip as a function of A for rric = 1.78 GeV. 



1.00 


10.44 


1.63 


3.68 


1.06 


9.10 


1.69 


3.41 


1.13 


8.02 


1.76 


3.17 


1.19 


7.12 


1.82 


2.95 


1.25 


6.37 


1.88 


2.76 


1.32 


5.74 


1.95 


2.59 


1.38 


5.20 


2.01 


2.43 


1.44 


4.74 


2.07 


2.29 


1.51 


4.34 


2.14 


2.16 


1.57 


3.99 


2.20 


2.04 



Table B.24: Normalisation for a gaussian form for J/ip as a function of A for = 1.87 GeV. 





N 


rrtc 


N 


1.60 


5.36 


1.74 


8.34 


1.61 


5.73 


1.76 


8.59 


1.63 


6.07 


1.77 


8.84 


1.64 


6.40 


1.78 


9.08 


1.66 


6.70 


1.80 


9.31 


1.67 


7.00 


1.81 


9.54 


1.69 


7.28 


1.83 


9.77 


1.70 


7.56 


1.84 


10.00 


1.71 


7.83 


1.86 


10.22 


1.73 


8.09 


1.87 


10.44 



Table B.26: Normalisation for a gaussian form for J/ip as a function of the c quark mass for 
A = 1.0 GeV. 



1.60 


3.09 


1.74 


4.26 


1.61 


3.25 


1.76 


4.36 


1.63 


3.38 


1.77 


4.45 


1.64 


3.51 


1.78 


4.54 


1.66 


3.63 


1.80 


4.63 


1.67 


3.75 


1.81 


4.71 


1.69 


3.86 


1.83 


4.79 


1.70 


3.97 


1.84 


4.88 


1.71 


4.07 


1.86 


4.96 


1.73 


4.17 


1.87 


5.04 



Table 8.28: Normalisation for a gaussian form for J/ip as a function of the c quark mass for 
A = 1.4 GeV. 



B.l Results relative to the normalisation 
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1.60 


2.08 


1.74 


2.66 


1.61 


2.16 


1.76 


2.70 


1.63 


2.23 


1.77 


2.74 


1.64 


2.29 


1.78 


2.79 


1.6C 


2.35 


1.80 


2.83 


1.67 


2.41 


1.81 


2.87 


1.69 


2.46 


1.83 


2.91 


1.70 


2.51 


1.84 


2.95 


1.71 


2.56 


1.86 


2.99 


1.73 


2.61 


1.87 


3.02 



Table B.30: Normalisation for a gaussian form for J /ip as a function of the c quark mass for 
A = 1.8 GeV. 



1.60 


1.52 


1.74 


1.84 


1.61 


1.50 


1.76 


1.87 


1.63 


1.60 


1.77 


1.89 


1.64 


1.64 


1.78 


1.91 


1.66 


1.67 


1.80 


1.94 


1.67 


1.70 


1.81 


1.96 


1.69 


1.73 


1.83 


1.98 


1.70 


1.76 


1.84 


2.00 


1.71 


1.79 


1.80 


2.02 


1.73 


1.82 


1.87 


2.04 



Table B.32: Normalisation for a gaussian form for J/ip as a function of the c quark mass for 
A = 2.2 GeV. 



A 


JV 


A 


JV 


3.00 


2.84 


4.71 


1.43 


3.17 


2.61 


4.88 


1.35 


3.34 


2.40 


5.05 


1.29 


3.51 


2.23 


5.22 


1.22 


3.68 


2.07 


5.39 


1.16 


3.86 


1.93 


5.57 


1.11 


4.03 


1.81 


5.74 


1.06 


4.20 


1.70 


5.91 


1.01 


4.37 


1.60 


6.08 


0.97 


4.54 


1.51 


6.25 


0.93 



Table B.34: Normalisation for a dipolar form for T(IS') as a function of A for m;, = 4.80 GeV. 



3.00 


4.06 


4.71 


1.80 


3.17 


3.66 


4.88 


1.69 


3.34 


3.33 


5.05 


1.59 


3.51 


3.04 


5.22 


1.50 


3.68 


2.79 


5.39 


1.42 


3.86 


2.57 


5.57 


1.35 


4.03 


2.38 


5.74 


1.28 


4.20 


2.21 


5.91 


1.21 


4.37 


2.06 


6.08 


1.16 


4.54 


1.92 


6.25 


1.10 



Table B.36: Normalisation for a dipolar form for T(IS') as a function of A for rrib = 5.00 GeV. 
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Normalisation 



3.00 


4.88 


4.71 


2.04 


3.17 


4.37 


4.88 


1.91 


3.34 


3.94 


5.05 


1.79 


3. .51 


3.57 


5.22 


1.68 


3.68 


3.26 


5.39 


1.58 


3.86 


2.98 


5.57 


1.49 


4.03 


2.75 


5.74 


1.41 


4.20 


2.54 


5.91 


1.34 


4.37 


2.35 


6.08 


1.27 


4.54 


2.19 


6.25 


1.21 



Table B.38: Normalisation for a dipolar form for T(IS') as a function of A for nib = 5.28 GeV. 



3.00 


5.56 


4.71 


2.23 


3.17 


4.95 


4.88 


2.08 


3.34 


4.44 


5.05 


1.94 


3.51 


4.01 


5.22 


1.82 


3.68 


3.64 


5.39 


1.71 


3.86 


3.32 


5.57 


1.61 


4.03 


3.04 


5.74 


1.52 


4.20 


2.80 


5.91 


1.44 


4.37 


2.59 


6.08 


1.36 


4.54 


2.40 


6.25 


1.29 



Table B.40: Normalisation for a dipolar form for T(IS') as a function of A for nih = 5.50 GeV. 





N 




N 


4.80 


2.84 


5.17 


4.56 


4.84 


3.10 


5.21 


4.68 


4.87 


3.31 


5.24 


4.80 


4.91 


3.51 


5.28 


4.92 


4.95 


3.69 


5.32 


5.03 


4.98 


3.85 


5.35 


5.14 


5.02 


4.01 


5.39 


5.25 


5.06 


4.16 


5.43 


5.36 


5.09 


4.30 


5.46 


5.46 


5.13 


4.43 


5.50 


5.56 



Table B.42: Normalisation for a dipolar form for T(IS') as a function of the b quark mass for 
A = 3.00 GeV. 



4.80 


1.77 


5.17 


2.53 


4.84 


1.89 


5.21 


2.59 


4.87 


1.99 


5.24 


2.64 


4.91 


2.08 


5.28 


2.69 


4.95 


2.16 


5.32 


2.74 


4.98 


2.23 


5.35 


2.78 


5.02 


2.30 


5.39 


2.83 


5.06 


2.36 


5.43 


2.87 


5.09 


2.42 


5.46 


2.92 


5.13 


2.48 


5.50 


2.96 



Table 8.44: Normalisation for a dipolar form for T(15') as a function of the b quark mass for 
A = 4.10 GeV. 



B.l Results relative to the normalisation 
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4.80 


1.24 


5.17 


1.64 


4.84 


1.31 


5.21 


1.67 


4.87 


1.3G 


5.24 


1.70 


4.91 


1.41 


5.28 


1.72 


4.95 


1.45 


5.32 


1.75 


4.98 


1.49 


5.35 


1.77 


5.02 


1.52 


5.39 


1.80 


5.06 


1.55 


5.43 


1.82 


5.09 


1.59 


5.46 


1.84 


5.13 


1.62 


5.50 


1.86 



Table B.46: Normalisation for a dipolar form for T(IS') as a function of the b quark mass for 
A = 5.20 GeV. 



4.80 


0.93 


5.17 


1.17 


4.84 


0.97 


5.21 


1.18 


4.87 


1.00 


5.24 


1.20 


4.91 


1.03 


5.28 


1.21 


4.95 


1.05 


5.32 


1.23 


4.98 


1.07 


5.35 


1.24 


5.02 


1.10 


5.39 


1.25 


5.06 


1.11 


5.43 


1.27 


5.09 


1.13 


5.46 


1.28 


5.13 


1.15 


5.50 


1.29 



Table B.48: Normalisation for a dipolar form for T(15') as a function of the b quark mass for 
A = 6.25 GeV. 



A 


JV 


A 


JV 


3.00 


2.53 


4.71 


1.31 


3.17 


2.33 


4.88 


1.24 


3.34 


2.16 


5.05 


1.18 


3.51 


2.00 


5.22 


1.12 


3.68 


1.87 


5.39 


1.07 


3.86 


1.75 


5.57 


1.02 


4.03 


1.64 


5.74 


0.98 


4.20 


1.54 


5.91 


0.94 


4.37 


1.46 


6.08 


0.90 


4.54 


1.38 


6.25 


0.87 



Table B.50: Normalisation for a gaussian form for T(15') as a function of A for rrif, = 4.80 GeV. 



3.00 


3.72 


4.71 


1.66 


3.17 


3.36 


4.88 


1.57 


3.34 


3.05 


5.05 


1.48 


3.51 


2.79 


5.22 


1.39 


3.68 


2.56 


5.39 


1.32 


3.86 


2.36 


5.57 


1.25 


4.03 


2.19 


5.74 


1.19 


4.20 


2.03 


5.91 


1.13 


4.37 


1.90 


6.08 


1.08 


4.54 


1.77 


6.25 


1.03 



Table B.52: Normalisation for a gaussian form for T(15') as a function of A for rrib = 5.00 GeV. 
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3.00 


4.60 


4.71 


1.91 


3.17 


4.11 


4.88 


1.79 


3.34 


3.70 


5.05 


1.68 


3.51 


3.35 


5.22 


1.58 


3.68 


3.05 


5.39 


1.49 


3.86 


2.79 


5.57 


1.40 


4.03 


2.57 


5.74 


1.33 


4.20 


2.37 


5.91 


1.26 


4.37 


2.20 


6.08 


1.20 


4.54 


2.05 


6.25 


1.14 



Table B.54: Normalisation for a gaussian form for T(IS') as a function of A for nib = 5.28 GeV. 



3.00 


5.37 


4.71 


2.11 


3.17 


4.76 


4.88 


1.97 


3.34 


4.26 


5.05 


1.84 


3.51 


3.83 


5.22 


1.73 


3.68 


3.47 


5.39 


1.62 


3.86 


3.16 


5.57 


1.53 


4.03 


2.89 


5.74 


1.44 


4.20 


2.66 


5.91 


1.37 


4.37 


2.46 


6.08 


1.29 


4.54 


2.27 


6.25 


1.23 



Table B.56: Normalisation for a gaussian form for T(15') as a function of A for nib = 5.50 GeV. 





N 




N 


4.80 


2.53 


5.17 


4.25 


4.84 


2.77 


5.21 


4.38 


4.87 


2.98 


5.24 


4.51 


4.91 


3.17 


5.28 


4.64 


4.95 


3.35 


5.32 


4.77 


4.98 


3.51 


5.35 


4.89 


5.02 


3.67 


5.39 


5.01 


5.06 


3.82 


5.43 


5.13 


5.09 


3.97 


5.46 


5.25 


5.13 


4.11 


5.50 


5.37 



Table B.58: Normalisation for a gaussian form for T(IS') as a function of the b quark mass for 
A = 3.00 GeV. 



4.80 


1.61 


5.17 


2.35 


4.84 


1.72 


5.21 


2.41 


4.87 


1.81 


5.24 


2.46 


4.91 


1.89 


5.28 


2.52 


4.95 


1.97 


5.32 


2.57 


4.98 


2.04 


5.35 


2.62 


5.02 


2.11 


5.39 


2.67 


5.06 


2.17 


5.43 


2.72 


5.09 


2.24 


5.46 


2.77 


5.13 


2.30 


5.50 


2.81 



Table 8.60: Normalisation for a gaussian form for T(15') as a function of the b quark mass for 
A = 4.10 GeV. 



B.2 Determination of parameters: another approach 
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4.80 


1.14 


5.17 


1.54 


4.84 


1.20 


5.21 


1.56 


4.87 


1.25 


5.24 


1.59 


4.91 


1.30 


5.28 


1.62 


4.95 


1.34 


5.32 


1.64 


4.98 


1.37 


5.35 


1.67 


5.02 


1.41 


5.39 


1.69 


5.06 


1.44 


5.43 


1.72 


5.09 


1.47 


5.4C 


1.74 


5.13 


1.51 


5.50 


1.77 



Table B.62: Normalisation for a gaussian form for T(IS') as a function of the b quark mass for 
A = 5.20 GeV. 



4.80 


0.87 


5.17 


1.10 


4.84 


0.90 


5.21 


1.11 


4.87 


0.93 


5.24 


1.13 


4.91 


0.96 


5.28 


1.14 


4.95 


0.98 


5.32 


1.16 


4.98 


1.00 


5.35 


1.17 


5.02 


1.02 


5.39 


1.19 


5.06 


1.04 


5.43 


1.20 


5.09 


1.06 


5.46 


1.21 


5.13 


1.08 


5.50 


1.23 



Table B.64: Normalisation for a gaussian form for T(IS') as a function of the b quark mass for 
A = 6.25 GeV. 



B.2 Determination of parameters: another approach 

Usually, the complete procedure to be followed is to obtain from the gap equation^ a relation 
between the quark mass and the bound-state mass; A - the paramater size- is determined 
from the leptonic decay width, whereas the normalisation is determined by the fact that 
the residue of the bound-state propagator be 1 at the bound-state mass. The value obtained 
for the quark mass from the gap equation depends unfortunately on A. This imposes the 
use of an iterative method to solve these coupled equations. 

In the following, we limit ourselves to the calculation of the leptonic decay width, to 
single out the few differences with our calculation but also to point out where there is a 
problem with the Wick rotation. 

We present here the calculation of the decay width for the dipole vertex function. The 
same procedure can be applied for the gaussian; instead of using Feynman-parameter method 
to gather the denominators, one has therefore to use Schwinger ones, or equally the properties 
of the Laplace transform. 



which is another way to write down the BSE. 
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B.2.1 Leptonic decay width 

From the Feynman rules and the vertex functions for the bound state, we have 



-l)e.3/^Trf(.rfch^).^-^----^-^ 



(2vr) 



m 



{k + y)" - m: 



-Sec 



d^k 2 .g^''{M'^ + ^m'^ -^k'^) + ^k^k'' -2P^'P'' 



(2vr) 



rr(p; 



rel) 



cq is the heavy quark charge, —1 comes for the fermionic loop, 3 is the colour factor. 
Since Prei = k, for the dipolar form of the vertex functions, we shall have 



r(p: 



rel) 



k^\2 
A2) 



(B.2) 



Bearing in mind that k is g^uk^k'^ and defining: 

d^k o . 1 



/ = 



d^k 



{{k - - ml){{k + - ml) 



■rip. 



ki'k'' 



(B.3) 



we have 



A^" = -2eiVA^ 



(g'^'iM^ + 4m') - SP^'P")! - Ag^"" g ^, P""' + 81"^ (B.4) 



Using the Feynman parametrisation for rational fractions, we have the following relation 

6x 



1 

—TT-. — ^= / dxdydz6(x + y + z — 1)- — — 
.42^2^3 Jo {xAi 



+ yA2 + zA^y 



(B.5) 



we may rewrite all the denominators as: 

1 



(F - A2)2((A; - 1)2 _ ml){{k + f )2 - m; 



— = y dxdydz5{x + y + z - I) — , (B.6) 



with 



D = x{k^ + A') - y{{k - - ml) + z{{k + - ml) + {y + z)ie 

M' 

= fc' + (z - y)k.P - xt^ + (?/ + ^)(^ m^) + (y + z)%e 



(B.7) 



B.2 Determination of parameters: another approach 
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The aim now is to get a complete square. To achieve this, we make use of the following 
change of variable 

i = k+ ^-^P ^k = i- 

2 2 

=^ fc2 ^ £2 ^ il^J^M^ -{z- y)LP (B.8) 

^ k.P = i.P-^^M'^, 
2 

which cancels linear terms in the variable. Indeed, we have 

(z - v)'^ 

D = f- L_^m2 - xA' + {y + z){—- m") + {y + z)ie ^^ ^^ 

= ^2 - A + (y + z)ie, 



with 



A = (2; - yY— + xA^ + {y + z) (m^ 



>o 



(B.IO) 



>0 if m?>^ 



Now we have. 



dxdydz5{..) 



d^k k^k''T{p'^^^i) A A A u \ [ '^^^ ^x^'T 



(27r)4 ((A; - f )2 _ nil){{k + f )2 -ml) Jo " ' ' J - A)^ 



dxdydz5{..) 



(B.ll) 



The previous lines implicitly require to make the Wick rotation one cannot trust as already 
explained. 

The integration over d'^l is carried out using these relations [102]: 



dH 1 I 2\ f iV i / 1 ^ ^ 



(27r)4 (£2 _ A)4 (47r)2 3!VAy 3(47r)2 ^A 
d^l -1 .^^'^ 1! / n -i g^" f I 



(27r)4 (£2 _ A)4 (47r)2 2 3! VA; 3(47r)2 4 
We are thus left with two types of integrals: 



(B.12) 



/■I X 

/ dxdydzS{..)— = A, 
Jo A 

X 

J dxdydzS{..)— = B, 



[B.13) 
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in terms of which A'^'^ reads 



Af"" = -2eiVA^ (3(^) ^^^''^^^^ + ^^^^ ~ 2P^P'^)5 



= [^^"((M'+4m')P + 2A) - 2P^P''B]. 

The integration on z is done using the 5-function, which gives z = 1 — x — y and the 
remaining domain of integration on {x,y) is then restricted byO<l— x — i/<0, that is 
< X < 1 and < y < 1 — x. 

First, let us concentrate on A. Defining 7^ = m? — we have 

'-^ ixdx f^~''' dy Axdx /"^"^ dy 



^ xdx /"^"^' (i?/ 



^I/J' Jo _ yl^l -x) + ^—J- + T72-^) 



M2 

(B.15) 



In order to simphfy the denominator, we make the following change of variable: 

1 — X 



y' = y 



2 ' 

1/0 = 0- (B.16) 

y^ = [I-X) — = — — . 



We thus obtain: 

1 — x -I 1 — x 

xdx dy f xdx dy 



;^ xcix M f (1- x)M 

A = 2 —— — arrtan 



(B.17) 



/o M2 ^72 + 3,(^2 _ ^2) y2v/72 + x(A2-72)_ 

Acting in the same way with B, we may also write 

xdx dy' 



/'^ T.li.T. /' 

B = 2 



^2^0 ^^,2 + £A!±gz£h!)^" ^^-^^^ 



which, integrating by part gives B = Bi + B2 with: 

^'^ 2x(l — x)dx 



(72 + X(A2 - 72)) (4(72 + x(A2 - 72)) + M2(l - X)2) 

1 xdx ( {l-x)M 



(B.19) 



B2 = — I arctan 



M Jo + a;(A2-72))l \2^-f^ + x{A^--f^) 



B.2 Determination of parameters: another approach 
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The latter three integrals are not computable analytically. We shall therefore use nu- 
merical method to obtain their value as function of A and M^. In our approach, we recall 
that we have choosen to take m as the mass of the lowest-lying meson composed of a heavy 
quark Q. 

On the other hand, we simply have 



(B.20) 



From the Feynman rules, we have 



Ue 



(B.21) 



where a and h are the momenta of the lepton and antilepton (a + 6 = P), and their mass. 
As is much smaller than M and m, we shall neglect it. 

To perform the integration on the two-particle phase space, we shall use the following 
well known relations [90]^ 



24 
24 



M2 



a2 ^2 (a2-62)2 

1 J I 2-^ — 

M2 M2 M4 



g^f^ + 2- 



pa p/3 



d2{PS)a.h = j d2{PS)a''bf^g 
Therefore, we have 

ttM^ 



M2 

«/3 = ^(4 + 2) 



fB.22) 



B 



{ley 



nM^g 



24 



pppp 



[le 



,pp 



pp pp 



M2 



(B.23) 



App' 



Now that all quantities in Eq. (4.20) are determined, we can put everything together, 
which gives'^: 



He?^M' 1 X 



X [gP'HM^ + Am^)B + 2A)]g,,APP' gp,,,[g'''^'' {{M^ + Am^)B + 2A)] 



187rM2 



(M^ + 4m2)P + 2AYA^^,AP''' 



((M2 + 4m2)P + 2A)2 
(B.24) 



2 A = A(l,aVM2,6VM2) = A{1, m^/M^, m^/Af^) ^ i _ 2mj/AP ~ 1. 
Recall the projector property of A^;^: A^^P'^ — A^^P^ — 0. 
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Eq. (4.19) then gives the width and provide with a relation between A and the masses. 

Combining this relation with the normalisation condition, it is therefore possible to fully 
determine the vertex-function parameters. The result obtained by Costa et al. [91] are given 
in Tab. (4.2). 
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